Neutral space analysis for a Boolean network model of the fission yeast cell cycle network

Background

Interactions between genes and their products give rise to complex circuits known as gene regulatory networks (GRN) that enable cells to process information and respond to external stimuli. Several important processes for life, depend of an accurate and context-specific regulation of gene expression, such as the cell cycle, which can be analyzed through its GRN, where deregulation can lead to cancer in animals or a directed regulation could be applied for biotechnological processes using yeast. An approach to study the robustness of GRN is through the neutral space. In this paper, we explore the neutral space of a Schizosaccharomyces pombe (fission yeast) cell cycle network through an evolution strategy to generate a neutral graph, composed of Boolean regulatory networks that share the same state sequences of the fission yeast cell cycle.

Results

Through simulations it was found that in the generated neutral graph, the functional networks that are not in the wildtype connected component have in general a Hamming distance more than 3 with the wildtype, and more than 10 between the other disconnected functional networks. Significant differences were found between the functional networks in the connected component of the wildtype network and the rest of the network, not only at a topological level, but also at the state space level, where significant differences in the distribution of the basin of attraction for the G₁ fixed point was found for deterministic updating schemes.

Conclusions

In general, functional networks in the wildtype network connected component, can mutate up to no more than 3 times, then they reach a point of no return where the networks leave the connected component of the wildtype. The proposed method to construct a neutral graph is general and can be used to explore the neutral space of other biologically interesting networks, and also formulate new biological hypotheses studying the functional networks in the wildtype network connected component.     

老年慢性心力衰竭患者肿瘤坏死因子、脑钠肽水平变化及其相关性分析*

摘要目的：探讨老年慢性心力衰竭（CHF）肿瘤坏死因子琢（TNF-α）以及脑利钠肽水平的变化,并对其进行相关性分析。方法：选 择老年CHF患者59 例,根据心功能分级,Ⅱ级19 例,Ⅲ级25 例,Ⅳ级15 例;另选心功能正常的老年人36例。所有入选者予以血 清BNP、TNF-α测定,并采用彩色多普勒超声测定左心室舒张末内径并计算左室射血分数（LVEF）。结果：与对照组比较,CHF组 心功能Ⅲ、Ⅳ级患者TNα-琢水平明显升高（P＜0.05）,心功能Ⅱ、Ⅲ、Ⅳ级患者血BNP水平明显升高（P＜0.05）,并随心功能分级加 重而增加。血TNF-α水平与BNP 呈正相关（r=0.57,P＜0.05）。CHF患者中血清TNF-α和BNP 水平与LVEF比较均呈负相关 （r=-0.48, r=-0.64, P＜0.05）。结论：老年CHF患者血TNF-α显著升高,并与BNP及LVEF关系密切,是反应CHF患者心功能恶化 的重要预测指标。     

A Framework to Optimize Biodiversity Restoration Efforts Based on Habitat Amount and Landscape Connectivity

The effectiveness of ecological restoration actions toward biodiversity conservation depends on both local and landscape constraints. Extensive information on local constraints is already available, but few studies consider the landscape context when planning restoration actions. We propose a multiscale framework based on the landscape attributes of habitat amount and connectivity to infer landscape resilience and to set priority areas for restoration. Landscapes with intermediate habitat amount and where connectivity remains sufficiently high to favor recolonization were considered to be intermediately resilient, with high possibilities of restoration effectiveness and thus were designated as priority areas for restoration actions. The proposed method consists of three steps: (1) quantifying habitat amount and connectivity; (2) using landscape ecology theory to identify intermediate resilience landscapes based on habitat amount, percolation theory, and landscape connectivity; and (3) ranking landscapes according to their importance as corridors or bottlenecks for biological flows on a broader scale, based on a graph theory approach. We present a case study for the Brazilian Atlantic Forest (approximately 150 million hectares) in order to demonstrate the proposed method. For the Atlantic Forest, landscapes that present high restoration effectiveness represent only 10% of the region, but contain approximately 15 million hectares that could be targeted for restoration actions (an area similar to today's remaining forest extent). The proposed method represents a practical way to both plan restoration actions and optimize biodiversity conservation efforts by focusing on landscapes that would result in greater conservation benefits .     

Aphid activities during sieve element punctures

Aphid salivation in sieve elements and phloem sap ingestion were linked to waveforms in the Electrical Penetration Graph (EPG). Non-viruliferousRhopalosiphum padi (L.) (Hemiptera, Aphididae) on barley yellow dwarf virus (BYDV) infected wheat could acquire the virus, which was used as an indication for phloem sap ingestion, whereas virus inoculation by viruliferous aphids on healthy plants was associated with salivation in sieve elements or other phloem cells. Probing was monitored and the waveforms recorded were related to ELISA results of test plants. The EPG patterns A, B, and C are indicative of the stylet pathway phase, whereas patterns E1 and E2 reflect the phloem (sieve element) phase with an unknown activity (E1) or with ingestion and concurrent salivation (E2). Aphids showing pathway and E1 rarely acquired virus, suggesting that little or no phloem sap ingestion can occur during these patterns, whereas those showing additionally pattern E2 did so substantially, indicating phloem sap ingestion. The main pattern related to virus inoculation was E1, although some aphids were able to inoculate plants during pathway. Pattern E1 clearly reflects the most important salivation into sieve elements. Pattern E2 had no clear contribution to virus inoculation, supporting the present hypothesis that during this pattern the saliva is mixed with the phloem sap in the single canal at the stylet tips and ingested immediately, without reaching the plant tissue. Sustained sap ingestion did not affect virus inoculation. So, BYDV inoculation mainly occurs during the first period of a sieve element puncture which is always formed by E1. Implications on persistent virus transmission are discussed.     

Evolutionary graph theory: breaking the symmetry between interaction and replacement

We study evolutionary dynamics in a population whose structure is given by two graphs: the interaction graph determines who plays with whom in an evolutionary game; the replacement graph specifies the geometry of evolutionary competition and updating. First, we calculate the fixation probabilities of frequency dependent selection between two strategies or phenotypes. We consider three different update mechanisms: birth-death, death-birth and imitation. Then, as a particular example, we explore the evolution of cooperation. Suppose the interaction graph is a regular graph of degree h, the replacement graph is a regular graph of degree g and the overlap between the two graphs is a regular graph of degree l. We show that cooperation is favored by natural selection if b/c>hg/l. Here, b and c denote the benefit and cost of the altruistic act. This result holds for death-birth updating, weak-selection and large population size. Note that the optimum population structure for cooperators is given by maximum overlap between the interaction and the replacement graph (g=h=l), which means that the two graphs are identical. We also prove that a modified replicator equation can describe how the expected values of the frequencies of an arbitrary number of strategies change on replacement and interaction graphs: the two graphs induce a transformation of the payoff matrix.     

Co-localisation of host plant resistance QTLs affecting the performance and feeding behaviour of the aphid Myzus persicae in the peach tree

The architecture and action of quantitative trait loci (QTL) contributing to plant resistance mechanisms against aphids, the largest group of phloem-feeding insects, are not well understood. Comparative mapping of several components of resistance to the green peach aphid (Myzus persicae) was undertaken in Prunus davidiana, a wild species related to peach. An interspecific F(1) population of Prunus persica var. Summergrand × P. davidiana clone P1908 was scored for resistance (aphid colony development and foliar damage) and 17 aphid feeding behaviour traits monitored by means of the electrical penetration graph technique. Seven resistance QTLs were detected, individually explaining 6.1-43.1% of the phenotypic variation. Consistency was shown over several trials. Nine QTLs affecting aphid feeding behaviour were identified. All resistance QTLs except one co-located with QTLs underlying aphid feeding behaviour. A P. davidiana resistance allele at the major QTL was associated with drastic reductions in phloem sap ingestion by aphids, suggesting a phloem-based resistance mechanism. Resistance was also positively correlated with aphid salivation into sieve elements, suggesting an insect response to restore the appropriate conditions for ingestion after phloem occlusion. No significant QTL was found for traits characterising aphid mouthpart activity in plant tissues other than phloem vessels. Two QTLs with effects on aphid feeding behaviour but without effect on resistance were identified. SSR markers linked to the main QTLs involved in resistance are of potential use in marker-assisted selection for aphid resistance. Linking our results with the recent sequencing of the peach genome may help clarify the physiological resistance mechanisms.     

Weak seed banks influence the signature and detectability of selective sweeps

Seed banking (or dormancy) is a widespread bet-hedging strategy, generating a form of population overlap, which decreases the magnitude of genetic drift. The methodological complexity of integrating this trait implies it is ignored when developing tools to detect selective sweeps. But, as dormancy lengthens the ancestral recombination graph (ARG), increasing times to fixation, it can change the genomic signatures of selection. To detect genes under positive selection in seed banking species it is important to (1) determine whether the efficacy of selection is affected, and (2) predict the patterns of nucleotide diversity at and around positively selected alleles. We present the first tree sequence-based simulation program integrating a weak seed bank to examine the dynamics and genomic footprints of beneficial alleles in a finite population. We find that seed banking does not affect the probability of fixation and confirm expectations of increased times to fixation. We also confirm earlier findings that, for strong selection, the times to fixation are not scaled by the inbreeding effective population size in the presence of seed banks, but are shorter than would be expected. As seed banking increases the effective recombination rate, footprints of sweeps appear narrower around the selected sites and due to the scaling of the ARG are detectable for longer periods of time. The developed simulation tool can be used to predict the footprints of selection and draw statistical inference of past evolutionary events in plants, invertebrates, or fungi with seed banks.     

Fine structure of aphid stylet routes in plant tissues in correlation with EPG signals

Abstract. Plant penetration by Aphis fabae (Scopoli) was recorded by the electrical penetration graph (EPG) technique and followed by stylectomy during wave-forms that were suspected of indicating sieve element punctures. The severed stylets in the plant tissue were subsequently processed for transmission electron microscopy (TEM) and sectioned either transverse or longitudinal to the stylets. Two completely serially sectioned probes from the epidermis to the phloem were reconstructed.
In one probe the stylet pathway went to a sieve element and showed many empty branches of salivary sheath material. Breaks in cell walls filled with sheath material demonstrated that the majority of cells bordering the track had been punctured, which supports earlier evidence from EPGs. All types of cells showed punctures and the highest number was found inside the vascular bundle. Very few cells died, which would appear to be important for virus transmission, and in others cellular reactions remained limited to some callose formation. The route of the stylets was intercellular and passed through the secondary wall material. The role of pectinase in intercellular penetration, and previous evidence for intracellular tracks are discussed. Most sieve elements had been punctured but only one was eventually accepted. Thus, reaching a sieve element in a host plant does not automatically imply its acceptance though the reason remains unclear. Gelation of phloem proteins was shown in the stylet canal.
In a second probe, plant cytological and morphological correlations with the EPG were emphasized. Probes by other aphid-plant combinations showed great similarity.     

Lethal and sublethal effect of heat shock on Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae)

Temperature is an important abiotic environmental factor, and is responsible for various kinds of behavioral and physiological changes in living organisms. Induced heat shock is associated with feeding behaviour, reproduction and reactive oxygen species (ROS) generation that causes oxidative damage. In this experiment, we examined the lethal and sublethal effects of heat shock on reproduction, feeding behaviour and antioxidant enzymes, including catalase (CAT), superoxide dismutase (SOD) and peroxidases (POD) in P. solenopsis. Results showed that males were highly susceptible to heat shock treatments than females, as LTemp₅₀ values were 43.8 °C for males and 45.11 °C for females. Heat shock events non-significantly affected the fecundity in female only treated adults and significantly affected the both sexes heat treated adults, it increased the xylem feeding duration, percentage of xylem feeding adults and reduce the phloem feeding duration and percentage of phloem feeding adults. Similarly it alter the antioxidant enzymes activities, an increase of CAT, SOD and POD activities were noticed in response to highest intensity of heat shock while a reduction of CAT and SOD activity were noticed in response to lowest intensity of heat shock compared to control (30 °C). These results suggest that heat shock may result in loss of body water and induce oxidative stress in P. solenopsis. However, antioxidant enzymes play a significant role in overcoming the oxidative damage.     

On adaptation: A reduction of the Kauffman-Levin model to a problem in graph theory and its consequences

It is shown that complex adaptations are best modelled as discrete processes represented on directed weighted graphs. Such a representation captures the idea that problems of adaptation in evolutionary biology are problems in a discrete space, something that the conventional representations using continuous adaptive landscapes does not. Further, this representation allows the utilization of well-known algorithms for the computation of several biologically interesting results such as the accessibility of one allele from another by a specified number of point mutations, the accessibility of alleles at a local maximum of fitness, the accessibility of the allele with the globally maximum fitness, etc. A reduction of a model due to Kauffman and Levin to such a representation is explicitly carried out and it is shown how this reduction clarifies the biological questions that are of interest.Thanks are due to William Wimsatt, James F. Crow, and the referees for Biology and Philosophy for comments on an earlier version of this paper. Remarks by members of the audience, especially Abner Shimony, of a seminar at Boston University, February 19, 1988, were also very helpful. The diagrams were prepared with the assistance of Tracy Lubas.