Tapping the wealth of microbial data in high-throughput metabolic model reconstruction

Genome-scale metabolic models are valuable tools in the metabolic engineering process, based on the ability of these models to integrate diverse sources of data to produce global predictions of organism behavior. At the most basic level, these models require only a genome sequence to construct, and once built, they may be used to predict essential genes, culture conditions, pathway utilization, and the modifications required to enhance a desired organism behavior. In this chapter, we address two key challenges associated with the reconstruction of metabolic models: (a) leveraging existing knowledge of microbiology, biochemistry, and available omics data to produce the best possible model; and (b) applying available tools and data to automate the reconstruction process. We consider these challenges as we progress through the model reconstruction process, beginning with genome assembly, and culminating in the integration of constraints to capture the impact of transcriptional regulation. We divide the reconstruction process into ten distinct steps: (1) genome assembly from sequenced reads; (2) automated structural and functional annotation; (3) phylogenetic tree-based curation of genome annotations; (4) assembly and standardization of biochemistry database; (5) genome-scale metabolic reconstruction; (6) generation of core metabolic model; (7) generation of biomass composition reaction; (8) completion of draft metabolic model; (9) curation of metabolic model; and (10) integration of regulatory constraints. Each of these ten steps is documented in detail.

Docking, virtual high throughput screening and in silico fragment-based drug design.

ABSTRACT The drug discovery process has been profoundly changed recently by the adoption of computational methods helping the design of new drug candidates more rapidly and at lower costs. In silico drug design consists of a collection of tools helping to make rational decisions at the different steps of the drug discovery process, such as the identification of a biomolecular target of therapeutical interest, the selection or the design of new lead compounds and their modification to obtain better affinities, as well as pharmacokinetic and pharmacodynamic properties. Among the different tools available, a particular emphasis is placed in this review on molecular docking, virtual high throughput screening and fragment-based ligand design.

Expression of the gene encoding the high-Km glucose transporter 2 by the early postimplantation mouse embryo is essential for neural tube defects associated with diabetic embryopathy.

AIMS/HYPOTHESIS: Excess glucose transport to embryos during diabetic pregnancy causes congenital malformations. The early postimplantation embryo expresses the gene encoding the high-Km GLUT2 (also known as SLC2A2) glucose transporter. The hypothesis tested here is that high-Km glucose transport by GLUT2 causes malformations resulting from maternal hyperglycaemia during diabetic pregnancy. MATERIALS AND METHODS: Glut2 mRNA was assayed by RT-PCR. The Km of embryo glucose transport was determined by measuring 0.5-20 mmol/l 2-deoxy[3H]glucose transport. To test whether the GLUT2 transporter is required for neural tube defects resulting from maternal hyperglycaemia, Glut2+/- mice were crossed and transient hyperglycaemia was induced by glucose injection on day 7.5 of pregnancy. Embryos were recovered on day 10.5, and the incidence of neural tube defects in wild-type, Glut2+/- and Glut2-/- embryos was scored. RESULTS: Early postimplantation embryos expressed Glut2, and expression was unaffected by maternal diabetes. Moreover, glucose transport by these embryos showed Michaelis-Menten kinetics of 16.19 mmol/l, consistent with transport mediated by GLUT2. In pregnancies made hyperglycaemic on day 7.5, neural tube defects were significantly increased in wild-type embryos, but Glut2+/- embryos were partially protected from neural tube defects, and Glut2-/- embryos were completely protected from these defects. The frequency of occurrence of wild-type, Glut2+/- and Glut2-/- embryos suggests that the presence of Glut2 alleles confers a survival advantage in embryos before day 10.5. CONCLUSIONS/INTERPRETATIONS: High-Km glucose transport by the GLUT2 glucose transporter during organogenesis is responsible for the embryopathic effects of maternal diabetes.

Characterization of the murine high Km glucose transporter GLUT2 gene and its transcriptional regulation by glucose in a differentiated insulin-secreting cell line.

In pancreatic beta-cells, the high Km glucose transporter GLUT2 catalyzes the first step in glucose-induced insulin secretion by glucose uptake. Expression of the transporter has been reported to be modulated by glucose either at the protein or mRNA levels. In this study we used the differentiated insulinoma cell line INS-1 which expresses high levels of GLUT2 and show that the expression of GLUT2 is regulated by glucose at the transcriptional level. By run-on transcription assays we showed that glucose induced GLUT2 gene transcription 3-4-fold in INS-1 cells which was paralleled by a 1.7-2.3-fold increase in cytoplasmic GLUT2 mRNA levels. To determine whether glucose regulatory sequences were present in the promoter region of GLUT2, we cloned and characterized a 1.4-kilobase region of mouse genomic DNA located 5' of the translation initiation site. By RNase protection assays and primer extension, we determined that multiple transcription initiation sites were present at positions -55, -64, and -115 from the first coding ATG and which were identified in liver, intestine, kidney, and beta-cells mRNAs. Plasmids were constructed with the mouse promoter region linked to the reporter gene chloramphenicol acetyltransferase (CAT), and transiently and stably transfected in the INS-1 cells. Glucose induced a concentration-dependent increase in CAT activity which reached a maximum of 3.6-fold at 20 mM glucose. Similar CAT constructs made of the human GLUT2 promoter region and the CAT gene displayed the same glucose-dependent increase in transcriptional activity when transfected into INS-1 cells. Comparison of the mouse and human promoter regions revealed sequence identity restricted to a few stretches of sequences which suggests that the glucose responsive element(s) may be conserved in these common sequences.

Isolation of genetic mutation leading to abnormal phenotype in human through ultra-high-throughput sequencing (UHTS)

The introduction of Next Generation Sequencing (NGS) facilitated the task of localizing DNA variation and identifying the genetic cause of yet unsolved Mendelian disorders. Using Whole Exome Capture method and NGS, we identified the causative genetic aberration responsible for a number of monogenic disorders previously undetermined. Due to the novelty of the NGS method we benchmarked different algorithms to assess their merits and defects. This allowed us to establish a pipeline that we successfully used to pinpoint genes responsible for a form of West's syndrome, a Complex Intellectual Disability syndrome associated with patellar dislocation and celiac disease, and correcting some erroneous molecular diagnosis of Alport's syndrome in a Saudi Arabian family.

Analysis of HIV-1 expression level and sense of transcription by high-throughput sequencing of the infected cell.

Next-generation sequencing offers an unprecedented opportunity to jointly analyze cellular and viral transcriptional activity without prerequisite knowledge of the nature of the transcripts. SupT1 cells were infected with a vesicular stomatitis virus G envelope protein (VSV-G)-pseudotyped HIV vector. At 24 h postinfection, both cellular and viral transcriptomes were analyzed by serial analysis of gene expression followed by high-throughput sequencing (SAGE-Seq). Read mapping resulted in 33 to 44 million tags aligning with the human transcriptome and 0.23 to 0.25 million tags aligning with the genome of the HIV-1 vector. Thus, at peak infection, 1 transcript in 143 is of viral origin (0.7%), including a small component of antisense viral transcription. Of the detected cellular transcripts, 826 (2.3%) were differentially expressed between mock- and HIV-infected samples. The approach also assessed whether HIV-1 infection modulates the expression of repetitive elements or endogenous retroviruses. We observed very active transcription of these elements, with 1 transcript in 237 being of such origin, corresponding on average to 123,123 reads in mock-infected samples (0.40%) and 129,149 reads in HIV-1-infected samples (0.45%) mapping to the genomic Repbase repository. This analysis highlights key details in the generation and interpretation of high-throughput data in the setting of HIV-1 cellular infection.

High throughput 3D super-resolution microscopy reveals Caulobacter crescentus in vivo Z-ring organization.

We created a high-throughput modality of photoactivated localization microscopy (PALM) that enables automated 3D PALM imaging of hundreds of synchronized bacteria during all stages of the cell cycle. We used high-throughput PALM to investigate the nanoscale organization of the bacterial cell division protein FtsZ in live Caulobacter crescentus. We observed that FtsZ predominantly localizes as a patchy midcell band, and only rarely as a continuous ring, supporting a model of "Z-ring" organization whereby FtsZ protofilaments are randomly distributed within the band and interact only weakly. We found evidence for a previously unidentified period of rapid ring contraction in the final stages of the cell cycle. We also found that DNA damage resulted in production of high-density continuous Z-rings, which may obstruct cytokinesis. Our results provide a detailed quantitative picture of in vivo Z-ring organization.

Ultra high throughput sequencing in human DNA variation detection: a comparative study on the NDUFA3-PRPF31 region.

BACKGROUND: Ultra high throughput sequencing (UHTS) technologies find an important application in targeted resequencing of candidate genes or of genomic intervals from genetic association studies. Despite the extraordinary power of these new methods, they are still rarely used in routine analysis of human genomic variants, in part because of the absence of specific standard procedures. The aim of this work is to provide human molecular geneticists with a tool to evaluate the best UHTS methodology for efficiently detecting DNA changes, from common SNPs to rare mutations. METHODOLOGY/PRINCIPAL FINDINGS: We tested the three most widespread UHTS platforms (Roche/454 GS FLX Titanium, Illumina/Solexa Genome Analyzer II and Applied Biosystems/SOLiD System 3) on a well-studied region of the human genome containing many polymorphisms and a very rare heterozygous mutation located within an intronic repetitive DNA element. We identify the qualities and the limitations of each platform and describe some peculiarities of UHTS in resequencing projects. CONCLUSIONS/SIGNIFICANCE: When appropriate filtering and mapping procedures are applied UHTS technology can be safely and efficiently used as a tool for targeted human DNA variations detection. Unless particular and platform-dependent characteristics are needed for specific projects, the most relevant parameter to consider in mainstream human genome resequencing procedures is the cost per sequenced base-pair associated to each machine.

Exploring the environmental diversity of kinetoplastid flagellates in the high-throughput DNA sequencing era

The class Kinetoplastea encompasses both free-living and parasitic species from a wide range of hosts. Several representatives of this group are responsible for severe human diseases and for economic losses in agriculture and livestock. While this group encompasses over 30 genera, most of the available information has been derived from the vertebrate pathogenic genera Leishmaniaand Trypanosoma. Recent studies of the previously neglected groups of Kinetoplastea indicated that the actual diversity is much higher than previously thought. This article discusses the known segment of kinetoplastid diversity and how gene-directed Sanger sequencing and next-generation sequencing methods can help to deepen our knowledge of these interesting protists.

IROme, a New High-Throughput Molecular Tool for the Diagnosis of Inherited Retinal Dystrophies-A Price Comparison with Sanger Sequencing.

The molecular diagnosis of retinal dystrophies (RD) is difficult because of genetic and clinical heterogeneity. Previously, the molecular screening of genes was done one by one, sometimes in a scheme based on the frequency of sequence variants and the number of exons/length of the candidate genes. Payment for these procedures was complicated and the sequential billing of several genes created endless paperwork. We therefore evaluated the costs of generating and sequencing a hybridization-based DNA library enriched for the 64 most frequently mutated genes in RD, called IROme, and compared them to the costs of amplifying and sequencing these genes by the Sanger method. The production cost generated by the high-throughput (HT) sequencing of IROme was established at CHF 2,875.75 per case. Sanger sequencing of the same exons cost CHF 69,399.02. Turnaround time of the analysis was 3 days for IROme. For Sanger sequencing, it could only be estimated, as we never sequenced all 64 genes in one single patient. Sale cost for IROme calculated on the basis of the sale cost of one exon by Sanger sequencing is CHF 8,445.88, which corresponds to the sale price of 40 exons. In conclusion, IROme is cheaper and faster than Sanger sequencing and therefore represents a sound approach for the diagnosis of RD, both scientifically and economically. As a drop in the costs of HT sequencing is anticipated, target resequencing might become the new gold standard in the molecular diagnosis of RD.

HTPSELEX--a database of high-throughput SELEX libraries for transcription factor binding sites.

HTPSELEX is a public database providing access to primary and derived data from high-throughput SELEX experiments aimed at characterizing the binding specificity of transcription factors. The resource is primarily intended to serve computational biologists interested in building models of transcription factor binding sites from large sets of binding sequences. The guiding principle is to make available all information that is relevant for this purpose. For each experiment, we try to provide accurate information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, assembled clone sequences (concatemers) and complete sets of in vitro selected protein-binding tags. In addition, we offer in-house derived binding sites models. HTPSELEX also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols. The FTP site contains the trace archives and database flatfiles. The web server offers user-friendly interfaces for viewing individual entries and quality-controlled download of SELEX sequence libraries according to a user-defined sequencing quality threshold. HTPSELEX is available from ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex/ and http://www.isrec.isb-sib.ch/htpselex.

Development and Validation of Single Nucleotide Polymorphisms (SNPs) Markers from Two Transcriptome 454-Runs of Turbot (Scophthalmus maximus) Using High-Throughput Genotyping

The turbot (Scophthalmus maximus) is a commercially valuable flatfish and one of the most promising aquaculture species in Europe. Two transcriptome 454-pyrosequencing runs were used in order to detect Single Nucleotide Polymorphisms (SNPs) in genesrelated to immune response and gonad differentiation. A total of 866 true SNPs were detected in 140 different contigs representing 262,093 bp as a whole. Only one true SNP was analyzed in each contig. One hundred and thirteen SNPs out of the 140 analyzed were feasible (genotyped), while Ш were polymorphic in a wild population. Transition/transversion ratio (1.354) was similar to that observed in other fish studies. Unbiased gene diversity (He) estimates ranged from 0.060 to 0.510 (mean = 0.351), minimum allele frequency (MAF) from 0.030 to 0.500 (mean = 0.259) and all loci were in Hardy-Weinberg equilibrium after Bonferroni correction. A large number of SNPs (49) were located in the coding region, 33 representing synonymous and 16 non-synonymous changes. Most SNP-containing genes were related to immune response and gonad differentiation processes, and could be candidates for functional changes leading to phenotypic changes. These markers will be useful for population screening to look for adaptive variation in wild and domestic turbot

Development of a high throughput PCR to detect Coxiella burnetii and its application in a diagnostic laboratory over a 7-year period.

Q fever is a worldwide zoonotic infectious disease due to Coxiella burnetii. The clinical presentation may be acute (pneumonia and/or hepatitis) or chronic (most commonly endocarditis). Diagnosis mainly relies on serology and PCR. We therefore developed a quantitative real-time PCR. We first tested blindly its performance on various clinical samples and then, when thoroughly validated, we applied it during a 7-year period for the diagnosis of both acute and persistent C. burnetii infection. Analytical sensitivity (< 10 copies/PCR) was excellent. When tested blindly on 183 samples, the specificity of the PCR was 100% (142/142) and the sensitivity was 71% (29/41). The sensitivity was 88% (7/8) on valvular samples, 69% (20/29) on blood samples and 50% (2/4) on urine samples. This new quantitative PCR was then successfully applied for the diagnosis of acute Q fever and endovascular infection due to C. burnetii, allowing the diagnosis of Q fever in six patients over a 7-year period. During a local small cluster of cases, the PCR was also applied to blood from 1355 blood donors; all were negative confirming the high specificity of this test. In conclusion, we developed a highly specific method with excellent sensitivity, which may be used on sera for the diagnosis of acute Q fever and on various samples such as sera, valvular samples, aortic specimens, bone and liver, for the diagnosis of persistent C. burnetii infection.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Two high throughput technologies to detect segmental aneuploidies identify new Williams-Beuren syndrome patients with atypical deletions.

OBJECTIVE: To develop and compare two new technologies for diagnosing a contiguous gene syndrome, the Williams-Beuren syndrome (WBS). METHODS: The first proposed method, named paralogous sequence quantification (PSQ), is based on the use of paralogous sequences located on different chromosomes and quantification of specific mismatches present at these loci using pyrosequencing technology. The second exploits quantitative real time polymerase chain reaction (QPCR) to assess the relative quantity of an analysed locus. RESULTS: A correct and unambiguous diagnosis was obtained for 100% of the analysed samples with either technique (n = 165 and n = 155, respectively). These methods allowed the identification of two patients with atypical deletions in a cohort of 182 WBS patients. Both patients presented with mild facial anomalies, mild mental retardation with impaired visuospatial cognition, supravalvar aortic stenosis, and normal growth indices. These observations are consistent with the involvement of GTF2IRD1 or GTF2I in some of the WBS facial features. CONCLUSIONS: Both PSQ and QPCR are robust, easy to interpret, and simple to set up. They represent a competitive alternative for the diagnosis of segmental aneuploidies in clinical laboratories. They have advantages over fluorescence in situ hybridisation or microsatellites/SNP genotyping for detecting short segmental aneuploidies as the former is costly and labour intensive while the latter depends on the informativeness of the polymorphisms.