Heuristics for the Critical Node Detection Problem in Large Complex Networks

Complex networks have recently attracted a significant amount of research attention due to their ability to model real world phenomena. One important problem often encountered is to limit diffusive processes spread over the network, for example mitigating pandemic disease or computer virus spread. A number of problem formulations have been proposed that aim to solve such problems based on desired network characteristics, such as maintaining the largest network component after node removal. The recently formulated critical node detection problem aims to remove a small subset of vertices from the network such that the residual network has minimum pairwise connectivity. Unfortunately, the problem is NP-hard and also the number of constraints is cubic in number of vertices, making very large scale problems impossible to solve with traditional mathematical programming techniques. Even many approximation algorithm strategies such as dynamic programming, evolutionary algorithms, etc. all are unusable for networks that contain thousands to millions of vertices. A computationally efficient and simple approach is required in such circumstances, but none currently exist. In this thesis, such an algorithm is proposed. The methodology is based on a depth-first search traversal of the network, and a specially designed ranking function that considers information local to each vertex. Due to the variety of network structures, a number of characteristics must be taken into consideration and combined into a single rank that measures the utility of removing each vertex. Since removing a vertex in sequential fashion impacts the network structure, an efficient post-processing algorithm is also proposed to quickly re-rank vertices. Experiments on a range of common complex network models with varying number of vertices are considered, in addition to real world networks. The proposed algorithm, DFSH, is shown to be highly competitive and often outperforms existing strategies such as Google PageRank for minimizing pairwise connectivity.

Fast regridding of large, complex geospatial datasets

In the earth sciences, data are commonly cast on complex grids in order to model irregular domains such as coastlines, or to evenly distribute grid points over the globe. It is common for a scientist to wish to re-cast such data onto a grid that is more amenable to manipulation, visualization, or comparison with other data sources. The complexity of the grids presents a significant technical difficulty to the regridding process. In particular, the regridding of complex grids may suffer from severe performance issues, in the worst case scaling with the product of the sizes of the source and destination grids. We present a mechanism for the fast regridding of such datasets, based upon the construction of a spatial index that allows fast searching of the source grid. We discover that the most efficient spatial index under test (in terms of memory usage and query time) is a simple look-up table. A kd-tree implementation was found to be faster to build and to give similar query performance at the expense of a larger memory footprint. Using our approach, we demonstrate that regridding of complex data may proceed at speeds sufficient to permit regridding on-the-fly in an interactive visualization application, or in a Web Map Service implementation. For large datasets with complex grids the new mechanism is shown to significantly outperform algorithms used in many scientific visualization packages.

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

Method to assess and optimise dependability of complex macro-systems: application to a railway signalling system

Achieving a high degree of dependability in complex macro-systems is challenging. Because of the large number of components and numerous independent teams involved, an overview of the global system performance is usually lacking to support both design and operation adequately. A functional failure mode, effects and criticality analysis (FMECA) approach is proposed to address the dependability optimisation of large and complex systems. The basic inductive model FMECA has been enriched to include considerations such as operational procedures, alarm systems. environmental and human factors, as well as operation in degraded mode. Its implementation on a commercial software tool allows an active linking between the functional layers of the system and facilitates data processing and retrieval, which enables to contribute actively to the system optimisation. The proposed methodology has been applied to optimise dependability in a railway signalling system. Signalling systems are typical example of large complex systems made of multiple hierarchical layers. The proposed approach appears appropriate to assess the global risk- and availability-level of the system as well as to identify its vulnerabilities. This enriched-FMECA approach enables to overcome some of the limitations and pitfalls previously reported with classical FMECA approaches.

Analysis and validation of CFD wind farm models in complex terrain. Effects induced by topography and wind turbines

Wind farms have been extensively simulated through engineering models for the estimation of wind speed and power deficits inside wind farms. These models were designed initially for a few wind turbines located in flat terrain. Other models based on the parabolic approximation of Navier Stokes equations were developed, making more realistic and feasible the operational resolution of big wind farms in flat terrain and offshore sites. These models have demonstrated to be accurate enough when solving wake effects for this type of environments. Nevertheless, few analyses exist on how complex terrain can affect the behaviour of wind farm wake flow. Recent numerical studies have demonstrated that topographical wakes induce a significant effect on wind turbines wakes, compared to that on flat terrain. This circumstance has recommended the development of elliptic CFD models which allow global simulation of wind turbine wakes in complex terrain. An accurate simplification for the analysis of wind turbine wakes is the actuator disk technique. Coupling this technique with CFD wind models enables the estimation of wind farm wakes preserving the extraction of axial momentum present inside wind farms. This paper describes the analysis and validation of the elliptical wake model CFDWake 1.0 against experimental data from an operating wind farm located in complex terrain. The analysis also reports whether it is possible or not to superimpose linearly the effect of terrain and wind turbine wakes. It also represents one of the first attempts to observe the performance of engineering models compares in large complex terrain wind farms.

Low density lipoprotein receptor-negative mice expressing human apolipoprotein B-100 develop complex atherosclerotic lesions on a chow diet: No accentuation by apolipoprotein(a)

We have generated mice with markedly elevated plasma levels of human low density lipoprotein (LDL) and reduced plasma levels of high density lipoprotein. These mice have no functional LDL receptors [LDLR−/−] and express a human apolipoprotein B-100 (apoB) transgene [Tg(apoB+/+)] with or without an apo(a) transgene [Tg(apoa+/−)]. Twenty animals (10 males and 10 females) of each of the following four genotypes were maintained on a chow diet: (i) LDLR−/−, (ii) LDLR−/−;Tg(apoa+/−), (iii) LDLR−/−;Tg(apoB+/+), and (iv)LDLR−/−;Tg(apoB+/+);Tg(apo+/−). The mice were killed at 6 mo, and the percent area of the aortic intimal surface that stained positive for neutral lipid was quantified. Mean percent areas of lipid staining were not significantly different between the LDLR−/− and LDLR−/−;Tg(apoa+/−) mice (1.0 ± 0.2% vs. 1.4 ± 0.3%). However, the LDLR−/−;Tg(apoB+/+) mice had ≈15-fold greater mean lesion area than the LDLR−/− mice. No significant difference was found in percent lesion area in the LDLR−/−;Tg(apoB+/+) mice whether or not they expressed apo(a) [18.5 ± 2.5%, without lipoprotein(a), Lp(a), vs. 16.0 ± 1.7%, with Lp(a)]. Histochemical analyses of the sections from the proximal aorta of LDLR−/−;Tg(apoB+/+) mice revealed large, complex, lipid-laden atherosclerotic lesions that stained intensely with human apoB-100 antibodies. In mice expressing Lp(a), large amounts of apo(a) protein colocalized with apoB-100 in the lesions. We conclude that LDLR−/−; Tg(apoB+/+) mice exhibit accelerated atherosclerosis on a chow diet and thus provide an excellent animal model in which to study atherosclerosis. We found no evidence that apo(a) increased atherosclerosis in this animal model.

Complex Patterns of Alternative Splicing Mediate the Spatial and Temporal Distribution of Perlecan/UNC-52 in Caenorhabditis elegans

The unc-52 gene encodes the nematode homologue of mammalian perlecan, the major heparan sulfate proteoglycan of the extracellular matrix. This is a large complex protein with regions similar to low-density lipoprotein receptors, laminin, and neural cell adhesion molecules (NCAMs). In this study, we extend our earlier work and demonstrate that a number of complex isoforms of this protein are expressed through alternative splicing. We identified three major classes of perlecan isoforms: a short form lacking the NCAM region and the C-terminal agrin-like region; a medium form containing the NCAM region, but still lacking the agrin-like region; and a newly identified long form that contains all five domains present in mammalian perlecan.  Using region-specific antibodies and unc-52 mutants, we reveal a complex spatial and temporal expression pattern for these UNC-52 isoforms. As well, using a series of mutations affecting different regions and thus different isoforms of UNC-52, we demonstrate that the medium NCAM-containing isoforms are sufficient for myofilament lattice assembly in developing nematode body-wall muscle. Neither short isoforms nor isoforms containing the C-terminal agrin-like region are essential for sarcomere assembly or muscle cell attachment, and their role in development remains unclear.

The Cellular Level of PR500, a Protein Complex Related to the 19S Regulatory Particle of the Proteasome, Is Regulated in Response to Stresses in Plants

In Arabidopsis seedlings and cauliflower florets, Rpn6 (a proteasome non-ATPase regulatory subunit) was found in two distinct protein complexes of ∼800 and 500 kDa, respectively. The large complex likely represents the proteasome 19S regulator particle (RP) because it displays the expected subunit composition and all characteristics. The small complex, designated PR500, shares at least three subunits with the “lid” subcomplex of 19S RP and is loosely associated with an hsp70 protein. In Arabidopsis COP9 signalosome mutants, PR500 was specifically absent or reduced to an extent that correlates with the severity of the mutations. Furthermore, PR500 was also diminished in response to potential protein-misfolding stresses caused by the heat shock and canavanine treatment. Immunofluorescence studies suggest that PR500 has a distinct localization pattern and is enriched in specific nuclear foci. We propose that PR500 may be evolved in higher plants to cope with the frequently encountered environmental stresses.

Two new genera, Provitellus and Ovipusillus, and four new species of Monorchiidae (Digenea) from carangid fishes of queensland, Australia

Two new genera and four new species of monorchiid digeneans are described from the Great Barrier Reef and Moreton Bay, Queensland. Provitellus turrum n. g., n. sp. from Pseudocaranx dentex and Trachinotus coppingeri is characterised by the presence of vitelline follicles in the forebody, a single testis, a unipartite terminal organ and filamented eggs. Ovipusillus mayu n. g., n. sp. from Gnathanodon speciosus is characterised by the presence of two testes, vitelline follicles overlapping the ventral sucker and a large, complex cirrus-sac that contains a coiled eversible ejaculatory duct joined by the pars prostatica halfway along its length. Paramonorcheides pseudocaranxi n. sp. from Pseudocaranx dentex differs from other species described in this genus in the longer flatter forebody, entire ovary and the well-developed cirrus-sac. Chrisomon gaigai n. sp. from Trachinotus coppingeri and T botla is characterised by the unflattened forebody and transversely oval pharynx. Chrisomon is redefined to include species of Lasiotocus with a vitellarium composed of clusters of tubular acini, creating the following new combinations: C. albulae n. comb. for L. albulae Overstreet, 1969, C. ulua n. comb, for L. ulua Yamaguti, 1970 and C. weke n. comb, for L. weke Yamaguti, 1970. The diagnosis of Lasiotocus is amended accordingly and the new combinations, L. polynemi n. comb. and L. sunderbanensis n. comb., are created for C.polynemi Dutta, Hafeezullah & Manna, 1994 and C. sunderbanensis Dutta, Hafeezullah B Manna, 1994, respectively. Extrapolation of our collection data suggests that there may be as many as 80 species of monorchiids infecting carangid fishes in Australia and 180 species infecting carangids in all oceans of the world. The latter figure greatly exceeds the number of monorchiids described from all host families to date.

Clinical and molecular genetics of myoclonic-astatic epilepsy and severe myoclonic epilepsy in infancy (Dravet syndrome)

The majority of severe epileptic encephalopathies of early childhood are symptomatic where a clear etiology is apparent. There is a small subgroup, however, where no etiology is found on imaging and metabolic studies, and genetic factors are important. Myoclonic-astatic epilepsy (MAE) and severe myoclonic epilepsy in infancy (SMEI), also known as Dravet syndrome, are epileptic encephalopathies where multiple seizure types begin in the first few years of life associated with developmental slowing. Clinical and molecular genetic studies of the families of probands with MAE and SMEI suggest a genetic basis. MAE was originally identified as part of the genetic epilepsy syndrome generalized epilepsy with febrile seizures plus (GEFS(+)). Recent clinical genetic studies suggest that SMEI forms the most severe end of the spectrum of the GEFS(+). GEF(+) has now been associated with molecular defects in three sodium channel subunit genes and a GABA subunit gene. Molecular defects of these genes have been identified in patients with MAE and SMEI. Interestingly, the molecular defects in MAE have been found in the setting of large GEFS(+) pedigrees, whereas, more severe truncation mutations arising de novo have been identified in patients with SMEI. It is likely that future molecular studies will shed light on the interaction of a number of genes, possibly related to the same or different ion channels, which result in a severe phenotype such as MAE and SMEI. (C) 2001 Elsevier Science B.V. All rights reserved.

Deletions within COL11A1 in Type 2 Stickler Syndrome Detected by Multiplex Ligation - Syndrome Detected by Multiplex Ligation Dependent Probe Amplification (MLPA)

Background: COL11A1 is a large complex gene around 250 kb in length and consisting of 68 exons. Pathogenic mutations in the gene can result in Stickler syndrome, Marshall syndrome or Fibrochondrogenesis. Many of the mutations resulting in either Stickler or Marshall syndrome alter splice sites and result in exon skipping, which because of the exon structure of collagen genes usually leaves the message in-frame. The mutant protein then exerts a dominant negative effect as it co-assembles with other collagen gene products. To date only one large deletion of 40 kb in the COL11A1, which was detected by RT-PCR, has been characterized. However, commonly used screening protocols, utilizing genomic amplification and exon sequencing, are unlikely to detect such large deletions. Consequently the frequency of this type of mutation is unknown. Case presentations: We have used Multiplex Ligation-Dependent Probe Amplification (MLPA) in conjunction with exon amplification and sequencing, to analyze patients with clinical features of Stickler syndrome, and have detected six novel deletions that were not found by exon sequencing alone. Conclusion: Exon deletions appear to represent a significant proportion of type 2 Stickler syndrome. This observation was previously unknown and so diagnostic screening of COL11A1 should include assays capable of detecting both large and small deletions, in addition to exon sequencing.

Evaluation of the determinism of the actuator sensor interface (ASI) on programmable logic controllers (PLC), personal computers (PC) and real-time operating systems (RTOS)

The purpose of this study was to evaluate the determinism of the AS-lnterface network and the 3 main families of control systems, which may use it, namely PLC, PC and RTOS. During the course of this study the PROFIBUS and Ethernet field level networks were also considered in order to ensure that they would not introduce unacceptable latencies into the overall control system. This research demonstrated that an incorrectly configured Ethernet network introduces unacceptable variable duration latencies into the control system, thus care must be exercised if the determinism of a control system is not to be compromised. This study introduces a new concept of using statistics and process capability metrics in the form of CPk values, to specify how suitable a control system is for a given control task. The PLC systems, which were tested, demonstrated extremely deterministic responses, but when a large number of iterations were introduced in the user program, the mean control system latency was much too great for an AS-I network. Thus the PLC was found to be unsuitable for an AS-I network if a large, complex user program Is required. The PC systems, which were tested were non-deterministic and had latencies of variable duration. These latencies became extremely exaggerated when a graphing ActiveX was included in the control application. These PC systems also exhibited a non-normal frequency distribution of control system latencies, and as such are unsuitable for implementation with an AS-I network. The RTOS system, which was tested, overcame the problems identified with the PLC systems and produced an extremely deterministic response, even when a large number of iterations were introduced in the user program. The RTOS system, which was tested, is capable of providing a suitable deterministic control system response, even when an extremely large, complex user program is required.

Protein-binding site at the immunoglobulin mu membrane polyadenylylation signal: possible role in transcription termination.

mRNAs specifying immunoglobulin mu and delta heavy chains are encoded by a single large, complex transcription unit (mu + delta gene). The transcriptional activity of delta gene segments in terminally differentiated, IgM-secreting B lymphocytes is 10-20 times lower than in earlier B-lineage cells expressing delta mRNA. We find that transcription of the mu + delta gene in IgM-secreting murine myeloma cells terminates within a region of 500-1000 nucleotides immediately following the mu membrane (mu m) polyadenylylation site. Transcription decreases only minimally through this region in murine cell lines representative of earlier stages in B-cell development. A DNA fragment containing the mu m polyadenylylation signal gives protein-DNA complexes with different mobilities in gel retardation assays with nuclear extracts from myeloma cells than with nuclear extracts from earlier B-lineage cells. However, using a recently developed "footprinting" procedure in which protein-DNA complexes resolved in gel retardation assays are subjected to nucleolytic cleavage while still in the polyacrylamide gel, we find that the DNA sequences protected by factors from the two cell types are indistinguishable. The factor-binding site on the DNA is located 5' of the mu m polyadenylylation signal AATAAA and includes the 15-nucleotide-long A + T-rich palindrome CTGTAAACAAATGTC. This type of palindromic binding site exhibits orientation-dependent activity consistent with the reported properties of polymerase II termination signals. This binding site is followed by two sets of directly repeated DNA sequences with different helical conformation as revealed by their reactivity with the chemical nuclease 1,10-phenanthroline-copper. The close proximity of these features to the signals for mu m mRNA processing may reflect a linkage of the processes of developmentally regulated mu m polyadenylylation and transcription termination.

Information systems for water resources management improvement in developing countries

Water resources management, as also water service provision projects in developing countries have difficulties to take adequate decisions due to scarce reliable information, and a lack of proper information managing. Some appropriate tools need to be developed in order to improve decision making to improve water management and access of the poorest, through the design of Decision Support Systems (DSS). On the one side, a DSS for developing co-operation projects on water access improvement has been developed. Such a tool has specific context constrains (structure of the system, software requirements) and needs (Logical Framework Approach monitoring, organizational-learning, accountability and evaluation) that shall be considered for its design. Key aspects for its successful implementation have appeared to be a participatory design of the system and support of the managerial positions at the inception phase. A case study in Tanzania was conducted, together with the Spanish NGO ONGAWA – Ingeniería para el Desarrollo. On the other side, DSS are required also to improve decision making on water management resources in order to achieve a sustainable development that not only improves the living conditions of the population in developing countries, but that also does not hinder opportunities of the poorest on those context. A DSS made to fulfil these requirements shall be using information from water resources modelling, as also on the environment and the social context. Through the research, a case study has been conducted in the Central Rift Valley of Ethiopia, an endhorreic basin 160 km south of Addis Ababa. There, water has been modelled using ArcSWAT, a physically based model which can assess the impact of land management practices on large complex watersheds with varying soils, land use and management conditions over long periods of time. Moreover, governance on water and environment as also the socioeconomic context have been studied.

Anàlisi del sistema productiu i comercial del vi de la Tarraconensis a partir de la combinació de metodologies

Des del segon quart del s. I aC i, especialment, durant el regnat d’August, es va desenvolupar a l’antiga província Tarraconensis un sistema productiu centrat en l’explotació agrària vitivinícola amb una finalitat clarament comercial. La majoria d’assentament vitivinícoles es troben emplaçats al litoral català, associats de vegades a figlinae que fabricaven les àmfores per al transport i comerç de l’excedent vinícola. No obstant, a l’àrea del Vallès Occidental i del Baix Llobregat es troben una sèrie de vil•les vinculades a la producció de vi i a la fabricació d’àmfores que han proporcionat restes molt significatives sobre la contribució d’aquesta zona a l’expansió econòmica de la província. La caracterització arqueològica i arqueomètrica d’un gran nombre d’àmfores procedents de diversos tallers ceràmics situats al Vallès Occidental i al Baix Llobregat, utilitzant diverses tècniques d’anàlisi química, mineralògica i petrogràfica, ha portat a establir quins tipus d’àmfores es van fabricar a cada taller i de quina manera. S’han identificat alguns dels processos tecnològics de la cadena operativa: la selecció i processat de les matèries primeres per conformar la pasta procedents, generalment, de l’àrea on es troba cada centre de producció, el modelatge, l’assecat i la cocció de les peces. En alguns dels casos analitzats, s’ha identificat quins tipus de contenidors van ser importants a l’establiment i la seva provinença. La integració d’aquests resultats en la base de dades analítica que disposa l’ERAAUB ha permès avaluar el grau d’estandardització dels processos tecnològics en aquesta àrea. La contrastació final amb les dades històriques i arqueològiques contribueix al coneixement arqueològic de les àmfores vinàries de la Tarraconensis i, a través d’elles, al coneixement de les societats que les van fabricar, comercialitzar i utilitzar.