Search

In this paper we describe applications of our ASAR package to functional, taxonomic and pathways analysis of metagenomes and propose future plans and perspectives. To illustrate an analytical potential of ASAR, we discuss outcomes of several projects. The main focus is made on metabolic plasticity of electrochemically active microbial communities and a potential role of integrated symbiotic bacterial interactions; antipathogenic properties of BES, manifested in its capacity to remove some pathogens from waste streams; and medical applications of this technology. We present ASAR-based metagenome analysis of evolving bacterial community from distillery waste over period of 36 months in BES environment as an example. Application of ASAR to personalised analyses of gut microbiome (GM) and the data interpretation based on publically available association studies are also discussed in this publication.

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.

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.

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.

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 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.

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.

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.

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.

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.