Complex networks: the key to systems biology

Though introduced recently, complex networks research has grown steadily because of its potential to represent, characterize and model a wide range of intricate natural systems and phenomena. Because of the intrinsic complexity and systemic organization of life, complex networks provide a specially promising framework for systems biology investigation. The current article is an up-to-date review of the major developments related to the application of complex networks in biology, with special attention focused on the more recent literature. The main concepts and models of complex networks are presented and illustrated in an accessible fashion. Three main types of networks are covered: transcriptional regulatory networks, protein-protein interaction networks and metabolic networks. The key role of complex networks for systems biology is extensively illustrated by several of the papers reviewed.

Xylella fastidiosa gene expression analysis by DNA microarrays

Xylella fastidiosa genome sequencing has generated valuable data by identifying genes acting either on metabolic pathways or in associated pathogenicity and virulence. Based on available information on these genes, new strategies for studying their expression patterns, such as microarray technology, were employed. A total of 2,600 primer pairs were synthesized and then used to generate fragments using the PCR technique. The arrays were hybridized against cDNAs labeled during reverse transcription reactions and which were obtained from bacteria grown under two different conditions (liquid XDM2 and liquid BCYE). All data were statistically analyzed to verify which genes were differentially expressed. In addition to exploring conditions for X. fastidiosa genome-wide transcriptome analysis, the present work observed the differential expression of several classes of genes (energy, protein, amino acid and nucleotide metabolism, transport, degradation of substances, toxins and hypothetical proteins, among others). The understanding of expressed genes in these two different media will be useful in comprehending the metabolic characteristics of X. fastidiosa, and in evaluating how important certain genes are for the functioning and survival of these bacteria in plants.

Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

GOOAL: A routine based on Genetic Algorithms to track harmonics in electric power systems

Voltage and current waveforms of a distribution or transmission power system are not pure sinusoids. There are distortions in these waveforms that can be represented as a combination of the fundamental frequency, harmonics and high frequency transients. This paper presents a novel approach to identifying harmonics in power system distorted waveforms. The proposed method is based on Genetic Algorithms, which is an optimization technique inspired by genetics and natural evolution. GOOAL, a specially designed intelligent algorithm for optimization problems, was successfully implemented and tested. Two kinds of representations concerning chromosomes are utilized: binary and real. The results show that the proposed method is more precise than the traditional Fourier Transform, especially considering the real representation of the chromosomes.

Macrogeographic Genetic Variation in Broad-Snouted Caiman (Caiman latirostris)

Broad-snouted caiman`s (Caiman latirostris) geographic distribution comprises one of the widest latitudinal ranges among all crocodilians. In this study we analyzed the relationship between geographic distance (along the species latitudinal range) and genetic differentiation using DNA microsatellite loci developed for C. latirostris and Alligator mississippiensis. The results suggest that there is a consistent relationship between geographic distance and genetic differentiation; however, other biogeographical factors seem to be relevant. The Atlantic Chain (Serra do Mar) seems to be an effective geographic barrier, as well as the relatively narrow (<= 1.5 km) sea channel between Cardoso Island and the continent. In addition, coastal populations seem to have been well connected in recent geological time (Pleistocene 16,000 years ago) all along the eastern Brazilian coast. Further studies should focus on the Sao Francisco River drainage, which is still poorly known for this species. J. Exp. Zool. 30.9A:628-636, 2008. (C) 2008 Wiley-Liss, Inc.

The-2518 bp promoter polymorphism at CCL2/MCP1 influences susceptibility to mucosal but not localized cutaneous leishmaniasis in Brazil

Mucosal leishmaniasis (ML) follows localized cutaneous leishmaniasis (CL) caused by Leishmania braziliensis. Proinflammatory responses mediate CL self-healing but are exaggerated in ML Proinflammatory monocyte chemoattractant protein 1 (MCP-1; encoded by CCL2) is associated with CL We explore its role in CL/ML through analysis of the regulatory CCL2 -2518 bp promoter polymorphism in CL/ML population samples and families from Brazil. Genotype frequencies were compared among ML/CL cases and control groups using logistic regression and the family-based association test (FBAT). MCP-1 was measured in plasma and macrophages. The GG recessive genotype at CCL2 -2518 bp was more common in patients with ML (N = 67) than in neighborhood control (NC; N = 60) subjects (OR 1.78; 95% Cl 1.01-3.14; P = 0.045), than in NC combined with leishmanin skin-test positive (N = 60) controls (OR 4.40; 95% CI 1.42-13.65; P = 0.010), and than in controls combined with CL (N = 60) patients (OR 2.78; 95% CI 1.13-6.85; P = 0.045). No associations were observed for CL compared to any groups. FBAT (91 ML and 223 CL cases in families) confirmed recessive association of ML with allele G (Z = 2.679; P = 0.007). Higher levels of MCP-1 occurred in plasma (P = 0.03) and macrophages (P < 0.0001) from GG compared to AA individuals. These results suggest that high MCP-1 increases risk of ML (C) 2010 Elsevier B.V. All rights reserved.

Molecular characterization of the Hepatitis B virus genotypes in Colombia: A Bayesian inference on the genotype F

Hepatitis B is a worldwide health problem affecting about 2 billion people and more than 350 million are chronic carriers of the virus. Nine HBV genotypes (A to I) have been described. The geographical distribution of HBV genotypes is not completely understood due to the limited number of samples from some parts of the world. One such example is Colombia, in which few studies have described the HBV genotypes. In this study, we characterized HBV genotypes in 143 HBsAg-positive volunteer blood donors from Colombia. A fragment of 1306 bp partially comprising HBsAg and the DNA polymerase coding regions (S/POL) was amplified and sequenced. Bayesian phylogenetic analyses were conducted using the Markov Chain Monte Carlo (MCMC) approach to obtain the maximum clade credibility (MCC) tree using BEAST v.1.5.3. Of all samples, 68 were positive and 52 were successfully sequenced. Genotype F was the most prevalent in this population (77%) - subgenotypes F3 (75%) and Fib (2%). Genotype G (7.7%) and subgenotype A2 (15.3%) were also found. Genotype G sequence analysis suggests distinct introductions of this genotype in the country. Furthermore, we estimated the time of the most recent common ancestor (TMRCA) for each HBV/F subgenotype and also for Colombian F3 sequences using two different datasets: (i) 77 sequences comprising 1306 bp of S/POL region and (ii) 283 sequences comprising 681 bp of S/POL region. We also used two other previously estimated evolutionary rates: (i) 2.60 x 10(-4) s/s/y and (ii) 1.5 x 10(-5) s/s/y. Here we report the HBV genotypes circulating in Colombia and estimated the TMRCA for the four different subgenotypes of genotype F. (C) 2010 Elsevier B.V. All rights reserved.

Dynamics of Hepatitis D (delta) virus genotype 3 in the Amazon region of South America

Hepatitis delta virus (HDV) is widely distributed and associated with fulminant hepatitis epidemics in areas with high prevalence of HBV. Several studies performed in the 1980s showed data on HDV infection in South America, but there are no studies on the viral dynamics of this virus. The aim of this study was to conduct an evolutionary analysis of hepatitis delta genotype 3 (HDV/3) prevalent in South America: estimate its nucleotide substitution rate, determine the time of most recent ancestor (TMRCA) and characterize the epidemic history and evolutionary dynamics. Furthermore, we characterized the presence of HBV/HDV infection in seven samples collected from patients who died due to fulminant hepatitis from Amazon region in Colombia and included them in the evolutionary analysis. This is the first study reporting HBV and HDV sequences from the Amazon region of Colombia. Of the seven Colombian patients, five were positive for HBV-DNA and HDV-RNA. Of them, two samples were successfully sequenced for HBV (subgenotypes F3 and Fib) and the five samples HDV positive were classified as HDV/3. By using all HDV/3 available reference sequences with sampling dates (n = 36), we estimated the HDV/3 substitution rate in 1.07 x 10(-3) substitutions per site per year (s/s/y), which resulted in a time to the most recent common ancestor (TMRCA) of 85 years. Also, it was determined that HDV/3 spread exponentially from early 1950s to the 1970s in South America. This work discusses for the first time the viral dynamics for the HDV/3 circulating in South America. We suggest that the measures implemented to control HBV transmission resulted in the control of HDV/3 spreading in South America, especially after the important raise in this infection associated with a huge mortality during the 1950s up to the 1970s. The differences found among HDV/3 and the other HDV genotypes concerning its diversity raises the hypothesis of a different origin and/or a different transmission route. (C) 2011 Elsevier B.V. All rights reserved.

Serotonergic neurotransmission in the dorsal raphe nucleus recruits in situ 5-HT2A/2C receptors to modulate the post-ictal antinociception

The post-ictal immobility syndrome is followed by a significant increase in the nociceptive thresholds in animals and humans. The aim of this study was to assess the involvement of the dorsal raphe nucleus (DRN) in the post-ictal antinociception. The second aim was to study the role of serotonergic intrinsic mechanisms of the DRN in this hypo-algesic phenomenon. Pentylenetetrazole (PTZ), an ionophore GABA-mediated Cl- influx antagonist, was peripherally used to induce tonic-clonic seizures in Wistar rats. The nociceptive threshold was measured by the tail-flick test. Neurochemical lesions of the DRN, performed with microinjection of ibotenic acid (1.0 mu g/0.2 mu L), caused a significant decrease of tonic-clonic seizure-induced antinociception, suggesting the involvement of this nucleus in this antinociceptive Process. Microinjections of methysergide (1.0 and 5.0 mu g/0.2 mu L), a non-selective serotonergic receptor antagonist, into DRN caused a significant decrease in the post-ictal antinociception in seizing animals, compared to controls, in all post-ictal periods Presently studied. These findings were corroborated by microinjections of ketanserin (at 1.0 and 5.0 mu g/0.2 mu L) into DRN. Ketanserin is an antagonist with large affinity for 5-HT2A/2C serotonergic receptors, which, in this Case, Caused a significant decrease in the tail-flick latencies in seizing animals, compared to controls after the first 20 min following tonic-clonic convulsive reactions. These results indicate that serotonergic neurotransmission of the DRN neuronal clusters is involved in the organization of the post-ictal hypo-algesia. The 5-HT2A/2C receptors of DRN neurons seem to be critically involved in the increase of nociceptive thresholds following tonic-clonic seizures. (c) 2008 Elsevier Inc, All rights reserved.

Low genetic diversities of rabies virus populations within different hosts in Brazil

The low rates of nonsynonymous evolution observed in natural rabies virus (RABV) isolates are suggested to have arisen in association with the structural and functional constraints operating on the virus protein and the infection strategies employed by RABV within infected hosts to avoid strong selection by the immune response. In order to investigate the relationship between the genetic characteristics of RABV populations within hosts and the virus evolution, the present study examined the genetic heterogeneities of RABV populations within naturally infected dogs and foxes in Brazil, as well as those of bat RABV populations that were passaged once in suckling mice. Sequence analyses of complete RABV glycoprotein (G) genes showed that RABV populations within infected hosts were genetically highly homogeneous whether they were infected naturally or experimentally (nucleotide diversities of 0-0.95 x 10(-3)). In addition, amino acid mutations were randomly distributed over the entire region of the G protein, and the nonsynonymous/synonymous rate ratios (d(N)/d(S)) for the G protein gene were less than 1. These findings suggest that the low genetic diversities of RABV populations within hosts reflect the stabilizing selection operating on the virus, the infection strategies of the virus, and eventually, the evolutionary patterns of the virus. (C) 2009 Elsevier B.V. All rights reserved.

European 1: A globally important clonal complex of Mycobacterium bovis

We have identified a globally important clonal complex of Mycobacterium bovis by deletion analysis of over one thousand strains from over 30 countries. We initially show that over 99% of the strains of M. bovis, the cause of bovine tuberculosis, isolated from cattle in the Republic of Ireland and the UK are closely related and are members of a single clonal complex marked by the deletion of chromosomal region RDEu1 and we named this clonal complex European 1 (Eu1). Eu1 strains were present at less than 14% of French, Portuguese and Spanish isolates of M. bovis but are rare in other mainland European countries and Iran. However, strains of the Eu1 clonal complex were found at high frequency in former trading partners of the UK (USA, South Africa, New Zealand, Australia and Canada). The Americas, with the exception of Brazil, are dominated by the Eu1 clonal complex which was at high frequency in Argentina, Chile, Ecuador and Mexico as well as North America. Eu1 was rare or absent in the African countries surveyed except South Africa. A small sample of strains from Taiwan were non-Eu1 but, surprisingly, isolates from Korea and Kazakhstan were members of the Eu1 clonal complex. The simplest explanation for much of the current distribution of the Eu1 clonal complex is that it was spread in infected cattle, such as Herefords, from the UK to former trading partners, although there is evidence of secondary dispersion since. This is the first identification of a globally dispersed clonal complex M. bovis and indicates that much of the current global distribution of this important veterinary pathogen has resulted from relatively recent International trade in cattle. Crown Copyright (C) 2011 Published by Elsevier B.V. All rights reserved.

Evolutionary history of Dengue virus type 4: Insights into genotype phylodynamics

Dengue virus type 4 (DENV-4) circulates in tropical and subtropical countries from Asia and the Americas. Despite the importance of dengue virus distribution, little is known about the worldwide viral spread. Following a Bayesian phylogenetic approach we inferred the evolutionary history of 310 isolates sampled from 37 countries during the time period 1956-2008 and the spreading dynamics for genotypes I and II. The region (tropical rainforest biome) comprised by Malaysia-Thailand was the most likely ancestral area from which the serotype has originated and spread. Interestingly, cross-correlation analysis on demographic time series with the Asian sequences showed a statistically significant negative correlation that could be suggestive of competition among genotypes within the same serotype. (C) 2011 Elsevier B.V. All rights reserved.

Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58 Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates. (C) 2009 Elsevier B.V. All rights reserved.

Evolutionary dynamics of the immunodominant repeats of the Plasmodium vivax malaria-vaccine candidate circumsporozoite protein (CSP)

The circumsporozoite protein (CSP) of Plasmodium vivax, a major target for malaria vaccine development, has immunodominant B-cell epitopes mapped to central nonapeptide repeat arrays. To determine whether rearrangements of repeat motifs during mitotic DNA replication of parasites create significant CSP diversity under conditions of low effective meiotic recombination rates, we examined csp alleles from sympatric P. vivax isolates systematically sampled from an area of low malaria endemicity in Brazil over a period of 14 months. Nine unique csp types, comprising six different nona peptide repeats, were observed in 45 isolates analyzed. Identical or nearly identical repeats predominated in most arrays, consistent with their recent expansion. We found strong linkage disequilibrium at sites across the chromosome 8 segment flanking the csp locus, consistent with rare meiotic recombination in this region. We conclude that CSP repeat diversity may not be severely constrained by rare meiotic recombination in areas of low malaria endemicity. New repeat variants may be readily created by nonhomologous recombination even when meiotic recombination is rare, with potential implications for CSP-based vaccine development. (C) 2010 Elsevier B.V. All rights reserved.