Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Computational fluid dynamics of the right ventricular outflow tract and of the pulmonary artery: a bench model of flow dynamics.

OBJECTIVES: The reconstruction of the right ventricular outflow tract (RVOT) with valved conduits remains a challenge. The reoperation rate at 5 years can be as high as 25% and depends on age, type of conduit, conduit diameter and principal heart malformation. The aim of this study is to provide a bench model with computer fluid dynamics to analyse the haemodynamics of the RVOT, pulmonary artery, its bifurcation, and left and right pulmonary arteries that in the future may serve as a tool for analysis and prediction of outcome following RVOT reconstruction. METHODS: Pressure, flow and diameter at the RVOT, pulmonary artery, bifurcation of the pulmonary artery, and left and right pulmonary arteries were measured in five normal pigs with a mean weight of 24.6 ± 0.89 kg. Data obtained were used for a 3D computer fluid-dynamics simulation of flow conditions, focusing on the pressure, flow and shear stress profile of the pulmonary trunk to the level of the left and right pulmonary arteries. RESULTS: Three inlet steady flow profiles were obtained at 0.2, 0.29 and 0.36 m/s that correspond to the flow rates of 1.5, 2.0 and 2.5 l/min flow at the RVOT. The flow velocity profile was constant at the RVOT down to the bifurcation and decreased at the left and right pulmonary arteries. In all three inlet velocity profiles, low sheer stress and low-velocity areas were detected along the left wall of the pulmonary artery, at the pulmonary artery bifurcation and at the ostia of both pulmonary arteries. CONCLUSIONS: This computed fluid real-time model provides us with a realistic picture of fluid dynamics in the pulmonary tract area. Deep shear stress areas correspond to a turbulent flow profile that is a predictive factor for the development of vessel wall arteriosclerosis. We believe that this bench model may be a useful tool for further evaluation of RVOT pathology following surgical reconstructions.

Type of prey influences biology and consumption rate of Orius insidiosus (Say) (Hemiptera, Anthocoridae)

Generalist predators are capable of consuming different types of prey, and as each prey may have distinct nutritional values, each may have a distinct impact on the biology of the predator. Therefore, the aim of this study was to determine how the consumption of different prey influences certain biological characteristics and the predatory capacity of Orius insidiosus (Say). The investigation was performed in climatic chamber at 25 ±1 ºC, RH 70 ± 10% and fotophase 12. Eggs of Anagasta kuehniella (Zeller), adults of Caliothrips phaseoli (Hood) and nymphs of Aphis gossypii Glover were used as prey and were provided daily ad libitum for all the mobile stages of the predator. The results showed that biological parameters of O. insidiosus are affected differently depending on the type of prey ingested. The development time of the nymphal stage was 13.1, 11.23 and 10.25 days for O. insidiosus feeding on eggs of A. kuehniella, nymphs of A. gossypii and adults of C. phaseoli, respectively. Longevity was five times larger for adults fed on eggs of A. kuehniella (56.25 days) compared to that of adults that preyed on nymphs of A. gossypii (11.44 days), and four times larger when the prey were adults of C. phaseoli (13.58 days). The consumption of eggs of A. kuehniella by predator females resulted in a shorter pre-oviposition period (3.2 days) and a longer oviposition period (44.4 days) when compared to the consumption of other types of prey. In addition, fecundity was increased with the consumption of eggs of A. kuehniella (195.25 eggs laid / female) when compared to feeding on the other prey, C. phaseoli (70.00 eggs laid / female) or A. gossypii (22.50 eggs laid / female). However, the consumption of aphids was larger (148.28 nymphs/ nymphal stage) than that of thrips (74.10 thrips / nymphal stage) or eggs of A. kuehniella (37.03 eggs /nymphal stage) for all of the nymphal stages of the predator. The results indicate that the eggs of A. kuehniella are the type of prey best suited for rearing this predator insect under laboratory conditions. Also fecundity was better with this prey even though the predator consumed during its nymphal stage a lesser quantity of eggs compared to the other prey.

Comparative biology of Cycloneda sanguinea (Linnaeus, 1763) and Hippodamia convergens Guérin-Méneville, 1842 (Coleoptera, Coccinellidae) focusing on the control of Cinara spp. (Hemiptera, Aphididae)

The giant conifer aphids Cinara pinivora (Wilson, 1919) and Cinara atlantica (Wilson, 1919) (Hemiptera, Aphididae) have been observed attacking Pinus spp. in Southern and Southeastern Brazil. The coccinellids, on the other hand, were found feeding on these aphids in the field, which can be regarded as potential biological control agents. The biological cycle and mortality rate of larvae of Cycloneda sanguinea (Linnaeus, 1763) and Hippodamia convergens Guérin-Méneville, 1842 (Coleoptera, Coccinellidae) were evaluated using twenty larvae of each predator species fed with nymphs of Cinara. The vials with the insects were kept under 15 ºC, 20 ºC and 25 ºC, with 12h photophase and 70 ± 10% relative humidity. The consumption was evaluated every 24 hours and the nymphs replaced. For C. sanguinea, the egg incubation time was 10.5, 5.0 and 4.0 days; the average larval development period was 33.3, 15.8 and 8.6 days and the larval mortality rate 20%,0% and 15%, respectively at 15 ºC, 20 ºC and 25 ºC. For H. convergens, the larval development time was 41.9, 19.3 and 10.9 days at 15 ºC, 20 ºC and 25 ºC, respectively. The larval mortality rate was 35%, 15% and 0% under the three temperatures. Both species developed adequately when fed nymphs of Cinara, however, C. sanguinea performed better than H. convergens, even at 15 ºC, at which temperature the biological cycles of the coccinellids are prolonged, but the temperature is favorable for the development of Cinara populations in the field.

Floral biology of Stachytarpheta maximiliani Scham. (Verbenaceae) and its floral visitors

This study describes the reproductive system of Stachytarpheta maximiliani (Verbenaceae), including its floral biology, nectar and pollen availability and insect foraging patterns, identifying whose species act as pollinators. It was carried out in a Brazilian Atlantic rain forest site. Observations on the pollination biology of the Verbenaceae S. maximiliani indicate that their flowering period extends from September through May. Anthesis occurs from 5:30 a.m. to 5:00 p.m. and nectar and pollen are available during all the anthesis. Many species of beetles, hemipterans, flies, wasps, bees and butterflies visit their flowers, but bees and butterflies are the most frequent visitors. The flowers are generally small, gathered in dense showy inflorescences. A complex of floral characteristcs, such as violet-blue color of flowers, long floral tubes, without scents, nectar not exposed, high concentration of sugar in nectar (about 32%), allowed identification of floral syndromes (melittophily and psicophily) and function for each visitor. The bees, Bombus morio, B. atratus, Trigonopedia ferruginea, Xylocopa brasilianorum and Apis mellifera and the butterflies Corticea mendica mendica, Corticea sp., Vehilius clavicula, Urbanus simplicius, U. teleus and Heraclides thoas brasiliensis, are the most important pollinators.

Nesting biology of the trap-nesting Neotropical wasp Trypoxylo n(Trypargilum) aurifrons Shuckard (Hymenoptera, Crabronidae)

The present study was carried out in three localities of the state of São Paulo, Brazil: Araras (Dec/03-Dec/06), São Carlos (Nov/04-Nov/06) and Rifaina (Jul/04-Dec/06). Trap-nests were distributed among sites in the sampling areas and were collected every 35 days. Data from 295 nests indicate that T. aurifrons is a multivoltine species, with higher rates of nest building and cell production in the warm, rainy season. The trap-nests used by the females ranged from 117 to 467 mm in length and 3.1 to 16.6 mm in diameter. All nests showed deep plugs and a vestibular cell was found in 37% of the complete nests. The number of cells per nest ranged from one to 12. Females were larger than males, emerged from longer cells and their cocoons were significantly larger. A secondary 1:1 sex ratio was found in Araras and Rifaina. No correlation was observed between the diameter of the trap-nest and sex ratio. Males were usually oviposited in the first brood cells. Male and female developmental time from egg to adult was longer in the cold, dry season. Trypoxylon aurifrons provisioned their nests mainly with orb-spiders from the family Araneidae. The most important mortality factor was the death of immature forms, probably due to development failure. The most important parasitoid was Melittobia sp.

Occurrence, biology and behavior of Liogenys fuscus Blanchard (Insecta, Coleoptera, Scarabaeidae) in Aquidauana, Mato Grosso do Sul, Brazil

Due to the importance of some Pleurosticti Scarabaeidae as agricultural pests allied to information absence on the species that occur in Brazilian Central-West region, on studies occurrence, biology and behavior on this group of scarabs were conducted. Biology and behavioral studies started with Liogenys fuscus Blanchard, 1850 (Melolonthinae), a very common species and were developed in Aquidauana, Mato Grosso do Sul. Adult beetles were collected from light traps from February 2005 to January 2007, at the experimental farm of the Universidade Estadual de Mato Grosso do Sul in Aquidauana (UEMS). In the laboratory adults were placed in plastic containers with soil with sprouts of Brachiaria decumbens Stapf (Poaceae). Eggs were transferred to a climatized chamber at 26 ± 1º C with a 12hourlight, 12hour darkness photoperiod cycle. Adult flight activity occurred in August and in September to December from 06:00 pm to 06:00 am, with the largest number of individuals flying from 07:00 to 10:00 pm. Eggs measured 1 x 1.5 mm and were laid individually or in groups in soil chambers; eggs were initially white and became yellow near hatching. The embryonic period lasted 14.3 days; first, second and third instars lasted 28.5, 48.8, and 68.2 days, respectively. The prepupal period lasted 120.2 days and the prepupa stayed inactive in soil. The mean duration of pupal stage was 27.5 days and the mean longevity of adults was 23.6 days. In laboratory the calling behavior between males and females was observed; copulation lasted, in mean, 25 minutes.

First remarks on the nesting biology of Hypodynerus andeus (Packard) (Hymenoptera, Vespidae, Eumeninae) in the Azapa valley, northern Chile

First remarks on the nesting biology of Hypodynerus andeus (Packard) (Hymenoptera, Vespidae, Eumeninae) in the Azapa valley, northern Chile. Some aspects about the nesting biology of the potter wasp Hypodynerus andeus (Packard, 1869) are reported for the first time. Observations were carried out at the Azapa valley, coastal desert of northern Chile. A total of sixty nests were collected and examined, each composed by 1-14 cells, most of them found attached to concrete lamp posts. The only preys recorded in the cells were Geometridae (Lepidoptera) caterpillars and the presence of the parasitoid Anthrax sp. (Diptera, Bombyliidae) was also recorded. A number of arthropods belonging to different groups, mainly spiders, were found occupying empty nests.

Developmental biology, polymorphism and ecological aspects of Stiretrus decemguttatus (Hemiptera, Pentatomidae), an important predator of cassidine beetles

Developmental biology, polymorphism and ecological aspects of Stiretrus decemguttatus (Hemiptera, Pentatomidae), an important predator of cassidine beetles. Stiretrus decemguttatus is an important predator of two species of cassidine beetles, Botanochara sedecimpustulata (Fabricius, 1781) and Zatrephina lineata (Fabricius, 1787) (Coleoptera, Cassidinae), on the Marajó Island, Brazil. It attacks individuals in all development stages, but preys preferentially on late-instar larvae. Its life cycle in the laboratory was 43.70 ± 1.09 days, with an egg incubation period of six days and duration from nymph and adult stages of 16.31 ± 0.11 and 22.10 ± 1.67 days, respectively. The duration of one generation (T) was 12.65 days and the intrinsic population growth rate (r) 0.25. These data reveal the adjustment of the life cycle of S. decemgutattus with those of the two preys, but suggest greater impact on Z. lineata. However, no preference over cassidine species was shown in the laboratory. Up to 17 different color patterns can be found in adults of S. decemguttatus, based on combinations of three basic sets of color markings. Some of them resemble the markings of chrysomelids associated with Ipomoea asarifolia (Convolvulaceae) and are possibly a mimetic ring. Three color patterns were identified in nymphs, none of which was associated with any specific adult color pattern.

On the nesting biology of Pirhosigma Giordani Soika (Hymenoptera, Vespidae, Eumeninae), with special reference to the use of vegetable matter

On the nesting biology of Pirhosigma Giordani Soika (Hymenoptera, Vespidae, Eumeninae), with special reference to the use of vegetable matter. The use of vegetable matter in nest building is not widespread among the Eumeninae, and is reported for the first time for the two species of potter wasps Pirhosigma superficiale and P. limpidum. These wasps make mostly spherical mud nests over which they attach small pieces of unmasticated plant matter. Use of plant fragments in this group of wasps is interpreted as camouflage behavior.