Soft Computing Approach for Optimization of siRNA Efficiency Prediction for Post-transcriptional Gene Silencing

Post-transcriptional gene silencing by RNA interference is mediated by small interfering RNA called siRNA. This gene silencing mechanism can be exploited therapeutically to a wide variety of disease-associated targets, especially in AIDS, neurodegenerative diseases, cholesterol and cancer on mice with the hope of extending these approaches to treat humans. Over the recent past, a significant amount of work has been undertaken to understand the gene silencing mediated by exogenous siRNA. The design of efficient exogenous siRNA sequences is challenging because of many issues related to siRNA. While designing efficient siRNA, target mRNAs must be selected such that their corresponding siRNAs are likely to be efficient against that target and unlikely to accidentally silence other transcripts due to sequence similarity. So before doing gene silencing by siRNAs, it is essential to analyze their off-target effects in addition to their inhibition efficiency against a particular target. Hence designing exogenous siRNA with good knock-down efficiency and target specificity is an area of concern to be addressed. Some methods have been developed already by considering both inhibition efficiency and off-target possibility of siRNA against agene. Out of these methods, only a few have achieved good inhibition efficiency, specificity and sensitivity. The main focus of this thesis is to develop computational methods to optimize the efficiency of siRNA in terms of “inhibition capacity and off-target possibility” against target mRNAs with improved efficacy, which may be useful in the area of gene silencing and drug design for tumor development. This study aims to investigate the currently available siRNA prediction approaches and to devise a better computational approach to tackle the problem of siRNA efficacy by inhibition capacity and off-target possibility. The strength and limitations of the available approaches are investigated and taken into consideration for making improved solution. Thus the approaches proposed in this study extend some of the good scoring previous state of the art techniques by incorporating machine learning and statistical approaches and thermodynamic features like whole stacking energy to improve the prediction accuracy, inhibition efficiency, sensitivity and specificity. Here, we propose one Support Vector Machine (SVM) model, and two Artificial Neural Network (ANN) models for siRNA efficiency prediction. In SVM model, the classification property is used to classify whether the siRNA is efficient or inefficient in silencing a target gene. The first ANNmodel, named siRNA Designer, is used for optimizing the inhibition efficiency of siRNA against target genes. The second ANN model, named Optimized siRNA Designer, OpsiD, produces efficient siRNAs with high inhibition efficiency to degrade target genes with improved sensitivity-specificity, and identifies the off-target knockdown possibility of siRNA against non-target genes. The models are trained and tested against a large data set of siRNA sequences. The validations are conducted using Pearson Correlation Coefficient, Mathews Correlation Coefficient, Receiver Operating Characteristic analysis, Accuracy of prediction, Sensitivity and Specificity. It is found that the approach, OpsiD, is capable of predicting the inhibition capacity of siRNA against a target mRNA with improved results over the state of the art techniques. Also we are able to understand the influence of whole stacking energy on efficiency of siRNA. The model is further improved by including the ability to identify the “off-target possibility” of predicted siRNA on non-target genes. Thus the proposed model, OpsiD, can predict optimized siRNA by considering both “inhibition efficiency on target genes and off-target possibility on non-target genes”, with improved inhibition efficiency, specificity and sensitivity. Since we have taken efforts to optimize the siRNA efficacy in terms of “inhibition efficiency and offtarget possibility”, we hope that the risk of “off-target effect” while doing gene silencing in various bioinformatics fields can be overcome to a great extent. These findings may provide new insights into cancer diagnosis, prognosis and therapy by gene silencing. The approach may be found useful for designing exogenous siRNA for therapeutic applications and gene silencing techniques in different areas of bioinformatics.

On the function of the Dictyostelium Argonaute A protein (AgnA) in epigenetic gene regulation

With molecular biology methods and bioinformatics, the Argonaute proteins in Dictyostelium discoideum were characterized, and the function of the AgnA protein in RNAi and DNA methylation was investigated, as well as cellular features. Also interaction partners of the PAZ-Piwi domain of AgnA (PAZ-PiwiAgnA) were discovered. The Dictyostelium genome encodes five Argonaute proteins, termed AgnA/B/C/D/E. The expression level of Argonaute proteins was AgnB/D/E > AgnA > AgnC. All these proteins contain the characteristic conserved of PAZ and Piwi domains. Fluorescence microscopy revealed that the overexpressed C-terminal GFP-fusion of PAZ-PiwiAgnA (PPWa-GFP) localized to the cytoplasm. Overexpression of PPWa-GFP leaded to an increased gene silencing efficiency mediated by RNAi but not by antisense RNA. This indicated that PAZ-PiwiAgnA is involved in the RNAi pathway, but not in the antisense pathway. An analysis of protein-protein interactions by a yeast-two-hybrid screen on a cDNA library from vegetatively grown Dictyostelium revealed that several proteins, such as EF2, EF1-I, IfdA, SahA, SamS, RANBP1, UAE1, CapA, and GpdA could interact with PAZ-PiwiAgnA. There was no interaction between PAZ-PiwiAgnA and HP1, HelF and DnmA detected by direct yeast-two-hybrid analysis. The fluorescence microscopy images showed that the overexpressed GFP-SahA or IfdA fusion proteins localized to both cytoplasm and nuclei, while the overexpressed GFP-SamS localized to the cytoplasm. The expression of SamS in AgnA knock down mutants was strongly down regulated on cDNA and mRNA level in, while the expression of SahA was only slightly down regulated. AgnA knock down mutants displayed defects in growth and phagocytosis, which suggested that AgnA affects also cell biological features. The inhibition of DNA methylation on DIRS-1 and Skipper retroelements, as well as the endogenous mvpB and telA gene, observed for the same strains, revealed that AgnA is involved in the DNA methylation pathway. Northern blot analysis showed that Skipper and DIRS-1 were rarely expressed in Ax2, but the expression of Skipper was upregulated in AgnA knock down mutants, while the expression of DIRS-1 was not changed. A knock out of the agnA gene failed even though the homologous recombination of the disruption construct occurred at the correct site, which indicated that there was a duplication of the agnA gene in the genome. The same phenomenon was also observed in ifdA knock out experiments.

Aberrant mRNA Transcripts and the Nonsense-Mediated Decay Proteins UPF2 and UPF3 Are Enriched in the Arabidopsis Nucleolus

The eukaryotic nucleolus is multifunctional and involved in the metabolism and assembly of many different RNAs and ribonucleoprotein particles as well as in cellular functions, such as cell division and transcriptional silencing in plants. We previously showed that Arabidopsis thaliana exon junction complex proteins associate with the nucleolus, suggesting a role for the nucleolus in mRNA production. Here, we report that the plant nucleolus contains mRNAs, including fully spliced, aberrantly spliced, and single exon gene transcripts. Aberrant mRNAs are much more abundant in nucleolar fractions, while fully spliced products are more abundant in nucleoplasmic fractions. The majority of the aberrant transcripts contain premature termination codons and have characteristics of nonsense-mediated decay (NMD) substrates. A direct link between NMD and the nucleolus is shown by increased levels of the same aberrant transcripts in both the nucleolus and in Up-frameshift (upf) mutants impaired in NMD. In addition, the NMD factors UPF3 and UPF2 localize to the nucleolus, suggesting that the Arabidopsis nucleolus is therefore involved in identifying aberrant mRNAs and NMD.

The role of satellite image-processing for national-scale estimates of gene flow from genetically modified crops: rapeseed in the UK as a model.

1. There is concern over the possibility of unwanted environmental change following transgene movement from genetically modified (GM) rapeseed Brassica napus to its wild and weedy relatives. 2. The aim of this research was to develop a remote sensing-assisted methodology to help quantify gene flow from crops to their wild relatives over wide areas. Emphasis was placed on locating sites of sympatry, where the frequency of gene flow is likely to be highest, and on measuring the size of rapeseed fields to allow spatially explicit modelling of wind-mediated pollen-dispersal patterns. 3. Remote sensing was used as a tool to locate rapeseed fields, and a variety of image-processing techniques was adopted to facilitate the compilation of a spatially explicit profile of sympatry between the crop and Brassica rapa. 4. Classified satellite images containing rapeseed fields were first used to infer the spatial relationship between donor rapeseed fields and recipient riverside B. rapa populations. Such images also have utility for improving the efficiency of ground surveys by identifying probable sites of sympatry. The same data were then also used for the calculation of mean field size. 5. This paper forms a companion paper to Wilkinson et al. (2003), in which these elements were combined to produce a spatially explicit profile of hybrid formation over the UK. The current paper demonstrates the value of remote sensing and image processing for large-scale studies of gene flow, and describes a generic method that could be applied to a variety of crops in many countries. 6. Synthesis and applications. The decision to approve or prevent the release of a GM cultivar is made at a national rather than regional level. It is highly desirable that data relating to the decision-making process are collected at the same scale, rather than relying on extrapolation from smaller experiments designed at the plot, field or even regional scale. It would be extremely difficult and labour intensive to attempt to carry out such large-scale investigations without the use of remote-sensing technology. This study used rapeseed in the UK as a model to demonstrate the value of remote sensing in assembling empirical information at a national level.

Triacylglycerol-rich lipoprotein-gene interactions in endothelial cells

Lipoproteins such as LDL (low-density lipoprotein) and oxidized LDL have potentially adverse effects on endothelial cells due to their ability to activate pro-inflammatory pathways regulated via the transcription factor NF-kappaB (nuclear factor kappaB). Triacylglycerol-rich lipoproteins (the chylomicrons, very-low-density lipoprotein and their respective remnant particles) have also been implicated in the induction of a pro-inflammatory phenotype and up-regulation of adhesion molecule expression. Although early studies supported the proposal that LPL (lipoprotein lipase)-mediated hydrolysis of TRLs (triglyceride-rich lipoproteins) at the endothelium could activate the NFkappaB pathway, more recent studies provide evidence of pro-and anti-inflammatory responses when cells are exposed to fatty acids of TRL particles. A large number of genes are up- and down-regulated when cells are exposed to TRL, with the net effect reflecting receptor- and nonreceptor-mediated pathways that are activated or inhibited depending on fatty acid type, the lipid and apolipoprotein composition of the TRL and the presence or absence of LPL. Early concepts of TRL particles as essentially pro-inflammatory stimuli to the endothelium provide an overly simplistic view of their impact on the vascular compartment.

Juvenile-mature wood transition in pine: correlation between wood properties and candidate gene expression profiles

There is a strong desire to exploit transcriptomics data from model species for the genetic improvement of non-model crops. Here, we use gene expression profiles from the commercial model Pinus taeda to identify candidate genes implicated in juvenile-mature wood transition in the non-model relative, P. sylvestris. Re-analysis of 'public domain' SAGE data from xylem tissues of P. taeda revealed 283 mature-abundant and 396 juvenile-abundant tags (P < 0.01), of which 70 and 137, respectively matched to genes with known function. Based on sequence similarity, we then isolated 16 putative homologues of genes that in P. taeda exhibited widest divergence in expression between juvenile and mature samples. Candidate expression levels in P. sylvestris were almost invariably differential between juvenile and mature woody tissue samples among two cohorts of five trees collected from the same seed source and selected for genetic uniformity by genetic distance analysis. However, the direction of differential expression was not always consistent with that described in the original P. taeda SAGE data. Correlation was observed between gene expression and juvenile-mature wood anatomical characteristics by OPLS analysis. Four candidates (alpha-tubulin, porin MIP1, lipid transfer protein and aquaporin like protein) apparently had greatest influence on the wood traits measured. Speculative function of these genes in relation to juvenile-mature wood transition is briefly explored. Thus, we demonstrate the feasibility of exploiting SAGE data from a model species to identify consistently differentially expressed candidates in a related non-model species.

Gene flow, risk assessment and the environmental release of transgenic plants

The release of genetically modified plants is governed by regulations that aim to provide an assessment of potential impact on the environment. One of the most important components of this risk assessment is an evaluation of the probability of gene flow. In this review, we provide an overview of the current literature on gene flow from transgenic plants, providing a framework of issues for those considering the release of a transgenic plant into the environment. For some plants gene flow from transgenic crops is well documented, and this information is discussed in detail in this review. Mechanisms of gene flow vary from plant species to plant species and range from the possibility of asexual propagation, short- or long-distance pollen dispersal mediated by insects or wind and seed dispersal. Volunteer populations of transgenic plants may occur where seed is inadvertently spread during harvest or commercial distribution. If there are wild populations related to the transgenic crop then hybridization and eventually introgression in the wild may occur, as it has for herbicide resistant transgenic oilseed rape (Brassica napus). Tools to measure the amount of gene flow, experimental data measuring the distance of pollen dispersal, and experiments measuring hybridization and seed survivability are discussed in this review. The various methods that have been proposed to prevent gene flow from genetically modified plants are also described. The current "transgenic traits'! in the major crops confer resistance to herbicides and certain insects. Such traits could confer a selective advantage (an increase in fitness) in wild plant populations in some circumstances, were gene flow to occur. However, there is ample evidence that gene flow from crops to related wild species occurred before the development of transgenic crops and this should be taken into account in the risk assessment process.

A rotating disc voltammetry study of the 1,8-dihydroxyanthraquinone mediated reduction of colloidal indigo

Colloidal indigo is reduced to an aqueous solution of leuco-indigo in a mediated two-electron process converting the water-insoluble dye into the water-soluble leuco form. The colloidal dye does not interact directly with the electrode surface, and to employ an electrochemical process for this reduction, the redox mediator 1,8-dihydroxyanthraquinone (1,8-DHAQ) is used to transfer electrons from the electrode to the dye. The mediated reduction process is investigated at a (500-kHz ultrasound-assisted) rotating disc electrode, and the quantitative analysis of voltammetric data is attempted employing the Digisim numerical simulation software package. At the most effective temperature, 353 K, the diffusion coefficient for 1,8-DHAQ is (0.84 +/- 0.08)x10(-9) m(2) s(-1), and it is shown that an apparently kinetically controlled reaction between the reduced form of the mediator and the colloidal indigo occurs within the diffusion layer at the electrode surface. The apparent bimolecular rate constant k (app)=3 mol m(-3) s(-1) for the rate law d[leuco-indigo]/dt = k(app) x [mediator] x [indigo] is determined and attributed to a mediator diffusion controlled dissolution of the colloid particles. The average particle size and the number of molecules per particles are estimated from the apparent bimolecular rate constant and confirmed by scanning electron microscopy.

The role of satellite image-processing for national-scale estimates of gene flow from genetically modified crops: rapeseed in the UK as a model

1.There is concern over the possibility of unwanted environmental change following transgene movement from genetically modified (GM) rapeseed Brassica napus to its wild and weedy relatives. 2. The aim of this research was to develop a remote sensing-assisted methodology to help quantify gene flow from crops to their wild relatives over wide areas. Emphasis was placed on locating sites of sympatry, where the frequency of gene flow is likely to be highest, and on measuring the size of rapeseed fields to allow spatially explicit modelling of wind-mediated pollen-dispersal patterns. 3. Remote sensing was used as a tool to locate rapeseed fields, and a variety of image-processing techniques was adopted to facilitate the compilation of a spatially explicit profile of sympatry between the crop and Brassica rapa. 4. Classified satellite images containing rapeseed fields were first used to infer the spatial relationship between donor rapeseed fields and recipient riverside B. rapa populations. Such images also have utility for improving the efficiency of ground surveys by identifying probable sites of sympatry. The same data were then also used for the calculation of mean field size. 5. This paper forms a companion paper to Wilkinson et al. (2003), in which these elements were combined to produce a spatially explicit profile of hybrid formation over the UK. The current paper demonstrates the value of remote sensing and image processing for large-scale studies of gene flow, and describes a generic method that could be applied to a variety of crops in many countries. 6.Synthesis and applications. The decision to approve or prevent the release of a GM cultivar is made at a national rather than regional level. It is highly desirable that data relating to the decision-making process are collected at the same scale, rather than relying on extrapolation from smaller experiments designed at the plot, field or even regional scale. It would be extremely difficult and labour intensive to attempt to carry out such large-scale investigations without the use of remote-sensing technology. This study used rapeseed in the UK as a model to demonstrate the value of remote sensing in assembling empirical information at a national level.

A genetic linkage map of microsatellite, gene-specific and morphological markers in diploid Fragaria

Diploid Fragaria provide a potential model for genomic studies in the Rosaceae. To develop a genetic linkage map of diploid Fragaria, we scored 78 markers (68 microsatellites, one sequence-characterised amplified region, six gene-specific markers and three morphological traits) in an interspecific F2 population of 94 plants generated from a cross of F.vesca f. semperflorens × F. nubicola. Co-segregation analysis arranged 76 markers into seven discrete linkage groups covering 448 cM, with linkage group sizes ranging from 100.3 cM to 22.9 cM. Marker coverage was generally good; however some clustering of markers was observed on six of the seven linkage groups. Segregation distortion was observed at a high proportion of loci (54%), which could reflect the interspecific nature of the progeny and, in some cases, the self-incompatibility of F. nubicola. Such distortion may also account for some of the marker clustering observed in the map. One of the morphological markers, pale-green leaf (pg) has not previously been mapped in Fragaria and was located to the mid-point of linkage group VI. The transferable nature of the markers used in this study means that the map will be ideal for use as a framework for additional marker incorporation aimed at enhancing and resolving map coverage of the diploid Fragaria genome. The map also provides a sound basis for linkage map transfer to the cultivated octoploid strawberry.

Global iron-dependent gene regulation in Escherichia coli - A new mechanism for iron homeostasis

Organisms generally respond to iron deficiency by increasing their capacity to take up iron and by consuming intracellular iron stores. Escherichia coli, in which iron metabolism is particularly well understood, contains at least 7 iron-acquisition systems encoded by 35 iron-repressed genes. This Fe-dependent repression is mediated by a transcriptional repressor, Fur ( ferric uptake regulation), which also controls genes involved in other processes such as iron storage, the Tricarboxylic Acid Cycle, pathogenicity, and redox-stress resistance. Our macroarray-based global analysis of iron- and Fur-dependent gene expression in E. coli has revealed several novel Fur-repressed genes likely to specify at least three additional iron- transport pathways. Interestingly, a large group of energy metabolism genes was found to be iron and Fur induced. Many of these genes encode iron- rich respiratory complexes. This iron- and Fur-dependent regulation appears to represent a novel iron-homeostatic mechanism whereby the synthesis of many iron- containing proteins is repressed under iron- restricted conditions. This mechanism thus accounts for the low iron contents of fur mutants and explains how E. coli can modulate its iron requirements. Analysis of Fe-55-labeled E. coli proteins revealed a marked decrease in iron- protein composition for the fur mutant, and visible and EPR spectroscopy showed major reductions in cytochrome b and d levels, and in iron- sulfur cluster contents for the chelator-treated wild-type and/or fur mutant, correlating well with the array and quantitative RT-PCR data. In combination, the results provide compelling evidence for the regulation of intracellular iron consumption by the Fe2+-Fur complex.

Dietary vitamin C down-regulates inflammatory gene expression in apoE4 smokers

The deleterious impact of cigarette smoking on cardiovascular health may be in part attributable to a free radical mediated proinflammatory response in circulating monocytes. In the current investigation, the impact of vitamin C supplementation on monocyte gene expression was determined in apoE4 smokers versus non-smokers. A total of 10 smokers and 11 non-smokers consumed 60 mg/day of vitamin C for four weeks and a fasting blood sample was taken at baseline and post-intervention for the determination of plasma vitamin C and monocyte gene expression profiles using cDNA array and real time PCR. In apoE4 smokers, supplementation resulted in a 43% increase in plasma vitamin C concentrations. Furthermore, a number of genes were differentially expressed more than 2-fold in response to treatment, including a downregulation of the proinflammatory mediators tumor necrosis factor (TNF) beta, TNF receptor, neurotrophin-3 growth factor receptor, and monocyte chemoattractant protein I receptor. The study has identified a number of molecular mechanisms underlying the benefit of vitamin C supplementation in smokers. (c) 2005 Elsevier Inc. All rights reserved.

Transfer of members of the genus Falcivibrio to the genus Mobiluncus, and emended description of the genus Mobiluncus

It has long been thought that the genera Mobiluncus and Falcivibrio contain the same organisms. Using a polyphasic approach, it was found that Mobiluncus curtisii and Mobiluncus mulieris were the same as Falcivibrio vaginalis and Falcivibrio grandis, respectively. As the genus name Mobiluncus takes precedence, it is proposed that F. vaginalis and F. grandis be transferred to the genus Mobiluncus. In agreement with previous studies, results from phenotypic tests did not support the separation of M. curtisii strains into its two subspecies, M. curtisii subsp. curtisii and M. curtisii subsp. holmesii. Phenotypic complexity within M. curtisii dictates that the species should be treated as a complex until more in-depth analyses of the species have been performed. Phylogenetic analyses, based on 16S rRNA gene sequences, demonstrated that the genus Mobiluncus was associated with Varibaculum cambriense and the two subspecies of Actinomyces neuii, and that A. neuii is only distantly related to Actinomyces sensu stricto.

Dietary long-chain n-3 PUFAs increase LPL gene expression in adipose tissue of subjects with an atherogenic lipoprotein phenotype

We sought to test the hypothesis that dietary long-chain n-3 PUFA (LC n-3 PUFA) in fish oil stimulate the gene expression of lipoprotein lipase (LPL) in human adipose tissue (AT). In a randomized, double blind, placebo-controlled, cross-over study, 51 male subjects expressing an atherogenic lipoprotein phenotype (ALP) had their diets supplemented with fish oil for 6 weeks. As we previously reported for this group, supplementation with LC n-3 PUFA produced a decrease in fasting plasma triglyceride (TG) (−35%, P < 0.05), attenuation of the postprandial TG response (area and incremental area under the curve; AUC and IAUC, P < 0.05), and a decrease in small, dense LDL. The present study extended these observations by showing that these changes were accompanied by a marked increase in the concentration of LPL mRNA in adipose tissue (AT-LPL mRNA, +55%, P = 0.003) and post-heparin LPL activity (PH-LPL, +31%, P = 0.036). There was also evidence of an association between LPL gene expression and polymorphism in the apolipoprotein E gene. We conclude that the favorable influence of dietary n-3 PUFA on the ALP may be mediated, in part, through an increase in the plasma activity and gene expression of lipoprotein lipase in human adipose tissue.

Variation in the FFAR1 gene modifies BMI, body composition and plasma lipids but not plasma insulin or Beta-cell function in overweight subjects

Background FFAR1 receptor is a long chain fatty acid G-protein coupled receptor which is expressed widely, but found in high density in the pancreas and central nervous system. It has been suggested that FFAR1 may play a role in insulin sensitivity, lipotoxicity and is associated with type 2 diabetes. Here we investigate the effect of three common SNPs of FFAR1 (rs2301151; rs16970264; rs1573611) on pancreatic function, BMI, body composition and plasma lipids. Methodology/Principal Findings For this enquiry we used the baseline RISCK data, which provides a cohort of overweight subjects at increased cardiometabolic risk with detailed phenotyping. The key findings were SNPs of the FFAR1 gene region were associated with differences in body composition and lipids, and the effects of the 3 SNPs combined were cumulative on BMI, body composition and total cholesterol. The effects on BMI and body fat were predominantly mediated by rs1573611 (1.06 kg/m2 higher (P = 0.009) BMI and 1.53% higher (P = 0.002) body fat per C allele). Differences in plasma lipids were also associated with the BMI-increasing allele of rs2301151 including higher total cholesterol (0.2 mmol/L per G allele, P = 0.01) and with the variant A allele of rs16970264 associated with lower total (0.3 mmol/L, P = 0.02) and LDL (0.2 mmol/L, P<0.05) cholesterol, but also with lower HDL-cholesterol (0.09 mmol/L, P<0.05) although the difference was not apparent when controlling for multiple testing. There were no statistically significant effects of the three SNPs on insulin sensitivity or beta cell function. However accumulated risk allele showed a lower beta cell function on increasing plasma fatty acids with a carbon chain greater than six. Conclusions/Significance Differences in body composition and lipids associated with common SNPs in the FFAR1 gene were apparently not mediated by changes in insulin sensitivity or beta-cell function.