Search

Functional coupling and comparative genomics analysis have been applied to study functional associations of orthologs of enterococcal cAD1 sex pheromone (P13268) known to be responsible for biofilm formation, conjugative plasmid transfer and spreading of bacterial antibiotics resistance. cAD1 peptide pheromone is released from the membrane lipoprotein with the peptide precursor encoded by a gene cad (tr|C2JQE7). Our analysis of genomic neighbourhood of cad and motifs of the encoded polypeptide and its orthologs suggests a close functional association between cAD1 and ApbE protein (Q82Z24), a FMN insertion and trafficking facilitator. The cad and apbE orthologs were coupled in the genomes and ApbE-specific motifs for FMN covalent attachment were identified in cad-encoded protein sequence and its orthologs. These findings suggest a potential role of FMN-based reductase function of the cAD1 lipoprotein precursor in its processing and release of the active sex pheromone peptide. They may lead to a new approach in prevention of antibiotic resistance spread via targeting sex pheromone processing chaperones or by suppression of the FMN availability and covalent binding. This methods can be also applied to a controlled evolution of bacterial pathogenicity in microbial fuel cells, as the findings suggest the crosstalk between bacterial pathogenicity and bacterial electro-activity.

Background Breast cancer metastasis is a highly prevalent cause of death for European females. DNA microarray analysis has established that primary tumors, which remain localized, differ in gene expression from those that metastasize. Crossanalysis of these studies allow to revile the differences that may be used as predictive in the disease prognosis and therapy. Objective The aim of the project was to validate suggested prognostic and therapeutic markers using meta-analysis of data on gene expression in metastatic and primary breast cancer tumors. Method Data on relative gene expression values from 12 studies on primary breast cancer and breast cancer metastasis were retrieved from Genevestigator (Nebion) database. The results of the data meta-analysis were compared with results of literature mining for suggested metastatic breast cancer markers and vectors and consistency of their reported differential expression. Results Our analysis suggested that transcriptional expression of the COX2 gene is significantly downregulated in metastatic tissue compared to normal breast tissue, but is not downregulated in primary tumors compared with normal breast tissue and may be used as a differential marker in metastatic breast cancer diagnostics. RRM2 gene expression decreases in metastases when compared to primary breast cancer and could be suggested as a marker to trace breast cancer evolution. Our study also supports MMP1, VCAM1, FZD3, VEGFC, FOXM1 and MUC1 as breast cancer onset markers, as these genes demonstrate significant differential expression in breast neoplasms compared with normal breast tissue. Conclusion COX2 and RRM2 are suggested to be prominent markers for breast cancer metastasis. The crosstalk between upstream regulators of genes differentially expressed in primary breast tumors and metastasis also suggests pathways involving p53, ER1, ERB-B2, TNF and WNT, as the most promising regulators that may be considered for new complex drug therapeutic interventions in breast cancer metastatic progression.

Next-generation sequencing (NGS), particularly RNA-sequencing (RNA-Seq) technique, allows detection and quantification of different RNA transcripts in a tissue sample, and in our case toxin transcripts from snake venom glands. Using this approach, novel toxin transcripts can be detected and abundancies of different isoforms of each toxin measured. The analytical pipeline can be briefly outlined as follows. Isolation of mRNA from tissue under RNase-free condition is essential to keep mRNA intact before sequencing. After mRNA fragmentation, the adapters are added to both ends of the fragments to synthesize complementary cDNAs. The obtained cDNA library is then sequenced on Illumina HiSeq 2000 platform. Quality of millions of reads produced from the NGS is checked and the sequences corresponding to the adapters and low-quality reads are removed. Subsequently, the NGS data are subjected to the workflow of de novo assembly, quantification of expression levels, annotation of transcripts, and identification of ORFs, signal peptides, structurally conserved domains, and functional motifs. In this report we describe the listed methodological steps and techniques in details and refer to the platforms and software that may be adopted for similar studies.

Next-generation sequencing (NGS), particularly RNA-sequencing (RNA-Seq) technique, allows detection and quantification of different RNA transcripts in a tissue sample, and in our case toxin transcripts from snake venom glands. Using this approach, novel toxin transcripts can be detected and abundancies of different isoforms of each toxin measured. The analytical pipeline can be briefly outlined as follows. Isolation of mRNA from tissue under RNase-free condition is essential to keep mRNA intact before sequencing. After mRNA fragmentation, the adapters are added to both ends of the fragments to synthesize complementary cDNAs. The obtained cDNA library is then sequenced on Illumina HiSeq 2000 platform. Quality of millions of reads produced from the NGS is checked and the sequences corresponding to the adapters and low-quality reads are removed. Subsequently, the NGS data are subjected to the workflow of de novo assembly, quantification of expression levels, annotation of transcripts, and identification of ORFs, signal peptides, structurally conserved domains, and functional motifs. In this report we describe the listed methodological steps and techniques in details and refer to the platforms and software that may be adopted for similar studies.

SUMMARYBio-electrochemical systems can generate electricity by virtue of mature microbial consortia that gradually and spontaneously optimize performance. To evaluate selective enrichment of these electrogenic microbial communities, five, 3-electrode reactors were inoculated with microbes derived from rice wash wastewater and incubated under a range of applied potentials. Reactors were sampled over a 12-week period and DNA extracted from anodal, cathodal, and planktonic bacterial communities was interrogated using a custom-made bioinformatics pipeline that combined 16S and metagenomic samples to monitor temporal changes in community composition. Some genera that constituted a minor proportion of the initial inoculum dominated within weeks following inoculation and correlated with applied potential. For instance, the abundance ofGeobacterincreased from 423-fold to 766-fold between −350 mV and −50 mV, respectively. Full metagenomic profiles of bacterial communities were obtained from reactors operating for 12 weeks. Functional analyses of metagenomes revealed metabolic changes between different species of the dominant genus,Geobacter, suggesting that optimal nutrient utilization at the lowest electrode potential is achieved via genome rearrangements and a strong inter-strain selection, as well as adjustment of the characteristic syntrophic relationships. These results reveal a certain degree of metabolic plasticity of electrochemically active bacteria and their communities in adaptation to adverse anodic and cathodic environments.

The heterogeneity of human populations is a common consideration in describing and validating their various age-related features. Heterogeneity, in particular, amongst other factors, is used to explain the variability of mortality rates across the lifespan and deviations from an exponential growth at young and very old ages. A mathematical model that combines the population heterogeneity with the assumption that the mortality of each constituent subpopulation increases exponentially with age, has recently been shown to successfully reproduce the entire mortality pattern across the lifespan as well as its evolution over time. Furthermore, the analysis of time-evolution of the mortality pattern, performed by fitting the model to actual data of consecutive periods, confirms the applicability of the compensation law of mortality to each subpopulation and concludes on the evolution of the population towards homogenisation. In this work we aim to show that the heterogeneity of human populations is not only a convenient consideration for fitting mortality data but is indeed the actual structure of the population as reflected by the dynamics of its mortality over age and time. In particular, we demonstrate that the model of heterogeneous populations fits mortality data better than most of the other models if the data are taken for the entire lifespan and better than all other models if we consider only old ages. Also, we show that the model can reproduce seemingly contradicting observations in late-life mortality dynamics namely deceleration, levelling-off and mortality decline. Assuming that heterogeneity is reflected in genetic variations within the population, using Swedish mortality data for 20th century we show that the homogenisation of the population, observed in the model fits, can be associated with the evolution of allele frequencies.

Objective:\ud \ud Many different gene families are currently being investigated for their potential role in epilepsy and in the response to antiepileptic drugs. A common research challenge is identifying the members of a gene family that are most significantly dysregulated within the human epileptic focus, before taking them forward for resource‐intensive functional studies. Published data about transcriptomic changes within the human epileptic focus remains incomplete. A need exists for an accurate in silico system for the prediction of dysregulated genes within the epileptic focus. We present such a bioinformatic system. We demonstrate the validity of our approach by applying it to the solute carrier (SLC) gene family. There are >400 known SLCs. SLCs have never been systematically studied in epilepsy.\ud Methods:\ud \ud Using our in silico system, we predicted the SLCs likely to be dysregulated in the epileptic focus. We validated our in silico predictions by identifying ex vivo the SLCs dysregulated in epileptic foci, and determining the overlap between our in silico and ex vivo results. For the ex vivo analysis, we used a custom oligonucleotide microarray containing exon probes for all known SLCs to analyze 24 hippocampal samples obtained from surgery for pharmacoresistant mesial temporal lobe epilepsy and 24 hippocampal samples from normal postmortem controls.\ud Results:\ud \ud There was a highly significant (p < 9.99 × 10−7) overlap between the genes identified by our in silico and ex vivo strategies. The SLCs identified were either metal ion exchangers or neurotransmitter transporters, which are likely to play a part in epilepsy by influencing neuronal excitability.\ud Significance:\ud \ud The identified SLCs are most likely to mediate pharmacoresistance in epilepsy by enhancing the intrinsic severity of epilepsy, but further functional work will be needed to fully evaluate their role. Our successful in silico strategy can be adapted in order to prioritize genes relevant to epilepsy from other gene families.

Snake venom metalloproteinases (SVMPs) are the key enzymes in Russell's viper (RV) venom which target all important components of haemostasis, such as clotting factors, platelets, endothelial cells and basement membrane. The structural diversity of SVMPs contributes to the broad spectrum of biological activities. The aim of the study was to investigate the SVMP transcript profile to gain better insights into the characteristic clinical manifestations of the Myanmar Russell's viper (MRV) bites that distinguish it from the RVs of other habitats. Next generation sequencing (RNA-Seq) of mRNA from MRV venom glands (2 males and 1 female) was performed on an Illumina HiSeq2000 platform and then de novo assembled using Trinity software. A total of 59 SVMP contigs were annotated through a Blastn search against the serpent nucleotide database from NCBI. Among them, disintegrins were the most abundant transcripts (75%) followed by the P-III class SVMPs (25%). The P-II SVMPs were scarce (0.002%), while no P-I SVMPs were detectable in the transcriptome. For detailed structural analysis, contigs were conceptually translated and compared with amino acid sequences from other RVs and other vipers using Clustal Omega. The RTS-disintegrin (jerdostatin homolog) was the most abundant among transcripts corresponding to 5 disintegrin isoforms. From 10 isoforms of SVMPs, RVV-X, and Vipera lebetina apoptosis-inducing protease (VLAIP) homolog, hereby termed Daboia siamensis AIP (DSAIP), were found to be highly expressed. Venom protein analysis using SDS-PAGE followed by mass spectrometry revealed that the disintegrin was scarce, while the latter two SVMPs were abundant. These two proteins can contribute to severe clinical manifestations caused by MRV envenomation.

n-3 Fatty acids, flavonoids and resveratrol are well publicised for their beneficial effects on human health and wellbeing. Identifying common, underlying biological mechanisms targeted by these functional foods would therefore be informative for the public health sector for advising on nutritional health and disease, food and drug product development and consumer interest. The aim of this study was to explore the potential effects of gene expression changes associated withn-3 fatty acids EPA and DHA, flavonoids and resveratrol on modifying biological systems and disease pathways. To test this, publicly available human microarray data for significant gene expression changes associated with dietary intervention with EPA/DHA, flavonoids and resveratrol was subjected to pathway analysis and significance testing for overlap with signals from genome-wide association studies (GWAS) for common non-communicable diseases and biological functions. There was an enrichment of genes implicated in immune responses and disease pathways which was common to all of the treatment conditions tested. Analysis of biological functions and disease pathways indicated anti-tumorigenic properties for EPA/DHA. In line with this, significance testing of the intersection of genes associated with these functional foods and GWAS hits for common biological functions (ageing and cognition) and non-communicable diseases (breast cancer, CVD, diabesity, neurodegeneration and psychiatric disorders) identified significant overlap between the EPA/DHA and breast cancer gene sets. Dietary intervention with EPA/DHA, flavonoids and resveratrol can target important biological and disease pathways suggesting a potentially important role for these bioactive compounds in the prevention and treatment of dietary-related diseases.

Introduction Invasive pneumococcal disease (IPD), caused by Streptococcus pneumoniae, is a leading cause of pneumonia, meningitis and septicaemia worldwide, with increased morbidity and mortality in HIV-infected children. Objectives We aimed to compare peripheral blood expression profiles between HIV-infected and uninfected children with pneumococcal meningitis and controls, and between survivors and non-survivors, in order to provide insight into the host inflammatory response leading to poorer outcomes. Design and setting Prospective case–control observational study in a tertiary hospital in Malawi Participants Children aged 2 months to 16 years with pneumococcal meningitis or pneumonia. Methods We used the human genome HGU133A Affymetrix array to explore differences in gene expression between cases with pneumococcal meningitis (n=12) and controls, and between HIV-infected and uninfected cases, and validated gene expression profiles for 34 genes using real-time quantitative PCR (RT-qPCR) in an independent set of cases with IPD (n=229) and controls (n=13). Pathway analysis was used to explore genes differentially expressed. Results Irrespective of underlying HIV infection, cases showed significant upregulation compared with controls of the following: S100 calcium-binding protein A12 (S100A12); vanin-1 (VNN1); arginase, liver (ARG1); matrix metallopeptidase 9 (MMP9); annexin A3 (ANXA3); interleukin 1 receptor, type II (IL1R2); CD177 molecule (CD177); endocytic adaptor protein (NUMB) and S100 calcium-binding protein A9 (S100A9), cytoskeleton-associated protein 4 (CKAP4); and glycogenin 1 (GYG1). RT-qPCR confirmed differential expression in keeping with microarray results. There was no differential gene expression in HIV-infected compared with HIV-uninfected cases, but there was significant upregulation of folate receptor 3 (FOLR3), S100A12 in survivors compared with non-survivors. Conclusion Children with IPD demonstrated increased expression in genes regulating immune activation, oxidative stress, leucocyte adhesion and migration, arginine metabolism, and glucocorticoid receptor signalling.

The SVA family of hominid-specific non-LTR retrotransposon comprises the youngest group of transposable elements in the human genome. The propagation of the most ancient SVA subfamily took place about 13.5 million years ago, and the youngest SVA subfamily appeared in the human genome after the human/chimpanzee divergence. Functional analysis of genes associated with SVA insertions demonstrated their link to multiple ontological categories, with one of the major categories being attributed to brain function. Further analysis of this subset demonstrated that SVA elements expanded their presence in the human genome at different stages of hominoid evolution and were associated with progressively evolving behavioral features that indicate a potential impact of SVA propagation on the cognitive ability of a modern human. Our analysis suggests a potential role of SVAs in the evolution of human central nervous system and especially in the emergence of functional trends relevant to social and parental behavior. Coevolution of behavioral features and reproductive functions are suggested by our analysis and discussed.

Idiopathic nephrotic syndrome (INS) is caused by renal diseases that increase the permeability of the glomerular filtration barrier without evidence of a specific systemic cause. The aim of the present work was to reveal inherent molecular features of INS in children using combined urinary proteomics and metabolomics profiling. In this study, label-free mass spectrometric analysis of urinary proteins and small molecule metabolites was carried out in 12 patients with INS versus 12 sex- and age-matched control subjects with normal renal function. Integration and biological interpretation of obtained results were carried out by Ingenuity IPA software. Validation of obtained proteomics data was carried out by Western blot method. Proteomics data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD000765. This study indicates for the first time that paediatric INS is associated with up-regulation of afamin, hydroxyphenylacetate and uridine, and concomitant down-regulation in glutamine and phenylalanine levels, and many of these molecular species were previously shown to be involved in oxidative stress. Further studies in larger patient population are underway to investigate the role of oxidative stress in renal injury in paediatric INS.

NAMPT, also known as PBEF and visfatin, can act extracellularly as a cytokine-like molecule or intracellularly as a NAMPT, regulating NAD biosynthesis in the NAD salvage pathway. Inhibitors of NAMPT have anti-inflammatory and anticancer activity and are finding use as therapeutic agents. In view of the importance of NAD metabolism in neutrophil function, we determined the effects of NAMPT inhibition on a variety of neutrophil functions associated with their role in host protection against infections. Incubation of human neutrophils with the NAMPT inhibitor APO866 decreased neutrophil NAD(P)/H levels in a dose- and time-dependent manner but without a concomitant change in cell viability. NAMPT inhibition did not affect the expression of a number of cell-surface receptors involved in adhesion and opsono-phagocytosis, but the respiratory burst was decreased significantly. Whereas opsono-phagocytosis of Staphylococcus aureus was unaffected by NAMPT inhibition, intraphagosomal oxidant production was decreased. However, the killing efficiency of neutrophils was unaffected. These data indicate that therapeutic NAMPT inhibition is unlikely to have deleterious effects on host protection against infections, in spite of this ability to down-regulate neutrophil respiratory burst activity significantly.

Anthropogenic endocrine disruptors now contaminate all environments globally, with concomitant deleterious effects across diverse taxa. While most studies on endocrine disruption (ED) have focused on vertebrates, the superimposition of male sexual characteristics in the female dogwhelk, Nucella lapillus (imposex), caused by organotins, provides one of the most clearcut ecological examples of anthropogenically induced ED in aquatic ecosystems. To identify the underpinning mechanisms of imposex for this 'nonmodel' species, we combined Roche 454 pyrosequencing with custom oligoarray fabrication inexpensively to both generate gene models and identify those responding to chronic tributyltin (TBT) treatment. The results supported the involvement of steroid, neuroendocrine peptide hormone dysfunction and retinoid mechanisms, but suggested additionally the involvement of putative peroxisome proliferator-activated receptor (PPAR) pathways. Application of rosiglitazone, a well-known vertebrate PPARγ ligand, to dogwhelks induced imposex in the absence of TBT. Thus, while TBT-induced imposex is linked to the induction of many genes and has a complex phenotype, it is likely also to be driven by PPAR-responsive pathways, hitherto not described in invertebrates. Our findings provide further evidence for a common signalling pathway between invertebrate and vertebrate species that has previously been overlooked in the study of endocrine disruption. The authors would like to thank Prof. Neil Hall, Dr Margaret Hughes and Dr Delphine Lallias for technical advice. We also acknowledge the long-term support of NERC at both Bangor and Liverpool. This work was supported by Fundação para a Ciência e a Tecnologia (FCT) Ph.D. grant SFRH/BD/27711/ 2006 awarded to Sonia Pascoal and by Natural Environment Research Council Molecular Genetics Facility (NERC MGF) grant 294 awarded to Simon Creer. published

Horizontal gene transfer (HGT) can radically alter the genomes of microorganisms, providing the capacity to adapt to new lifestyles, environments, and hosts. However, the extent of HGT between eukaryotes is unclear. Using whole-genome, gene-by-gene phylogenetic analysis we demonstrate an extensive pattern of cross-kingdom HGT between fungi and oomycetes. Comparative genomics, including the de novo genome sequence of Hyphochytrium catenoides , a free-living sister of the oomycetes, shows that these transfers largely converge within the radiation of oomycetes that colonize plant tissues. The repertoire of HGTs includes a large number of putatively secreted proteins; for example, 7.6% of the secreted proteome of the sudden oak death parasite Phytophthora ramorum has been acquired from fungi by HGT. Transfers include gene products with the capacity to break down plant cell walls and acquire sugars, nucleic acids, nitrogen, and phosphate sources from the environment. Predicted HGTs also include proteins implicated in resisting plant defense mechanisms and effector proteins for attacking plant cells. These data are consistent with the hypothesis that some oomycetes became successful plant parasites by multiple acquisitions of genes from fungi.

Some patients with pharmacoresistant epilepsy undergo therapeutic resection of the epileptic focus. At least 12 large-scale microarray studies on brain tissue from epilepsy surgery have been published over the last 10 years, but they have failed to make a significant impact upon our understanding of pharmacoresistance, because (1) doubts have been raised about their reproducibility, (2) only a small number of the gene expression changes found in each microarray study have been independently validated and (3) the results of different studies have not been integrated to give a coherent picture of the genetic changes involved in epilepsy pharmacoresistance. To overcome these limitations, we (1) assessed the reproducibility of the microarray studies by calculating the overlap between lists of differentially regulated genes from pairs of microarray studies and determining if this was greater than would be expected by chance alone, (2) used an inter-study cross-validation technique to simultaneously verify the expression changes of large numbers of genes and (3) used the combined results of the different microarray studies to perform an integrative analysis based on enriched gene ontology terms, networks and pathways. Using this approach, we respectively (1) demonstrate that there are statistically significant overlaps between the gene expression changes in different publications, (2) verify the differential expression of 233 genes and (3) identify the biological processes, networks and genes likely to be most important in the development of pharmacoresistant epilepsy. Our analysis provides novel biologically plausible candidate genes and pathways which warrant further investigation to assess their causal relevance.

Microbes can defend their host against virulent infections, but direct evidence for the adaptive origin of microbe-mediated protection is lacking. Using experimental evolution of a novel, tripartite interaction, we demonstrate that mildly pathogenic bacteria (Enterococcus faecalis) living in worms (Caenorhabditis elegans) rapidly evolved to defend their animal hosts against infection by a more virulent pathogen (Staphylococcus aureus), crossing the parasitism–mutualism continuum. Host protection evolved in all six, independently selected populations in response to within-host bacterial interactions and without direct selection for host health. Microbe-mediated protection was also effective against a broad spectrum of pathogenic S. aureus isolates. Genomic analysis implied that the mechanistic basis for E. faecalis-mediated protection was through increased production of antimicrobial superoxide, which was confirmed by biochemical assays. Our results indicate that microbes living within a host may make the evolutionary transition to mutualism in response to pathogen attack, and that microbiome evolution warrants consideration as a driver of infection outcome.

Background Carbamazepine (CBZ) therapy is associated with hypersensitivity reactions, which can be multisystemic and sometimes fatal. The pathogenesis of these reactions is incompletely understood but heterologous immunity generated to viral antigens has been proposed as a potential mechanism. Our aim was to investigate changes in gene expression profiles of peripheral blood mononuclear cells (PBMCs) from patients with CBZ hypersensitivity to characterise the pathways involved in immune activation. Methods PBMCs were isolated from five individuals with a history of hypersensitivity to CBZ. The PBMCs were incubated for 24 h with CBZ, carbamazepine 10,11-epoxide (CBZE), cell culture medium, or tetanus toxoid. Expression profiles for mRNA and microRNA were generated with microarrays (Affymetrix, Santa Clara, CA, USA). Differential gene expression was undertaken by limma analysis in the R statistical software. Ingenuity pathway analysis (IPA) was used to define the molecular mechanisms in development of CBZ hypersensitivity. Findings Incubation of PBMCs with CBZ, CBZE, and tetanus toxoid led to significant differential expression (log 2 -fold change >1 or COX8A, IFI35, IFIT3, PARP9, PARP12, PTGES, RSAD2, USP18, USP41, miR433, miR455, miR3194, miR4723 , and miR345 . These genes were highly interconnected in functional networks and defined five top functional categories—namely, viral infection, psoriasis, inflammation of body region or organ, and antiviral response. Interpretation Incubation of PBMCs from CBZ hypersensitive individuals with CBZ or CBZE led to treatment-specific changes in gene expression associated with activation of antiviral and inflammatory pathways. These CBZ-specific and CBZE-specific mRNA and microRNA changes could potentially act as novel diagnostic biomarkers and treatment targets for CBZ hypersensitivity reactions, but further validation and confirmation of specificity in larger numbers of patients are needed. Funding VY was supported by the North West England Medical Research Council Fellowship Scheme (clinical pharmacology and therapeutics), which is funded by the Medical Research Council (grant number G1000417/94909), ICON, GlaxoSmithKline, AstraZeneca, and the Medical Evaluation Unit.

Background: Pre-term birth (PTB) remains the leading cause of infant mortality and morbidity. Its etiology is multifactorial, with a strong genetic component. Genetic predisposition for the two subtypes, spontaneous PTB with intact membranes (sPTB) and preterm prelabor rapture of membranes (PPROM), and differences between them, have not yet been systematically summarised. Methods and findings: Our literature search identified 15 association studies conducted in 3,600 women on 2175 SNPs in 274 genes. We used Ingenuity software to impute gene pathways and networks related to sPTB and PPROM. Detailed insight in the defined functional ontologies clearly separated integrated datasets for sPTB and PPROM. Our analysis of upstream regulators of genes suggests that glucocorticoid receptor (NR3C1), peroxisome proliferator activated receptor c (PPARG) and interferon regulating factor 3 (IRF3) may be sPTB specific. PPROM-specific genes may be regulated by estrogen receptor2 (ESR2) and signal transducer and activator of transcription (STAT1). The inflammatory transcription factor NFkB is linked to both sPTB and PPROM, however, their inflammatory response is distinctly different. Conclusions: Based on our analyses, we propose an autoimmune/hormonal regulation axis for sPTB, whilst pathways implicated in the etiology of PPROM include hematologic/coagulation function disorder, collagen metabolism, matrix degradation and local inflammation. Our hypothesis generating study has identified new candidate genes in the pathogenesis of PPROM and sPTB, which should be validated in large cohorts.

The genome of Rhodobacter capsulatus has been completely sequenced. It consists of a single chromosome containing 3.5 Mb and a circular plasmid of 134 kb. This effort, started in 1992, began with a fine-structure restriction map of an overlapping set of cosmids that covered the genome. Cosmid sequencing led to a gapped genome that was filled by primer walking on the chromosome and by using lambda clones. Methods had to be developed to handle strong stops in the high GC (68%) inserts. Annotation was done with the ERGO system at Integrated Genomics, as was the reconstruction of the cell's metabolism. It was possible to recognize 3709 orfs of which functional assignments could be made with high confidence to 2392 (65%). Unusual features include the presence of numerous cryptic phage genomes embedded in the chromosome.

In multicellular organisms, tight regulation of gene expression ensures appropriate tissue and organismal growth throughout development. Reversible phosphorylation of the RNA Polymerase II (RNAPII) C-terminal domain (CTD) is critical for the regulation of gene expression states, but how phosphorylation is actively modified in a developmental context remains poorly understood. Protein phosphatase 1 (PP1) is one of several enzymes that has been reported to dephosphorylate the RNAPII CTD. However, PP1's contribution to transcriptional regulation during animal development and the mechanisms by which its activity is targeted to RNAPII have not been fully elucidated. Here we show that the Drosophila orthologue of the PP1 Nuclear Targeting Subunit (dPNUTS) is essential for organismal development and is cell autonomously required for growth of developing tissues. The function of dPNUTS in tissue development depends on its binding to PP1, which we show is targeted by dPNUTS to RNAPII at many active sites of transcription on chromosomes. Loss of dPNUTS function or specific disruption of its ability to bind PP1 results in hyperphosphorylation of the RNAPII CTD in whole animal extracts and on chromosomes. Consistent with dPNUTS being a global transcriptional regulator, we find that loss of dPNUTS function affects the expression of the majority of genes in developing 1st instar larvae, including those that promote proliferative growth. Together, these findings shed light on the in vivo role of the PNUTS-PP1 holoenzyme and its contribution to the control of gene expression during early Drosophila development., Author Summary During development, cells rely on appropriate patterns of gene expression to regulate metabolism in order to meet cellular demands and maintain rapid tissue growth. Conversely, dysregulation of gene expression is critical in various disease states, such as cancer, and during ageing. A key mechanism that is ubiquitously employed to control gene expression is reversible phosphorylation, a molecular switch that is used to regulate the activity of the transcriptional machinery. Here we identify an enzyme that binds to and regulates the phosphorylation state of RNA Polymerase II, a central component of the general transcription machinery. We also show that an essential role of this enzyme is to support normal patterns of gene expression that facilitate organismal growth. These findings are not only of relevance to the understanding of normal enzyme function but may also assist in the development of therapeutic strategies for the treatment of aberrant patterns of gene expression that occur during ageing and disease progression.

In recent decades, molecular and cellular biology has benefited from numerous fascinating developments in experimental technique, generating an overwhelming amount of data on various biological objects and processes. This, in turn, has led biologists to look for appropriate tools to facilitate systematic analysis of data. Thus, the need for mathematical techniques, which can be used to aid the classification and understanding of this ever-growing body of experimental data, is more profound now than ever before. Mathematical modelling is becoming increasingly integrated into biological studies in general and into developmental biology particularly. This review outlines some achievements of mathematics as applied to developmental biology and demonstrates the mathematical formulation of basic principles driving morphogenesis. We begin by describing a mathematical formalism used to analyse the formation and scaling of morphogen gradients. Then we address a problem of interplay between the dynamics of morphogen gradients and movement of cells, referring to mathematical models of gastrulation in the chick embryo. In the last section, we give an overview of various mathematical models used in the study of the developmental cycle ofDictyostelium discoideum, which is probably the best example of successful mathematical modelling in developmental biology.

The Library of Apicomplexan Metabolic Pathways (LAMP, http://www.llamp.net) is a web database that provides near complete mapping from genes to the central metabolic functions for some of the prominent intracellular parasites of the phylum Apicomplexa. This phylum includes the causative agents of malaria, toxoplasmosis and theileriosis—diseases with a huge economic and social impact. A number of apicomplexan genomes have been sequenced, but the accurate annotation of gene function remains challenging. We have adopted an approach called metabolic reconstruction, in which genes are systematically assigned to functions within pathways/networks for Toxoplasma gondii, Neospora caninum, Cryptosporidium and Theileria species, and Babesia bovis. Several functions missing from pathways have been identified, where the corresponding gene for an essential process appears to be absent from the current genome annotation. For each species, LAMP contains interactive diagrams of each pathway, hyperlinked to external resources and annotated with detailed information, including the sources of evidence used. We have also developed a section to highlight the overall metabolic capabilities of each species, such as the ability to synthesize or the dependence on the host for a particular metabolite. We expect this new database will become a valuable resource for fundamental and applied research on the Apicomplexa.

Dietary restriction (DR), limiting nutrient intake from diet without causing malnutrition, delays the aging process and extends lifespan in multiple organisms. The conserved life-extending effect of DR suggests the involvement of fundamental mechanisms, although these remain a subject of debate. To help decipher the life-extending mechanisms of DR, we first compiled a list of genes that if genetically altered disrupt or prevent the life-extending effects of DR. We called these DR–essential genes and identified more than 100 in model organisms such as yeast, worms, flies, and mice. In order for other researchers to benefit from this first curated list of genes essential for DR, we established an online database called GenDR (http://genomics.senescence.info/diet/). To dissect the interactions of DR–essential genes and discover the underlying lifespan-extending mechanisms, we then used a variety of network and systems biology approaches to analyze the gene network of DR. We show that DR–essential genes are more conserved at the molecular level and have more molecular interactions than expected by chance. Furthermore, we employed a guilt-by-association method to predict novel DR–essential genes. In budding yeast, we predicted nine genes related to vacuolar functions; we show experimentally that mutations deleting eight of those genes prevent the life-extending effects of DR. Three of these mutants (OPT2, FRE6, and RCR2) had extended lifespan under ad libitum, indicating that the lack of further longevity under DR is not caused by a general compromise of fitness. These results demonstrate how network analyses of DR using GenDR can be used to make phenotypically relevant predictions. Moreover, gene-regulatory circuits reveal that the DR–induced transcriptional signature in yeast involves nutrient-sensing, stress responses and meiotic transcription factors. Finally, comparing the influence of gene expression changes during DR on the interactomes of multiple organisms led us to suggest that DR commonly suppresses translation, while stimulating an ancient reproduction-related process., Author Summary Dietary restriction has been shown to extend lifespan in diverse, evolutionarily distant species, yet its underlying mechanisms remain unknown. We first constructed a database of genes essential for the life-extending effects of dietary restriction in various model organisms and then studied their interactions using a variety of network and systems biology approaches. This enabled us to predict novel genes related to dietary restriction, which we validated experimentally in yeast. By comparing large-scale data compilations (interactomes and transcriptomes) from multiple organisms, we were able to condense this -omics information to the most conserved essential elements, eliminating species-specific adaptive responses. These results lead us to the rather surprising conclusion that lifespan extension by a restricted diet commonly may exploit an ancient rejuvenation process derived from gametogenesis.

Clinically, our ability to predict disease outcome for patients with early stage lung cancer is currently poor. To address this issue, tumour specimens were collected at surgery from non-small cell lung cancer (NSCLC) patients as part of the European Early Lung Cancer (EUELC) consortium. The patients were followed-up for three years post-surgery and patients who suffered progressive disease (PD, tumour recurrence, metastasis or a second primary) or remained disease-free (DF) during follow-up were identified. RNA from both tumour and adjacent-normal lung tissue was extracted from patients and subjected to microarray expression profiling. These samples included 36 adenocarcinomas and 23 squamous cell carcinomas from both PD and DF patients. The microarray data was subject to a series of systematic bioinformatics analyses at gene, network and transcription factor levels. The focus of these analyses was 2-fold: firstly to determine whether there were specific biomarkers capable of differentiating between PD and DF patients, and secondly, to identify molecular networks which may contribute to the progressive tumour phenotype. The experimental design and analyses performed permitted the clear differentiation between PD and DF patients using a set of biomarkers implicated in neuroendocrine signalling and allowed the inference of a set of transcription factors whose activity may differ according to disease outcome. Potential links between the biomarkers, the transcription factors and the genes p21/CDKN1A and Myc, which have previously been implicated in NSCLC development, were revealed by a combination of pathway analysis and microarray meta-analysis. These findings suggest that neuroendocrine-related genes, potentially driven through p21/CDKN1A and Myc, are closely linked to whether or not a NSCLC patient will have poor clinical outcome.

Olga VasievaInstitute of Integrative Biology, University of Liverpool, Liverpool, United Kingdom; Fellowship for the Interpretation of Genomes, Burr Ridge, IL, USAAbstract: Shwachman-Bodian-Diamond syndrome (SBDS) is linked to a mutation in a single gene. The SBDS proinvolved in RNA metabolism and ribosome-associated functions, but SBDS mutation is primarily linked to a defect in polymorphonuclear leukocytes unable to orient correctly in a spatial gradient of chemoattractants. Results of data mining and comparative genomic approaches undertaken in this study suggest that SBDS protein is also linked to tRNA metabolism and translation initiation. Analysis of crosstalk between translation machinery and cytoskeletal dynamics provides new insights into the cellular chemotactic defects caused by SBDS protein malfunction. The proposed functional interactions provide a new approach to exploit potential targets in the treatment and monitoring of this disease.Keywords: Shwachman-Bodian-Diamond syndrome, wybutosine, tRNA, chemotaxis, translation, genomics, gene proximity

Glutathione transferase Pi (GST-pi, GSTP) is known to strongly affect human susceptibility to several cancers, asthma and neurodegenerative disorders. As with other glutathione transferases, it catalyses the addition of reduced glutathione to electrophilic species, and it is important in metabolite detoxification. It also was shown to bind proteins and compounds containing iron and nitric oxide. Some of these interactions have developed in the course of evolution into regulatory pathways that back up the GST's most ancient catalytic functions and provide precise and diverse responses to chemical and redox stresses. An aim of this review is to summarise recent knowledge on GSTP's complementary functions in crosstalking pathways of conventional glutathione transfer, nitric oxide and lipid metabolism and ASK1-dependent stress response. This review will describe how these complex interactions affect regulation of cell respiration, biosynthesis of lung surfactant, organism's immunity and circadian rhythms. Integration of the data leads to a new interpretation of the role of GSTP in normal human physiology, pathology and an organism's susceptibility to diseases.

The ERGOTM (http://ergo.integratedgenomics.com/ERGO/) genome analysis and discovery suite is an integration of biological data from genomics, biochemistry, high-throughput expression profiling, genetics and peer-reviewed journals to achieve a comprehensive analysis of genes and genomes. Far beyond any conventional systems that facilitate functional assignments, ERGO combines pattern-based analysis with comparative genomics by visualizing genes within the context of regulation, expression profiling, phylogenetic clusters, fusion events, networked cellular pathways and chromosomal neighborhoods of other functionally related genes. The result of this multifaceted approach is to provide an extensively curated database of the largest available integration of genomes, with a vast collection of reconstructed cellular pathways spanning all domains of life. Although access to ERGO is provided only under subscription, it is already widely used by the academic community. The current version of the system integrates 500 genomes from all domains of life in various levels of completion, 403 of which are available for subscription.

We present a complete DNA sequence and metabolic analysis of the dominant oral bacterium Fusobacterium nucleatum . Although not considered a major dental pathogen on its own, this anaerobe facilitates the aggregation and establishment of several other species including the dental pathogens Porphyromonas gingivalis and Bacteroides forsythus . The F. nucleatum strain ATCC 25586 genome was assembled from shotgun sequences and analyzed using the ERGO bioinformatics suite ( http://www.integratedgenomics.com ). The genome contains 2.17 Mb encoding 2,067 open reading frames, organized on a single circular chromosome with 27% GC content. Despite its taxonomic position among the gram-negative bacteria, several features of its core metabolism are similar to that of gram-positive Clostridium spp., Enterococcus spp., and Lactococcus spp. The genome analysis has revealed several key aspects of the pathways of organic acid, amino acid, carbohydrate, and lipid metabolism. Nine very-high-molecular-weight outer membrane proteins are predicted from the sequence, none of which has been reported in the literature. More than 137 transporters for the uptake of a variety of substrates such as peptides, sugars, metal ions, and cofactors have been identified. Biosynthetic pathways exist for only three amino acids: glutamate, aspartate, and asparagine. The remaining amino acids are imported as such or as di- or oligopeptides that are subsequently degraded in the cytoplasm. A principal source of energy appears to be the fermentation of glutamate to butyrate. Additionally, desulfuration of cysteine and methionine yields ammonia, H 2 S, methyl mercaptan, and butyrate, which are capable of arresting fibroblast growth, thus preventing wound healing and aiding penetration of the gingival epithelium. The metabolic capabilities of F. nucleatum revealed by its genome are therefore consistent with its specialized niche in the mouth.

Background The Apicomplexa is a large phylum of intracellular parasites that show great diversity and adaptability in the various ecological niches they occupy. They are the causative agents of human and animal infections including malaria, toxoplasmosis and theileriosis, which have a huge economic and social impact. A number of apicomplexan genomes have been sequenced and are publicly available. However, the prediction of gene models and annotation of gene functions remains challenging. Methods We have utilised an approach called ‘metabolic reconstruction’, in which genes are systematically assigned to functions within pathways/networks [1-4]. Functional annotation and metabolic reconstruction was carried out using a semi-automatic approach, integrating genomic information with biochemical evidence from the literature. The functions were automatically assigned using a sequence similarity-based approach and protein motif information. Experimental evidence was also

Background The ageing of the worldwide population means there is a growing need for research on the biology of ageing. DNA damage is likely a key contributor to the ageing process and elucidating the role of different DNA repair systems in ageing is of great interest. In this paper we propose a data mining approach, based on classification methods (decision trees and Naive Bayes), for analysing data about human DNA repair genes. The goal is to build classification models that allow us to discriminate between ageing-related and non-ageing-related DNA repair genes, in order to better understand their different properties. Results The main patterns discovered by the classification methods are as follows: (a) the number of protein-protein interactions was a predictor of DNA repair proteins being ageing-related; (b) the use of predictor attributes based on protein-protein interactions considerably increased predictive accuracy of attributes based on Gene Ontology (GO) annotations; (c) GO terms related to "response to stimulus" seem reasonably good predictors of ageing-relatedness for DNA repair genes; (d) interaction with the XRCC5 (Ku80) protein is a strong predictor of ageing-relatedness for DNA repair genes; and (e) DNA repair genes with a high expression in T lymphocytes are more likely to be ageing-related. Conclusions The above patterns are broadly integrated in an analysis discussing relations between Ku, the non-homologous end joining DNA repair pathway, ageing and lymphocyte development. These patterns and their analysis support non-homologous end joining double strand break repair as central to the ageing-relatedness of DNA repair genes. Our work also showcases the use of protein interaction partners to improve accuracy in data mining methods and our approach could be applied to other ageing-related pathways.

During 2015-2017 we have conducted experiments in Japan to test the capacity of microbial fuel cells (MFC) to treat different types of wastewaters (swine farm, domestic, yeast fermentation, winery, etc.) and concomitantly collecting DNA samples from MFC anodic and planktonic bacterial communities. W e analyzed these metagenomes in UK, using our new bioinformatics tool (ASAR) that allow integration of phylogenetic and functional data. Characteristic MFC communities and the associated functional signatures were shown to reflect effective waste water treatment. We also found that the fraction of opportunistic pathogenic bacteria DNA was reduced in metagenomes from MFC communities during swine waste treatment. The highest loss was recorded for Enterobacteriaceae family (such as Yersinia, Vibrio, and Shigella). The abundance of virulent genes responsible for adhesion, secretion systems, invasion and intracellular survival, as well as antibiotic resistance, associated with Firmicutes and Actinobacteria phyla of Gram-positive bacteria, also decreased in the MFC residential metagenomes. Key metabolic functions were redistributed among bacteria on the anode and archaea in plankton. We propose to use MFC, inoculated with electroactive bacterial communities, for waste disinfection, and potentially for development of novel antibacterial therapies. This approach promises to be effective and economically justified, especially in cases of epidemics of enterobacteria-associated diseases, and common residential hospital pathogens such as Enterococcus.

Rubisco is the essential enzyme mediating the fixation of atmospheric CO(2) during photosynthesis. In cyanobacteria, Rubisco enzymes are densely packed and encapsulated in a specialized organelle known as the carboxysome. Well-defined Rubisco assembly and carboxysome formation are pivotal for efficient CO(2) fixation. Numerous chaperone proteins, including RbcX, are essential for proper protein folding and Rubisco assembly. In this study, we investigated the in vivo function of RbcX in the cyanobacterium Synechococcus elongatus PCC 7942 (Syn7942) using molecular, biochemical, and live-cell fluorescence imaging approaches. Our results show that genetic deletion of the rbcX gene affects Rubisco abundance, as well as carboxysome formation and spatial distribution. Moreover, RbcX appears as one component of the carboxysome and shows a dynamic interaction with Rubisco enzymes. These in vivo observations provide insight into the role of RbcX from Syn7942 in mediating carboxysome assembly. Understanding the molecular mechanism underlying Rubisco assembly and carboxysome biogenesis will provide essential information required for engineering functional CO(2)-fixing complexes in heterogeneous organisms, especially plants, with the aim of boosting photosynthesis and agricultural productivity.

In this case study, we validate a personalized approach to gut microbiome (GM) analysis and interpretation based on published association studies. We apply our ASAR data annotation and clustering package to a series of 10 sequenced GM’s from individuals of different ethnical and geographical backgrounds, age and health groups. The differentially presented and detectable taxonomic and functional signatures in each GM metagenome are used to predict the hosts’ characteristics via correlations established in published studies, and the predictions are being validated by available individual-associated metadata. We also test sensitivity of the routine annotation and data clustering pipeline to an individual and family-linked signatures in GM structure and functionalities, when applied to a limited number of varying samples. The number of samples was sufficient to demonstrate 2 main types of GM composition, based on Bacteroides or Prevotella as main abundant genera, limitation of a variety of taxa as a result of antibiotics application, clustering of family members’ GM metagenomes both in taxonomic and in functional space, individual signatures related to chronic diseases and pharmacological interventions, and elements of ethnicityrelated characteristics in the metagenomes. The method and logical algorithms of the analysis applied here may be utilized in rather computational pipeline for a personalized microbiome analyses, and their potential useful outputs and limitations are being discussed.

Background Prediction of the role of one protein in a network of regulatory and metabolic interactions involves analysis of many complex data sets in their interference, as no experimental source along is completely and generally reliable. To discriminate the most promising pathways related to a protein of interest from background noise we integrate multiple bioinformatics data by means of a graphical designer/database of protein-protein interactions supported by data mining and logical modeling of resulted networks.

Russell’s viper bites are a major public health problem in tropical and subtropical regions. In Myanmar, a Russell’s viper (Daboia russellii siamensis) bite has a 60% morbidity rate and 8.2% fatality rate. Most victims encounter severe bleeding, renal failure and capillary leakage and the bite can possibly lead to death. Snake venom metalloproteinases (SVMPs) are the major components of the Viperidae venom and all mentioned lethal effects of the bites are attributed to these. The only available and partially effective agent for the treatment of the toxic effects is antivenom. Antivenom therapy is not always effective towards small toxins, however, and it can also provoke an anaphylactic response. The development of new therapeutic approaches is becoming increasingly important therefore. For the analysis of SVMP transcripts from Myanmar Russell’s viper, Next-Generation Sequencing (NGS) of mRNA from venom glands derived from 2 male snakes and 1 female snake was performed on an Illumina HiSeq2000 platform. De novo assembly of the reads was performed using Trinity software and the transcripts were annotated through Blastn against the collection of NCBI nucleotide sequences defined by the key-words (‘venom’ and ‘serpents’) search. Blastx hit results against the UniProtKB/Swiss-Prot (swissprot) database were also used for annotation of the transcripts. The abundance distribution (in term of FPKM value) of SVMPs toxin transcripts: disintegrin (75%), P-III SVMPs (25%) and P-II SVMPs (0.002 %), were the same for both male and female samples. No P-I SVMP transcripts were detected in the present analysis. A comparison of the contents of SVMP transcripts in adult male and female venom glands showed some gender-related differences. For example, a disintegrin transcript isoform (Dis 1b) was highly expressed only in the female venom gland, and some P-III SVMP isoforms (P-III 6, 7a, 7b) were only expressed at low level in the male venom glands. The P-II SVMP transcripts were expressed as different isoforms in male and female animals, which could reflect a sex-dimorphism of viper venom biological activities. This finding would support a requirement to use combined venoms of both sexes for preparation of antivenom. In addition to SVMP transcripts, mRNAs of novel tripeptide SVMP inhibitors (SVMPIs) were also discovered. These endogenous inhibitors have potential as a new treatment modality for neutralization of the effect of SVMP toxins. Two major snake SVMPs, RVV-X and daborhagin, were purified from Myanmar Russell’s viper venom using a new purification strategy. Moreover, the two novel endogenous tripeptides identified in transcript analysis, pERW and pEKW were identified and isolated from the crude venom. Both purified SVMPs showed caseinolytic activity. Additionally, RVV-X displayed specific proteolytic activity towards gelatin and Daborhagin showed potent fibrinogenolytic activity. These activities were inhibited by metal chelators. Notably, synthetic versions of the peptide inhibitors, pERW and pEKW, completely inhibited the gelatinolytic and fibrinogenolytic activities of the respective SVMPs when used at 5 mM concentration (estimated molar ratio of SVMP to tripeptide was 1:500). These complete inhibitory effects suggest that these tripeptides deserve further study as candidates for new therapeutic treatment against Russell’s viper envenomation.