Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

Discriminating Different Classes of Biological Networks by Analyzing the Graphs Spectra Distribution

The brain's structural and functional systems, protein-protein interaction, and gene networks are examples of biological systems that share some features of complex networks, such as highly connected nodes, modularity, and small-world topology. Recent studies indicate that some pathologies present topological network alterations relative to norms seen in the general population. Therefore, methods to discriminate the processes that generate the different classes of networks (e. g., normal and disease) might be crucial for the diagnosis, prognosis, and treatment of the disease. It is known that several topological properties of a network (graph) can be described by the distribution of the spectrum of its adjacency matrix. Moreover, large networks generated by the same random process have the same spectrum distribution, allowing us to use it as a "fingerprint". Based on this relationship, we introduce and propose the entropy of a graph spectrum to measure the "uncertainty" of a random graph and the Kullback-Leibler and Jensen-Shannon divergences between graph spectra to compare networks. We also introduce general methods for model selection and network model parameter estimation, as well as a statistical procedure to test the nullity of divergence between two classes of complex networks. Finally, we demonstrate the usefulness of the proposed methods by applying them to (1) protein-protein interaction networks of different species and (2) on networks derived from children diagnosed with Attention Deficit Hyperactivity Disorder (ADHD) and typically developing children. We conclude that scale-free networks best describe all the protein-protein interactions. Also, we show that our proposed measures succeeded in the identification of topological changes in the network while other commonly used measures (number of edges, clustering coefficient, average path length) failed.

Development and Evaluation of a Pseudoreference Pt//Ag/AgCl Electrode for Electrochemical Systems

The use of standard reference electrodes, such as Ag/AgCl or saturated calomel electrodes, in potentiometric and amperometric studies involving miniaturized electrochemical systems, or those operating under positive hydraulic pressure, is often impractical. Placement of the reference electrode in the direct vicinity of the working electrode is often prohibited by the dimensions or layout of the electrochemical cell, while the alternative strategy of locating the reference electrode in a separate compartment often leads to electrolyte leakage and contamination of the system. In the present study, we have investigated the functionality of a pseudoreference electrode comprising a platinum wire, one end of which was maintained in intimate contact with the internal solution of an Ag/AgCl reference electrode while the other was connected, via a BNC connector, to a platinum probe located within the electrochemical cell. Linear and cyclic voltammetric studies, involving both aqueous and nonaqueous electrolytes, were carried out using the pseudoreference electrode and an electrochemical cup-type cell with three electrodes or an electrochemical flow reactor. In all cases, the functionality of the Pt//Ag/AgCl system was similar to that of a conventional Ag/AgCl reference electrode. Variations in the electrolyte did not alter the potential or voltammetric profile recorded when using the pseudoreference system, although peak currents were generally improved and potential values shifted by approximately +350 mV in comparison with the Ag/AgCl electrode, therefore, the system pseudoreference can be applied in any electrochemical system due to the constant potential difference. It is concluded that the pseudoreference electrode can be used with advantage to obtain potentiometric and amperometric measurements in both simple and complex electrochemical systems.

Genetic and environmental risk factors in congenital heart disease functionally converge in protein networks driving heart development

Congenital heart disease (CHD) occurs in similar to 1% of newborns. CHD arises from many distinct etiologies, ranging from genetic or genomic variation to exposure to teratogens, which elicit diverse cell and molecular responses during cardiac development. To systematically explore the relationships between CHD risk factors and responses, we compiled and integrated comprehensive datasets from studies of CHD in humans and model organisms. We examined two alternative models of potential functional relationships between genes in these datasets: direct convergence, in which CHD risk factors significantly and directly impact the same genes and molecules and functional convergence, in which risk factors significantly impact different molecules that participate in a discrete heart development network. We observed no evidence for direct convergence. In contrast, we show that CHD risk factors functionally converge in protein networks driving the development of specific anatomical structures (e.g., outflow tract, ventricular septum, and atrial septum) that are malformed by CHD. This integrative analysis of CHD risk factors and responses suggests a complex pattern of functional interactions between genomic variation and environmental exposures that modulate critical biological systems during heart development.

Medium Light and Medium Roast Paper-Filtered Coffee Increased Antioxidant Capacity in Healthy Volunteers: Results of a Randomized Trial

We compared the effects of medium light roast (MLR) and medium roast (MR) paper-filtered coffee on antioxidant capacity and lipid peroxidation in healthy volunteers. In a randomized crossover study, 20 volunteers consumed 482 +/- 61 ml/day of MLR or MR for four weeks. Plasma total antioxidant status (TAS), oxygen radical absorbance capacity (ORAC), oxidized LDL and 8-epi-prostaglandin F2 alpha, erythrocyte superoxide dismutase (SOD), glutathione peroxidase (GPx), and catalase (CAT) activity were measured at baseline and after the interventions. MLR had higher chlorogenic acids-(CGA; 334 mg/150 mL) and less caffeine (231 mg/150 ml) than MR had (210 and 244 mg/150 ml, respectively). MLR also had fewer Maillard reaction products (MRP) than MR had. Compared with baseline, subjects had an increase of 21 and 26 % in TAS, 13 and 13 % in CAT, 52 and 75 % in SOD, and 62 and 49 % in GPx after MLR and MR consumption (P < 0.001), respectively. ORAC increased after MLR (P = 0.004). No significant alteration in lipid peroxidation biomarkers was observed. Both coffees had antioxidant effects. Although MLR contained more CGA, there were similar antioxidant effects between the treatments. MRP may have contributed as an antioxidant. These effects may be important in protecting biological systems and reducing the risk of diseases related to oxidative stress.

Genetic analyses of smoking initiation, persistence, quantity, and age-at-onset of regular cigarette use in Brazilian families: the Baependi Heart Study

Background: The purpose of this study was to estimate the genetic influences on the initiation of cigarette smoking, the persistence, quantity and age-at-onset of regular cigarette use in Brazilian families. Methods: The data set consisted of 1,694 individuals enrolled in the Baependi Heart Study. The heritability and the heterogeneity in genetic and environmental variance components by gender were estimated from variance components approaches, using the SOLAR (Sequential Oligogenic Linkage Analysis Routines) computer package. The mixed-effects Cox model was used for the genetic analysis of the age-at onset of regular cigarette use. Results: The heritability estimates were high (> 50%) for smoking initiation and were intermediate, ranging from 23.4 to 31.9%, for smoking persistence and quantity. Significant evidence for heterogeneity in variance components by gender was observed for smoking initiation and age-at-onset of regular cigarette use. Genetic factors play an important role in the interindividual variation of these phenotypes in females, while in males there is a predominant environmental component, which could be explained by greater social influences in the initiation of tobacco use. Conclusions: Significant heritabilities were observed in smoking phenotypes for both males and females from the Brazilian population. These data add to the literature and are concordant with the notion of significant biological determination in smoking behavior. Samples from the Baependi Heart Study may be valuable for the mapping of genetic loci that modulate this complex biological trait.

MODULARITY, NOISE, AND NATURAL SELECTION

Most biological systems are formed by component parts that are to some degree interrelated. Groups of parts that are more associated among themselves and are relatively autonomous from others are called modules. One of the consequences of modularity is that biological systems usually present an unequal distribution of the genetic variation among traits. Estimating the covariance matrix that describes these systems is a difficult problem due to a number of factors such as poor sample sizes and measurement errors. We show that this problem will be exacerbated whenever matrix inversion is required, as in directional selection reconstruction analysis. We explore the consequences of varying degrees of modularity and signal-to-noise ratio on selection reconstruction. We then present and test the efficiency of available methods for controlling noise in matrix estimates. In our simulations, controlling matrices for noise vastly improves the reconstruction of selection gradients. We also perform an analysis of selection gradients reconstruction over a New World Monkeys skull database to illustrate the impact of noise on such analyses. Noise-controlled estimates render far more plausible interpretations that are in full agreement with previous results.

Determination of the critical coupling of explosive synchronization transitions in scale-free networks by mean-field approximations

An explosive synchronization can be observed in scale-free networks when Kuramoto oscillators have natural frequencies equal to their number of connections. The present paper reports on mean-field approximations to determine the critical coupling of such explosive synchronization. It has been verified that the equation obtained for the critical coupling has an inverse dependence on the network average degree. This expression differs from those whose frequency distributions are unimodal and even. In this case, the critical coupling depends on the ratio between the first and second statistical moments of the degree distribution. Numerical simulations were also conducted to verify our analytical results.

A graph-theoretical approach in brain functional networks. Possible implications in EEG studies

Abstract Background Recently, it was realized that the functional connectivity networks estimated from actual brain-imaging technologies (MEG, fMRI and EEG) can be analyzed by means of the graph theory, that is a mathematical representation of a network, which is essentially reduced to nodes and connections between them. Methods We used high-resolution EEG technology to enhance the poor spatial information of the EEG activity on the scalp and it gives a measure of the electrical activity on the cortical surface. Afterwards, we used the Directed Transfer Function (DTF) that is a multivariate spectral measure for the estimation of the directional influences between any given pair of channels in a multivariate dataset. Finally, a graph theoretical approach was used to model the brain networks as graphs. These methods were used to analyze the structure of cortical connectivity during the attempt to move a paralyzed limb in a group (N=5) of spinal cord injured patients and during the movement execution in a group (N=5) of healthy subjects. Results Analysis performed on the cortical networks estimated from the group of normal and SCI patients revealed that both groups present few nodes with a high out-degree value (i.e. outgoing links). This property is valid in the networks estimated for all the frequency bands investigated. In particular, cingulate motor areas (CMAs) ROIs act as ‘‘hubs’’ for the outflow of information in both groups, SCI and healthy. Results also suggest that spinal cord injuries affect the functional architecture of the cortical network sub-serving the volition of motor acts mainly in its local feature property. In particular, a higher local efficiency El can be observed in the SCI patients for three frequency bands, theta (3-6 Hz), alpha (7-12 Hz) and beta (13-29 Hz). By taking into account all the possible pathways between different ROI couples, we were able to separate clearly the network properties of the SCI group from the CTRL group. In particular, we report a sort of compensatory mechanism in the SCI patients for the Theta (3-6 Hz) frequency band, indicating a higher level of “activation” Ω within the cortical network during the motor task. The activation index is directly related to diffusion, a type of dynamics that underlies several biological systems including possible spreading of neuronal activation across several cortical regions. Conclusions The present study aims at demonstrating the possible applications of graph theoretical approaches in the analyses of brain functional connectivity from EEG signals. In particular, the methodological aspects of the i) cortical activity from scalp EEG signals, ii) functional connectivity estimations iii) graph theoretical indexes are emphasized in the present paper to show their impact in a real application.

Preimplantation development regulatory pathway construction through a text-mining approach

Background: The integration of sequencing and gene interaction data and subsequent generation of pathways and networks contained in databases such as KEGG Pathway is essential for the comprehension of complex biological processes. We noticed the absence of a chart or pathway describing the well-studied preimplantation development stages; furthermore, not all genes involved in the process have entries in KEGG Orthology, important information for knowledge application with relation to other organisms. Results: In this work we sought to develop the regulatory pathway for the preimplantation development stage using text-mining tools such as Medline Ranker and PESCADOR to reveal biointeractions among the genes involved in this process. The genes present in the resulting pathway were also used as seeds for software developed by our group called SeedServer to create clusters of homologous genes. These homologues allowed the determination of the last common ancestor for each gene and revealed that the preimplantation development pathway consists of a conserved ancient core of genes with the addition of modern elements. Conclusions: The generation of regulatory pathways through text-mining tools allows the integration of data generated by several studies for a more complete visualization of complex biological processes. Using the genes in this pathway as “seeds” for the generation of clusters of homologues, the pathway can be visualized for other organisms. The clustering of homologous genes together with determination of the ancestry leads to a better understanding of the evolution of such process.

Morphological alterations and gene and protein expression profiling of bladder tumor cells after treatment with gemcitabine.

Chemical agents used in cancer therapy are associated with cell cycle arrest, activation or deactivation of mechanisms associated to DNA repair and apoptosis. However, due to the complexity of biological systems, the molecular mechanisms responsible for these activities are not fully understood. Thus, studies about gene and protein expression have shown promising results for understanding the mechanisms related to cellular responses and regression of cancer after chemotherapy. This study aimed to evaluate the gene and protein expression profiling in bladder transitional cell carcinoma (TCC) with different TP53 status after gemcitabine (1.56 μM) treatment. The RT4 (grade 1, TP53 wild type), 5637 (grade 2, TP53 mutated) and T24 (grade 3, TP53 mutated) cell lines were used. PCR arrays and mass spectrometry were used to analyze gene and protein expression, respectively. Morphological alterations were observed using scanning electron microscopy (SEM) and transmission electron microscopy (TEM). The results of PCR array showed that gemcitabine activity was mainly related to CDKN1A, GADD45A and SERTDA1 overexpression, and BAX overexpression only in the wild type TP53 cells. Mass spectrometry demonstrated that gemcitabine modulated the protein expression, especially those from genes related to apoptosis, transport of vesicles and stress response. Analyses using SEM and TEM showed changes in cell morphology independently on the cell line studied. The observed decreased number of microvillus suggests low contact among the cells and between cell and extracellular matrix; irregular forms might indicate actin cytoskeleton deregulation; and the reduction in the amount of organelles and core size might indicate reduced cellular metabolism. In conclusion, independently on TP53 status or grade of bladder tumor, gemcitabine modulated genes related to the cell cycle and apoptosis, that reflected in morphological changes indicative of future cell death.

Optical storage in azobenzene-containing epoxy polymers processed as Langmuir Blodgett films

In this study, azocopolymers containing different main-chain segments have been synthesized with diglycidyl ether of bisphenol A (DGEBA, DER 332, n=0.03) and the azochromophore Disperse Orange 3 (DO3) cured with twomonoamines, viz. benzylamine (BA) and m-toluidine (MT). The photoinduced birefringence was investigated in films produced with these azopolymers using the spin coating (SC) and Langmuir Blodgett (LB) techniques. In the LB films, birefringence increased with the content of azochromophore and the film thickness, as expected. The nanostructured nature of the LB films led to an enhanced birefringence and faster dynamics in the writing process, compared to the SC films. In summary, the combination of azocopolymers and the LBmethod may allow materials with tuned properties for various optical applications, including in biological systems were photoisomerization may be used to trigger actions such as drug delivery.