Positive Selection Results in Frequent Reversible Amino Acid Replacements in the G Protein Gene of Human Respiratory Syncytial Virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a ""flipflop'' phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background: Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor sigma(S) (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain. Results: Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms. Conclusions: The ppGpp and sigma(S) variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Identification of the Schistosoma mansoni TNF-Alpha Receptor Gene and the Effect of Human TNF-Alpha on the Parasite Gene Expression Profile

Background: Schistosoma mansoni is the major causative agent of schistosomiasis. The parasite takes advantage of host signals to complete its development in the human body. Tumor necrosis factor-alpha (TNF-alpha) is a human cytokine involved in skin inflammatory responses, and although its effect on the adult parasite's metabolism and egg-laying process has been previously described, a comprehensive assessment of the TNF-alpha pathway and its downstream molecular effects is lacking. Methodology/Principal Findings: In the present work we describe a possible TNF-alpha receptor (TNFR) homolog gene in S. mansoni (SmTNFR). SmTNFR encodes a complete receptor sequence composed of 599 amino acids, and contains four cysteine-rich domains as described for TNFR members. Real-time RT-PCR experiments revealed that SmTNFR highest expression level is in cercariae, 3.5 (+/- 0.7) times higher than in adult worms. Downstream members of the known human TNF-alpha pathway were identified by an in silico analysis, revealing a possible TNF-alpha signaling pathway in the parasite. In order to simulate parasite's exposure to human cytokine during penetration of the skin, schistosomula were exposed to human TNF-alpha just 3 h after cercariae-to-schistosomula in vitro transformation, and large-scale gene expression measurements were performed with microarrays. A total of 548 genes with significantly altered expression were detected, when compared to control parasites. In addition, treatment of adult worms with TNF-alpha caused a significantly altered expression of 1857 genes. Interestingly, the set of genes altered in adults is different from that of schistosomula, with 58 genes in common, representing 3% of altered genes in adults and 11% in 3 h-old early schistosomula. Conclusions/Significance: We describe the possible molecular elements and targets involved in human TNF-alpha effect on S. mansoni, highlighting the mechanism by which recently transformed schistosomula may sense and respond to this host mediator at the site of cercarial penetration into the skin.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Construct and Compare Gene Coexpression Networks with DAPfinder and DAPview

Background: DAPfinder and DAPview are novel BRB-ArrayTools plug-ins to construct gene coexpression networks and identify significant differences in pairwise gene-gene coexpression between two phenotypes. Results: Each significant difference in gene-gene association represents a Differentially Associated Pair (DAP). Our tools include several choices of filtering methods, gene-gene association metrics, statistical testing methods and multiple comparison adjustments. Network results are easily displayed in Cytoscape. Analyses of glioma experiments and microarray simulations demonstrate the utility of these tools. Conclusions: DAPfinder is a new friendly-user tool for reconstruction and comparison of biological networks.

Real-time PCR investigation on the expression of sboA and ituD genes in Bacillus spp

Aims: To investigate the expression of sboA and ituD genes among strains of Bacillus spp. at different pH and temperature. Methods and Results: Different Bacillus strains from the Amazon basin and Bacillus subtilis ATCC 19659 were investigated for the production of subtilosin A and iturin A by qRT-PCR, analysing sboA and ituD gene expression under different culture conditions. Amazonian strains presented a general gene expression level lower than B. subtilis ATCC 19659 for sboA. In contrast, when analysing the expression of ituD gene, the strains from the Amazon, particularly P40 and P45B, exhibited higher levels of expression. Changes in pH (6 and 8) and temperature (37 and 42 degrees C) caused a decrease in sboA expression, but increased ituD expression among strains from Amazonian environment. Conclusions: Temperature and pH have an important influence on the expression of genes sboA (subtilosin A) and ituD (iturin A) among Bacillus spp. The strains P40 and P45B can be useful for the production of antimicrobial peptide iturin A. Significance and Impact of the Study: Monitoring the expression of essential biosynthetic genes by qRT-PCR is a valuable tool for optimization of the production of antimicrobial peptides.

Characterisation of the effect of a simulated hydrocarbon spill on diazotrophs in mangrove sediment mesocosm

An analysis of the effect of an oil spill on mangrove sediments was carried out by contamination of mesocosms derived from two different mangroves, one with a history of contamination and one pristine. The association between N(2) fixers and hydrocarbon degradation was assessed using quantitative PCR (qPCR) for the genes rrs and nifH, nifH clone library sequencing and total petroleum hydrocarbon (TPH) quantification using gas chromatography. TPH showed that the microbial communities of both mangroves were able to degrade the hydrocarbons added; however, whereas the majority of oil added to the mesocosm derived from the polluted mangrove was degraded in the 75 days of the experiment, there was only partially degradation in the mesocosm derived from the pristine mangrove. qPCR showed that the addition of oil led to an increase in rrs gene copy numbers in both mesocosms, having almost no effect on the nifH copy numbers in the pristine mangrove. Sequencing of nifH clones indicated that the changes promoted by the oil in the polluted mangrove were greater than those observed in the pristine mesocosm. The main effect observed in the polluted mesocosm was the selection of a single phylotype which is probably adapted to the presence of petroleum. These results, together with previous reports, give hints about the relationship between N(2) fixation and hydrocarbon degradation in natural ecosystems.

EXPRESSION OF GENES RELATED TO MYOSTATIN SIGNALING DURING RAT SKELETAL MUSCLE LONGITUDINAL GROWTH

In this study we investigated the gene expression of proteins related to myostatin (MSTN) signaling during skeletal muscle longitudinal growth. To promote muscle growth, Wistar male rats were submitted to a stretching protocol for different durations (12, 24, 48, and 96 hours). Following this protocol, soleus weight and length and sarcomere number were determined. In addition, expression levels of the genes that encode MSTN, follistatin isoforms 288 and 315 (FLST288 and FLST315), follistatin-like 3 protein (FLST-L3), growth and differentiation factor-associated protein-1 (GASP-1), activin IIB receptor (ActIIB), and SMAD-7 were determined by real-time polymerase chain reaction. Prolonged stretching increased soleus weight, length, and sarcomere number. In addition, MSTN gene expression was increased at 12-24 hours, followed by a decrease at 96 hours when compared with baseline values. FLST isoforms, FLST-L3, and GASP-1 mRNA levels increased significantly over all time-points. ActIIB gene expression decreased quickly at 12-24 hours. SMAD-7 mRNA levels showed a late increase at 48 hours, which peaked at 96 hours. The gene expression pattern of inhibitory proteins related to MSTN signaling suggests a strong downregulation of this pathway in response to prolonged stretching. Muscle Nerve 40: 992-999, 2009

Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice

Lacchini S, Heimann AS, Evangelista FS, Cardoso L, Silva GJ, Krieger JE. Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice. Physiol Genomics 37: 225-230, 2009. First published March 3, 2009; doi:10.1152/physiolgenomics.90288.2008.-We tested the hypothesis that small changes in angiotensin I-converting enzyme (ACE) expression can alter the vascular response to injury. Male mice containing one, two, three, and four copies of the Ace gene with no detectable vascular abnormality or changes in blood pressure were submitted to cuff-induced femoral artery injury. Femoral thickening was higher in 3- and 4-copy mice (42.4 +/- 4.3% and 45.7 +/- 6.5%, respectively) compared with 1- and 2-copy mice (8.3 +/- 1.3% and 8.5 +/- 0.9%, respectively). Femoral ACE levels from control and injured vessels were assessed in 1- and 3-copy Ace mice, which represent the extremes of the observed response. ACE vascular activity was higher in 3- vs. 1-copy Ace mice (2.4-fold, P < 0.05) in the control uninjured vessel. Upon injury, ACE activity significantly increased in both groups [2.41-fold and 2.14-fold (P < 0.05) for 1- and 3-copy groups, respectively] but reached higher levels in 3- vs. 1-copy Ace mice (P < 0.05). Pharmacological interventions were then used as a counterproof and to indirectly assess the role of angiotensin II (ANG II) on this response. Interestingly, ACE inhibition (enalapril) and ANG II AT(1) receptor blocker (losartan) reduced intima thickening in 3-copy mice to 1-copy mouse values (P < 0.05) while ANG II treatment significantly increased intima thickening in 1-copy mice to 3-copy mouse levels (P < 0.05). Together, these data indicate that small physiologically relevant changes in ACE, not associated with basal vascular abnormalities or blood pressure levels, do influence the magnitude of cuff-induced neointima thickening in mice.

Endothelial nitric oxide synthase polymorphisms and adaptation of parasympathetic modulation to exercise training

SILVA, B. M., F. J. NEVES, M. V. NEGRÃO, C. R. ALVES, R. G. DIAS, G. B. ALVES, A. C. PEREIRA, M. Urbana A. RONDON, J. E. KRIEGER, C. E. NEGRÃO, and A. C. DA NOBREGA. Endothelial Nitric Oxide Synthase Polymorphisms and Adaptation of Parasympathetic Modulation to Exercise Training. Med. Sci. Sports Exerc., Vol. 43, No. 9, pp. 1611-1618, 2011. Purpose: There is a large interindividual variation in the parasympathetic adaptation induced by aerobic exercise training, which may be partially attributed to genetic polymorphisms. Therefore, we investigated the association among three polymorphisms in the endothelial nitric oxide gene (-786T>C, 4b4a, and 894G>T), analyzed individually and as haplotypes, and the parasympathetic adaptation induced by exercise training. Methods: Eighty healthy males, age 20-35 yr, were genotyped by polymerase chain reaction-restriction fragment length polymorphism analysis, and haplotypes were inferred using the software PHASE 2.1. Autonomic modulation (i.e., HR variability and spontaneous baroreflex sensitivity) and peak oxygen consumption ((V) over dotO(2peak)) were measured before and after training (running, moderate to severe intensity, three times per week, 60 min.day(-1), during 18 wk). Results: Training increased (V) over dotO(2peak) (P < 0.05) and decreased mean arterial pressure (P < 0.05) in the whole sample. Subjects with the -786C polymorphic allele had a significant reduction in baroreflex sensitivity after training (change: wild type (-786TT) = 2% +/- 89% vs polymorphic (-786TC/CC) = -28% +/- 60%, median +/- quartile range, P = 0.03), and parasympathetic modulation was marginally reduced in subjects with the 894T polymorphic allele (change: wild type (894GG) = 8% +/- 67% vs polymorphic (894GT/TT) = -18% +/- 59%, median +/- quartile range, P = 0.06). Furthermore, parasympathetic modulation percent change was different between the haplotypes containing wild-type alleles(-786T/4b/894G) and polymorphic alleles at positions -786 and 894 (-786C/4b/894T) (-6% +/- 56% vs -41% +/- 50%, median T quartile range, P = 0.04). Conclusions: The polymorphic allele at position -786 and the haplotype containing polymorphic alleles at positions -786 and 894 in the endothelial nitric oxide gene were associated with decreased parasympathetic modulation after exercise training.

Exercise training improves muscle vasodilatation in individuals with T786C polymorphism of endothelial nitric oxide synthase gene

Negrão M.V, Alves CR, Alves G.B, Pereira A.C, Dias R.G, Laterza M.C, Mota G.F, Oliveira E.M, Bassaneze V, Krieger J.E, Negrão C.E, Rondon M.U.P. Exercise training improves muscle vasodilatation in individuals with T786C polymorphism of endothelial nitric oxide synthase gene. Physiol Genomics 42A: 71-77, 2010. First published July 6, 2010; doi:10.1152/physiolgenomics.00145.2009.-Allele T at promoter region of the eNOS gene has been associated with an increase in coronary disease mortality, suggesting that this allele increases susceptibility for endothelial dysfunction. In contrast, exercise training improves endothelial function. Thus, we hypothesized that: 1) Muscle vasodilatation during exercise is attenuated in individuals homozygous for allele T, and 2) Exercise training improves muscle vasodilatation in response to exercise for TT genotype individuals. From 133 preselected healthy individuals genotyped for the T786C polymorphism, 72 participated in the study: TT (n = 37; age 27 +/- 1 yr) and CT + CC (n = 35; age 26 +/- 1 yr). Forearm blood flow (venous occlusion plethysmography) and blood pressure (oscillometric automatic cuff) were evaluated at rest and during 30% handgrip exercise. Exercise training consisted of three sessions per week for 18 wk, with intensity between anaerobic threshold and respiratory compensation point. Resting forearm vascular conductance (FVC, P = 0.17) and mean blood pressure (P = 0.70) were similar between groups. However, FVC responses during handgrip exercise were significantly lower in TT individuals compared with CT + CC individuals (0.39 +/- 0.12 vs. 1.08 +/- 0.27 units, P = 0.01). Exercise training significantly increased peak VO(2) in both groups, but resting FVC remained unchanged. This intervention significantly increased FVC response to handgrip exercise in TT individuals (P = 0.03), but not in CT + CC individuals (P = 0.49), leading to an equivalent FVC response between TT and CT + CC individuals (1.05 +/- 0.18 vs. 1.59 +/- 0.27 units, P = 0.27). In conclusion, exercise training improves muscle vasodilatation in response to exercise in TT genotype individuals, demonstrating that genetic variants influence the effects of interventions such as exercise training.

Molecular characterization of a miraculin-like gene differentially expressed during coffee development and coffee leaf miner infestation

The characterization of a coffee gene encoding a protein similar to miraculin-like proteins, which are members of the plant Kunitz serine trypsin inhibitor (STI) family of proteinase inhibitors (PIs), is described. PIs are important proteins in plant defence against insects and in the regulation of proteolysis during plant development. This gene has high identity with the Richadella dulcifica taste-modifying protein miraculin and with the tomato protein LeMir; and was named as CoMir (Coffea miraculin). Structural protein modelling indicated that CoMir had structural similarities with the Kunitz STI proteins, but suggested specific folding structures. CoMir was up-regulated after coffee leaf miner (Leucoptera coffella) oviposition in resistant plants of a progeny derived from crosses between C. racemosa (resistant) and C. arabica (susceptible). Interestingly, this gene was down-regulated during coffee leaf miner herbivory in susceptible plants. CoMir expression was up-regulated after abscisic acid application and wounding stress and was prominent during the early stages of flower and fruit development. In situ hybridization revealed that CoMir transcripts accumulated in the anther tissues that display programmed cell death (tapetum, endothecium and stomium) and in the metaxylem vessels of the petals, stigma and leaves. In addition, the recombinant protein CoMir shows inhibitory activity against trypsin. According to the present results CoMir may act in proteolytic regulation during coffee development and in the defence against L. coffeella. The similarity of CoMir with other Kunitz STI proteins and the role of CoMir in plant development and plant stress are discussed.

The Penicillium chrysogenum aclA gene encodes a broad-substrate-specificity acyl-coenzyme A ligase involved in activation of adipic acid, a side-chain precursor for cephem antibiotics

Activation of the cephalosporin side-chain precursor to the corresponding CoA-thioester is an essential step for its incorporation into the P-lactam backbone. To identify an acyl-CoA ligase involved in activation of adipate, we searched in the genome database of Penicillium chrysogenum for putative structural genes encoding acyl-CoA ligases. Chemostat-based transcriptome analysis was used to identify the one presenting the highest expression level when cells were grown in the presence of adipate. Deletion of the gene renamed aclA, led to a 32% decreased specific rate of adipate consumption and a threefold reduction of adipoyl-6-aminopenicillanic acid levels, but did not affect penicillin V production. After overexpression in Escherichia coli, the purified protein was shown to have a broad substrate range including adipate. Finally, protein-fusion with cyan-fluorescent protein showed co-localization with microbody-borne acyl-transferase. Identification and functional characterization of aclA may aid in developing future metabolic engineering strategies for improving the production of different cephalosporins. (C) 2009 Elsevier Inc. All rights reserved.