Designing minimal genomes using whole-cell models.
Ontology highlight
ABSTRACT:
SUBMITTER: Rees-Garbutt J
PROVIDER: S-EPMC7012841 | biostudies-literature | 2020 Feb
REPOSITORIES: biostudies-literature
Publications
Designing minimal genomes using whole-cell models.
Rees-Garbutt Joshua J Chalkley Oliver O Landon Sophie S Purcell Oliver O Marucci Lucia L Grierson Claire C
Nature communications 20200211 1
In the future, entire genomes tailored to specific functions and environments could be designed using computational tools. However, computational tools for genome design are currently scarce. Here we present algorithms that enable the use of design-simulate-test cycles for genome design, using genome minimisation as a proof-of-concept. Minimal genomes are ideal for this purpose as they have a simple functional assay whether the cell replicates or not. We used the first (and currently only publis ...[more]
Similar Datasets
Project description:Accurate reference genome sequences provide the foundation for modern molecular biology and genomics as the interpretation of sequence data to study evolution, gene expression, and epigenetics depends heavily on the quality of the genome assembly used for its alignment. Correctly organising sequenced fragments such as contigs and scaffolds in relation to each other is a critical and often challenging step in the construction of robust genome references. We previously identified misoriented regions in the mouse and human reference assemblies using Strand-seq, a single cell sequencing technique that preserves DNA directionality Here we demonstrate the ability of Strand-seq to build and correct full-length chromosomes by identifying which scaffolds belong to the same chromosome and determining their correct order and orientation, without the need for overlapping sequences. We demonstrate that Strand-seq exquisitely maps assembly fragments into large related groups and chromosome-sized clusters without using new assembly data. Using template strand inheritance as a bi-allelic marker, we employ genetic mapping principles to cluster scaffolds that are derived from the same chromosome and order them within the chromosome based solely on directionality of DNA strand inheritance. We prove the utility of our approach by generating improved genome assemblies for several model organisms including the ferret, pig, Xenopus, zebrafish, Tasmanian devil and the Guinea pig.
| S-EPMC8037727 | biostudies-literature
Project description:Ethnic differences in pharmacogenomic (PGx) variants have been well documented in literature and could significantly impact variability in response and adverse events to therapeutics. India is a large country with diverse ethnic populations of distinct genetic architecture. India's national genome sequencing initiative (IndiGen) provides a unique opportunity to explore the landscape of PGx variants using population-scale whole genome sequences. We have analyzed the IndiGen variation dataset (N = 1029 genomes) along with global population scale databases to map the most prevalent clinically actionable and potentially deleterious PGx variants among Indians. Differential frequencies for the known and novel variants were studied and interaction of the disrupted PGx genes affecting drug responses were analyzed by performing a pathway analysis. We have highlighted significant differences in the allele frequencies of clinically actionable PGx variants in Indians when compared to the global populations. We identified 134 mostly common (allele frequency [AF] > 0.1) potentially deleterious PGx variants that could alter or inhibit the function of 102 pharmacogenes in Indians. We also estimate that on, an average, each Indian individual carried eight PGx variants (single nucleotide variants) that have a direct impact on the choice of treatment or drug dosing. We have also highlighted clinically actionable PGx variants and genes for which preemptive genotyping is most recommended for the Indian population. The study has put forward the most comprehensive PGx landscape of the Indian population from whole genomes that could enable optimized drug selection and genotype-guided prescriptions for improved therapeutic outcomes and minimizing adverse events.
| S-EPMC9010271 | biostudies-literature
Project description:Cancer is one of the leading causes of death, but mortality can be reduced by detecting tumors earlier so that treatment is initiated at a less aggressive stage. The tradeoff between costs associated with screening and its benefit makes the decision of whom to screen and when a challenge. To enable comparisons across screening strategies for any cancer type, we demonstrate a mathematical modeling platform based on the theory of queuing networks designed for quantifying the benefits of screening strategies. Our methodology can be used to design optimal screening protocols and to estimate their benefits for specific patient populations. Our method is amenable to exact analysis, thus circumventing the need for simulations, and is capable of exactly quantifying outcomes given variability in the age of diagnosis, rate of progression, and screening sensitivity and intervention outcomes. We demonstrate the power of this methodology by applying it to data from the Surveillance, Epidemiology and End Results (SEER) program. Our approach estimates the benefits that various novel screening programs would confer to different patient populations, thus enabling us to formulate an optimal screening allocation and quantify its potential effects for any cancer type and intervention.
| S-EPMC9182689 | biostudies-literature
Project description:Two common sources of DNA for whole exome sequencing (WES) are whole blood (WB) and immortalized lymphoblastoid cell line (LCL). However, it is possible that LCLs have a substantially higher rate of mutation than WB, causing concern for their use in sequencing studies. We compared results from paired WB and LCL DNA samples for 16 subjects, using LCLs of low passage number (<5). Using a standard analysis pipeline we detected a large number of discordant genotype calls (approximately 50 per subject) that we segregated into categories of "confidence" based on read-level quality metrics. From these categories and validation by Sanger sequencing, we estimate that the vast majority of the candidate differences were false positives and that our categories were effective in predicting valid sequence differences, including LCLs with putative mosaicism for the non-reference allele (3-4 per exome). These results validate the use of DNA from LCLs of low passage number for exome sequencing.
| S-EPMC3812417 | biostudies-literature
Project description:Predictive models that accurately emulate complex scientific processes can achieve speed-ups over numerical simulators or experiments and at the same time provide surrogates for improving the subsequent analysis. Consequently, there is a recent surge in utilizing modern machine learning methods to build data-driven emulators. In this work, we study an often overlooked, yet important, problem of choosing loss functions while designing such emulators. Popular choices such as the mean squared error or the mean absolute error are based on a symmetric noise assumption and can be unsuitable for heterogeneous data or asymmetric noise distributions. We propose Learn-by-Calibrating, a novel deep learning approach based on interval calibration for designing emulators that can effectively recover the inherent noise structure without any explicit priors. Using a large suite of use-cases, we demonstrate the efficacy of our approach in providing high-quality emulators, when compared to widely-adopted loss function choices, even in small-data regimes.
| S-EPMC7648787 | biostudies-literature
Project description:UnlabelledBackgroundWith the progress of modern sequencing technologies a large number of complete genomes are now available. Traditionally the comparison of two related genomes is carried out by sequence alignment. There are cases where these techniques cannot be applied, for example if two genomes do not share the same set of genes, or if they are not alignable to each other due to low sequence similarity, rearrangements and inversions, or more specifically to their lengths when the organisms belong to different species. For these cases the comparison of complete genomes can be carried out only with ad hoc methods that are usually called alignment-free methods.MethodsIn this paper we propose a distance function based on subword compositions called Underlying Approach (UA). We prove that the matching statistics, a popular concept in the field of string algorithms able to capture the statistics of common words between two sequences, can be derived from a small set of "independent" subwords, namely the irredundant common subwords. We define a distance-like measure based on these subwords, such that each region of genomes contributes only once, thus avoiding to count shared subwords a multiple number of times. In a nutshell, this filter discards subwords occurring in regions covered by other more significant subwords.ResultsThe Underlying Approach (UA) builds a scoring function based on this set of patterns, called underlying. We prove that this set is by construction linear in the size of input, without overlaps, and can be efficiently constructed. Results show the validity of our method in the reconstruction of phylogenetic trees, where the Underlying Approach outperforms the current state of the art methods. Moreover, we show that the accuracy of UA is achieved with a very small number of subwords, which in some cases carry meaningful biological information.Availabilityhttp://www.dei.unipd.it/∼ciompin/main/underlying.html.
| S-EPMC3549825 | biostudies-literature
Project description:Industrial ecology has revolutionized our understanding of material stocks and flows in our economy and society. For this important discipline to have even deeper impact, we must understand the inherent nature of these materials in terms of human health and the environment. This paper focuses on methods to design synthetic chemicals to reduce their intrinsic ability to cause adverse consequence to the biosphere. Advances in the fields of computational chemistry and molecular toxicology in recent decades allow the development of predictive models that inform the design of molecules with reduced potential to be toxic to humans or the environment. The approach presented herein builds on the important work in quantitative structure-activity relationships by linking toxicological and chemical mechanistic insights to the identification of critical physical-chemical properties needed to be modified. This in silico approach yields design guidelines using boundary values for physiochemical properties. Acute aquatic toxicity serves as a model endpoint in this study. Defining value ranges for properties related to bioavailability and reactivity eliminates 99% of the chemicals in the highest concern for acute aquatic toxicity category. This approach and its future implementations are expected to yield very powerful tools for life cycle assessment practitioners and molecular designers that allow rapid assessment of multiple environmental and human health endpoints and inform modifications to minimize hazard.
| S-EPMC4443366 | biostudies-literature
Project description:BackgroundThe regulatory information encoded in the DNA of promoter regions usually enforces a minimal, non-zero distance between the coding regions of neighboring genes. However, the size of this minimal regulatory space is not generally known. In particular, it is unclear if minimal promoter size differs between species and between uni- and bi-directionally acting regulatory regions.ResultsAnalyzing the genomes of 11 yeasts, we show that the lower size limit on promoter-containing regions is species-specific within a relatively narrow range (80-255 bp). This size limit applies equally to regions that initiate transcription on one or both strands, indicating that bi-directional promoters and uni-directional promoters are constrained similarly. We further find that young, species-specific regions are on average much longer than older regions, suggesting either a bias towards deletions or selection for genome compactness in yeasts. While the length evolution of promoter-less intergenic regions is well described by a simplistic, purely neutral model, regions containing promoters typically show an excess of unusually long regions. Regions flanked by divergently transcribed genes have a bi-modal length distribution, with short lengths found preferentially among older regions. These old, short regions likely harbor evolutionarily conserved bi-directionally active promoters. Surprisingly, some of the evolutionarily youngest regions in two of the eleven species (S. cerevisiae and K. waltii) are shorter than the lower limit observed in older regions.ConclusionsThe minimal chromosomal space required for transcriptional regulation appears to be relatively similar across yeast species, and is the same for uni-directional and bi-directional promoters. New intergenic regions created by genome rearrangements tend to evolve towards the more narrow size distribution found among older regions.
| S-EPMC3128071 | biostudies-literature
Project description:The hippocampus is a heavily studied brain structure due to its involvement in learning and memory. Detailed models of excitatory, pyramidal cells in hippocampus have been developed using a range of experimental data. These models have been used to help us understand, for example, the effects of synaptic integration and voltage gated channel densities and distributions on cellular responses. However, these cellular outputs need to be considered from the perspective of the networks in which they are embedded. Using modeling approaches, if cellular representations are too detailed, it quickly becomes computationally unwieldy to explore large network simulations. Thus, simple models are preferable, but at the same time they need to have a clear, experimental basis so as to allow physiologically based understandings to emerge. In this article, we describe the development of simple models of CA1 pyramidal cells, as derived in a well-defined experimental context of an intact, whole hippocampus preparation expressing population oscillations. These models are based on the intrinsic properties and frequency-current profiles of CA1 pyramidal cells, and can be used to build, fully examine, and analyze large networks.
| S-EPMC4215760 | biostudies-literature
Project description:Importance:To date, there is no diagnostic test for Kawasaki disease (KD). Diagnosis is based on clinical features shared with other febrile conditions, frequently resulting in delayed or missed treatment and an increased risk of coronary artery aneurysms. Objective:To identify a whole-blood gene expression signature that distinguishes children with KD in the first week of illness from other febrile conditions. Design, Setting, and Participants:The case-control study comprised a discovery group that included a training and test set and a validation group of children with KD or comparator febrile illness. The setting was pediatric centers in the United Kingdom, Spain, the Netherlands, and the United States. The training and test discovery group comprised 404 children with infectious and inflammatory conditions (78 KD, 84 other inflammatory diseases, and 242 bacterial or viral infections) and 55 healthy controls. The independent validation group comprised 102 patients with KD, including 72 in the first 7 days of illness, and 130 febrile controls. The study dates were March 1, 2009, to November 14, 2013, and data analysis took place from January 1, 2015, to December 31, 2017. Main Outcomes and Measures:Whole-blood gene expression was evaluated using microarrays, and minimal transcript sets distinguishing KD were identified using a novel variable selection method (parallel regularized regression model search). The ability of transcript signatures (implemented as disease risk scores) to discriminate KD cases from controls was assessed by area under the curve (AUC), sensitivity, and specificity at the optimal cut point according to the Youden index. Results:Among 404 patients in the discovery set, there were 78 with KD (median age, 27 months; 55.1% male) and 326 febrile controls (median age, 37 months; 56.4% male). Among 202 patients in the validation set, there were 72 with KD (median age, 34 months; 62.5% male) and 130 febrile controls (median age, 17 months; 56.9% male). A 13-transcript signature identified in the discovery training set distinguished KD from other infectious and inflammatory conditions in the discovery test set, with AUC of 96.2% (95% CI, 92.5%-99.9%), sensitivity of 81.7% (95% CI, 60.0%-94.8%), and specificity of 92.1% (95% CI, 84.0%-97.0%). In the validation set, the signature distinguished KD from febrile controls, with AUC of 94.6% (95% CI, 91.3%-98.0%), sensitivity of 85.9% (95% CI, 76.8%-92.6%), and specificity of 89.1% (95% CI, 83.0%-93.7%). The signature was applied to clinically defined categories of definite, highly probable, and possible KD, resulting in AUCs of 98.1% (95% CI, 94.5%-100%), 96.3% (95% CI, 93.3%-99.4%), and 70.0% (95% CI, 53.4%-86.6%), respectively, mirroring certainty of clinical diagnosis. Conclusions and Relevance:In this study, a 13-transcript blood gene expression signature distinguished KD from other febrile conditions. Diagnostic accuracy increased with certainty of clinical diagnosis. A test incorporating the 13-transcript disease risk score may enable earlier diagnosis and treatment of KD and reduce inappropriate treatment in those with other diagnoses.
| S-EPMC6233768 | biostudies-literature