Dataset Information

Nègre2023 - Genome-Scale Metabolic Network of Penicillium rubens Wisconsin 54-1255

ABSTRACT: Using in silico modelling, such Genome-Scale Metabolic Network (GSMN), represents a promising approach to predicting and understanding the potential for specialised metabolite production in a given organism. To address these questions, we reconstructed a new high-quality GSMN for the Penicillium rubens Wisconsin 54-1255 strain, a commonly used model organism. Our reconstruction, iPrub22, adheres to current convention standards and quality criteria, incorporating updated functional annotations, orthology searches with different GSMN templates, data from previous reconstructions, and manual curation steps targeting basal and specialised metabolites. With a MEMOTE score of 74% and a metabolic coverage of 45%, iPrub22 includes 5,464 metabolites interconnected by 5,919 reactions, of which 5,033 are supported by at least one genomic sequence. Of the metabolites present in iPrub22, 13% are categorised as belonging to specialised metabolism.

SUBMITTER: Delphine Nègre

PROVIDER: MODEL2306150001 | BioModels | 2023-06-20

REPOSITORIES: BioModels

ACCESS DATA

Dataset's files

Source:

			Action	DRS
	MODEL2306150001?filename=FROG_report_iPrub22.omex	Other
	MODEL2306150001?filename=Features_reconstruction_and_model.html	Other
	MODEL2306150001?filename=MEMOTE_report_iPrub22.html	Other
	MODEL2306150001?filename=iPrub22.sbml	Other
	MODEL2306150001?filename=iPrub22_miniFROG.xlsx	Xlsx

Items per page:

1 - 5 of 5

Similar Datasets

Project description:Microarray technology provides a powerful tool for defining gene expression profiles of airway epithelium that lend insight into the pathogenesis of human airway disorders. The focus of this study was to establish rigorous quality control parameters to ensure that microarray assessment of the airway epithelium is not confounded by experimental artifact. Samples (total n=223) of trachea, large and small airway epithelium were collected by fiberoptic bronchoscopy of 144 individuals (42 healthy non-smokers, 49 healthy smokers, 11 symptomatic smokers, 22 smokers with lone emphysema with normal spirometry, and 20 smokers with COPD) were processed and hybridized to Affymetrix HG-U133 2.0 Plus microarrays. The pre- and post-chip quality control (QC) criteria established, included: (1) RNA quality, assessed by RNA Integrity Number (RIN) ≥7.0 using Agilent 2100 Bioanalyzer software; (2) cRNA transcript integrity, assessed by signal intensity ratio of GAPDH 3' to 5' probe sets ≤3.0; and (3) the multi-chip normalization scaling factor ≤10.0 Of the 223 samples, 213 (95.5%) passed the QC criteria. In a data set of 34 arrays (10 samples failing QC criteria, 24 randomly chosen samples passing QC criteria), correlation coefficients for pairwise comparisons of expression levels for 100 housekeeping genes in which at least one array failed the QC criteria were significantly lower (average Pearson r = 0.90 ± 0.04) and more broadly dispersed than correlation coefficients for pairwise comparisons between any two arrays that passed the QC criteria (average Pearson r = 0.97 ± 0.01). By using the QC cutoff criteria, the inter-array variability, as assessed by the coefficient of variation in the expression levels for 100 housekeeping genes, was reduced from 35.7% to 21.7%. Based on the aberrant housekeeping gene data generated from samples failing the established QC criteria, we propose that the QC criteria outlined in this study can accurately distinguish high quality from low quality data and can be used to delete poor quality microarray samples before proceeding to higher-order biological analyses and interpretation.

Project description:Next generation sequencing platforms have become essential tools for understanding DNA in a wide range of contexts. Their success heavily relies on the accuracy, sensitivity and specificity of methods used to discern differences between the reference genome and genomes under investigation. Here we compare the relative performances of five popular single nucleotide variant callers with and without their associated recommended hard filtering criteria. We compare: FreeBayes; the Genome Analysis Toolkit’s Haplotype Caller and Unified Genotyper; SAMtools; and VarScan. We tailor this comparison to suit smaller projects with modest sample numbers (n = 10) and coverage (~10X) to fill a current gap in the literature. Other comparison studies are generally applicable only to larger projects in model species, where there is access to large amounts of sequencing data and curated callsets for base and variant quality score recalibration. We estimated the accuracy, sensitivity and specificity of each pipeline according to the genotype concordance rate and number with the “truth” dataset for 10 canine samples. The truth dataset was defined as genotypes obtained from the CanineHD BeadChip array. Whole genome sequencing data was performed on the Illumina HiSeq2000 or HiSeq2500 platform as 100-101 base pair, paired end reads to an average sample coverage of 10.3X. Apart from GATK Haplotype Caller, applying recommended hard filters did not improve the performance of genotyping concordance at the tested levels of minimum coverage. The default VarScan pipeline with no additional filters applied (VarScan uses SAMtools mpileup, without base alignment quality computation) generally outperformed other callers in terms of accuracy, sensitivity and specificity. The results of this study demonstrate that hard filtering of variant calls from low-powered genome studies can impair accuracy, sensitivity and specificity of callsets and provides some benchmark performance metrics on a range of low coverage levels.

Project description:While the importance of random sequencing errors decreases at higher DNA or RNA sequencing depths, systematic sequencing errors (SSEs) dominate at high sequencing depths and can be difficult to distinguish from biological variants. These SSEs can cause base quality scores to underestimate the probability of error at certain genomic positions, resulting in false positive variant calls, particularly in mixtures such as samples with RNA editing, tumors, circulating tumor cells, bacteria, mitochondrial heteroplasmy, or pooled DNA. Most algorithms proposed for correction of SSEs require a training data set, which is typically either from a part of the data set being M-bM-^@M-^\recalibratedM-bM-^@M-^] (Genome Analysis ToolKit, or GATK) or from a separate data set with special characteristics (SysCall). Here, we combine the advantages of these approaches by adding synthetic RNA spike-in standards to human RNA, and use GATK to recalibrate base quality scores with reads mapped to the spike-in standards. Compared to conventional GATK recalibration that uses reads mapped to the genome, spike-ins improve the accuracy of Illumina base quality scores by a mean of 5 units, and by as much as 13 units M-BM- at CpG sites. In addition, since reads mapping to the genome are not used for recalibration, our method allows run-specific recalibration even for the many species without a comprehensive and accurate SNP database. We also use GATK with the spike-in standards to demonstrate that the Illumina RNA sequencing runs overestimate quality scores for AC, CC, GC, GG, and TC dinucleotides, while SOLiD has less dinucleotide SSEs but more SSEs for certain cycles. We conclude that using these DNA and RNA spike-in standards with GATK improves base quality score recalibration. Four human RNA samples with equimolar ERCC spike-in standards were sequenced on Illumina. Two human brain/liver/muscle RNA mixtures with dynamic range of ERCC spike-in standards were sequenced on SOLiD.

Dataset Information

Nègre2023 - Genome-Scale Metabolic Network of Penicillium rubens Wisconsin 54-1255

Dataset's files

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets