Deterministic evolution and stringent selection during preneoplasia.

The earliest events during human tumour initiation, although poorly characterized, may hold clues to malignancy detection and prevention1. Here we model occult preneoplasia by biallelic inactivation of TP53, a common early event in gastric cancer, in human gastric organoids. Causal relationships between this initiating genetic lesion and resulting phenotypes were established using experimental evolution in multiple clonally derived cultures over 2 years. TP53 loss elicited progressive aneuploidy, including copy number alterations and structural variants prevalent in gastric cancers, with evident preferred orders. Longitudinal single-cell sequencing of TP53-deficient gastric organoids similarly indicates progression towards malignant transcriptional programmes. Moreover, high-throughput lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programmes repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a marked degree of phenotypic convergence in premalignant epithelial organoids. These data imply predictability in the earliest stages of tumorigenesis and show evolutionary constraints and barriers to malignant transformation, with implications for earlier detection and interception of aggressive, genome-instable tumours.

Modeling RNA-Binding Protein Specificity In Vivo by Precisely Registering Protein-RNA Crosslink Sites.

RNA-binding proteins (RBPs) regulate post-transcriptional gene expression by recognizing short and degenerate sequence motifs in their target transcripts, but precisely defining their binding specificity remains challenging. Crosslinking and immunoprecipitation (CLIP) allows for mapping of the exact protein-RNA crosslink sites, which frequently reside at specific positions in RBP motifs at single-nucleotide resolution. Here, we have developed a computational method, named mCross, to jointly model RBP binding specificity while precisely registering the crosslinking position in motif sites. We applied mCross to 112 RBPs using ENCODE eCLIP data and validated the reliability of the discovered motifs by genome-wide analysis of allelic binding sites. Our analyses revealed that the prototypical SR protein SRSF1 recognizes clusters of GGA half-sites in addition to its canonical GGAGGA motif. Therefore, SRSF1 regulates splicing of a much larger repertoire of transcripts than previously appreciated, including HNRNPD and HNRNPDL, which are involved in multivalent protein assemblies and phase separation.

JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profiles.

JASPAR (https://jaspar.elixir.no/) is a widely-used open-access database presenting manually curated high-quality and non-redundant DNA-binding profiles for transcription factors (TFs) across taxa. In this 10th release and 20th-anniversary update, the CORE collection has expanded with 329 new profiles. We updated three existing profiles and provided orthogonal support for 72 profiles from the previous release's UNVALIDATED collection. Altogether, the JASPAR 2024 update provides a 20% increase in CORE profiles from the previous release. A trimming algorithm enhanced profiles by removing low information content flanking base pairs, which were likely uninformative (within the capacity of the PFM models) for TFBS predictions and modelling TF-DNA interactions. This release includes enhanced metadata, featuring a refined classification for plant TFs' structural DNA-binding domains. The new JASPAR collections prompt updates to the genomic tracks of predicted TF binding sites (TFBSs) in 8 organisms, with human and mouse tracks available as native tracks in the UCSC Genome browser. All data are available through the JASPAR web interface and programmatically through its API and the updated Bioconductor and pyJASPAR packages. Finally, a new TFBS extraction tool enables users to retrieve predicted JASPAR TFBSs intersecting their genomic regions of interest.

JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

Improving the productivity of Xinjiang cotton in heat-limited regions under two life history strategies.

Cotton is a major cash crop globally, playing a pivotal role in the textile sector. However, cotton growers in Xinjiang region are experiencing cotton yield penalty caused by limited heat environment. In this region, limited heat conditions strongly arrest cotton plant growth and development resulting in recued productivity. To counteract this problem, there is an urgent need to robustly identify efficient management strategies to improve plant performance and increase cotton yield under heat-limited situations. This will hold crucial implications for agricultural sustainability and global cotton supply. This review article identified challenges faced by cotton producers under heat limited environments with potential solutions to enhance cotton productivity. Specifically, we focused on the implementation of two life history strategies including planting early maturing and cold tolerant cultivars, and adjusting sowing date that can promote early maturity and increase cold stress tolerance. These strategies have shown promising results in protecting cotton plants from limited heat injury and consequently improved cotton productivity. By focusing on Xinjiang province unique climate and associated agronomic practices, valuable insights can be gained, which may have broader applications in other heat-limited cotton-growing regions globally. This comprehensive review endeavors to provide a foundation for future research and practical interventions aimed at boosting cotton yields under limited heat areas.

Genome wide identification of GDSL gene family explores a novel GhirGDSL26 gene enhancing drought stress tolerance in cotton.

BACKGROUND: Current climate change scenarios are posing greater threats to the growth and development of plants. Thus, significant efforts are required that can mitigate the negative effects of drought on the cotton plant. GDSL esterase/lipases can offer an imperative role in plant development and stress tolerance. However, thesystematic and functional roles of the GDSL gene family, particularly in cotton under water deficit conditions have not yet been explored. RESULTS: In this study, 103, 103, 99, 198, 203, 239, 249, and 215 GDSL proteins were identified in eight cotton genomes i.e., Gossypium herbaceum (A1), Gossypium arboretum (A2), Gossypium raimondii (D5), Gossypium hirsutum (AD1), Gossypium barbadense (AD2), Gossypium tomentosum (AD3), Gossypium mustelinum (AD4), Gossypium darwinii (AD5), respectively. A total of 198 GDSL genes of Gossypium hirsutum were divided into eleven clades using phylogenetic analysis, and the number of GhirGDSL varied among different clades. The cis-elements analysis showed that GhirGDSL gene expression was mainly related to light, plant hormones, and variable tense environments. Combining the results of transcriptome and RT-qPCR, GhirGDSL26 (Gh_A01G1774), a highly up-regulated gene, was selected for further elucidating its tole in drought stress tolerance via estimating physiological and biochemical parameters. Heterologous expression of the GhirGDSL26 gene in Arabidopsis thaliana resulted in a higher germination and survival rates, longer root lengths, lower ion leakage and induced stress-responsive genes expression under drought stress. This further highlighted that overexpressed plants had a better drought tolerance as compared to the wildtype plants. Moreover, 3, 3'-diaminobenzidine (DAB) and Trypan staining results indicated reduced oxidative damage, less cell membrane damage, and lower ion leakage in overexpressed plants as compared to wild type. Silencing of GhirGDSL26 in cotton via VIGS resulting in a susceptible phenotype, higher MDA and H2O2 contents, lower SOD activity, and proline content. CONCLUSION: Our results demonstrated that GhirGDSL26 plays a critical role in cotton drought stress tolerance. Current findings enrich our knowledge of GDSL genes in cotton and provide theoretical guidance and excellent gene resources for improving drought tolerance in cotton.

Rationale and design of the TUXEDO-2 India study: Ultra-Thin strUt Supraflex Cruz versus XiencE in a Diabetic pOpulation with multi-vessel disease-2.

BACKGROUND: The role of percutaneous coronary interventions (PCI) in patients with diabetes mellitus and multi-vessel disease has been questioned by the results of the FREEDOM trial, which showed superiority of coronary artery bypass graft(CABG) over first generation drug-eluting stents (DES) including a reduction in mortality. In the light of safer and more efficacious stents and significantly better medical management, those results that date back to 2012 need to be revisited. TUXEDO-2 is a study designed to compare two contemporary stents in Indian diabetic patients with multi-vessel disease. AIMS: The primary objective of the TUXEDO-2 study is to compare the clinical outcomes of PCI with ultra-thin Supraflex Cruz vs Xience when combined with contemporary optimal medical therapy (OMT) in diabetic patients with multi-vessel disease. The secondary objective is to compare clinical outcomes between a pooled cohort from both arms of the study (Supraflex Cruz + Xience; PCI arm) vs CABG based on a performance goal derived from the CABG arm of the FREEDOM trial (historical cohort). The tertiary objective is a randomized comparison of ticagrelor vs prasugrel in addition to aspirin for the composite of ischemic and bleeding events. METHODS: In this prospective, open-label, multi-centre, 2 × 2 factorial, randomized, controlled study, 1,800 patients with diabetes mellitus and multi-vessel disease (inclusion criteria similar to FREEDOM trial) with indication for coronary revascularization will be randomly assigned to Supraflex Cruz or Xience stents and also to ticagrelor- or prasugrel- based antiplatelet strategies. All patients will receive guideline directed OMT and optimal PCI including image- and physiology-guided complete revascularization where feasible. The patients will be followed through five years to assess their clinical status and major clinical events. The primary endpoint is a non-inferiority comparison of target lesion failure at one-year for Supraflex Cruz vs Xience (primary objective) with an expected event rate of 11% and a non-inferiority margin of 4.5%. For PCI vs CABG (secondary objective), the primary endpoint is major adverse cardiac events (MACE), defined as a composite of all cause death, nonfatal myocardial infarction, or stroke at one-year and yearly up to five years, with a performance goal of 21.6%. For ticagrelor vs prasugrel (tertiary objective), the primary endpoint is composite of death, myocardial infarction, stroke, and major bleeding as per the Bleeding Academic Research Consortium (BARC) at one-year with expected event rate of 15% and a non-inferiority margin of 5%. CONCLUSIONS: The TUXEDO-2 study is a contemporary study involving state-of-the-art PCI combined with guideline directed OMT in a complex subset of patients with diabetes mellitus and multi-vessel disease. The trial will answer the question as to whether a biodegradable polymer coated ultra-thin Supraflex Cruz stent is an attractive option for PCI in diabetic patients with multi-vessel disease. It will also help address the question whether the results of FREEDOM trial would have been different in the current era of safer and more efficacious stents and modern medical therapy. In addition, the comparative efficacy and safety of ticagrelor vs prasugrel in addition to aspirin will be evaluated. (CTRI/2019/11/022088).

Mitochondrial gene expression analysis reveals aberrant transcription of cox3 in Gossypium barbadense CMS line H276A.

Cotton cytoplasmic male sterility (CMS) has been extensively studied; however, information regarding its molecular mechanisms has not yet been disclosed. Therefore, to explore the molecular mechanism of pollen abortion of cotton CMS line H276A, transcript profiles of 30 mitochondrial protein-encoding genes at tetrad stage were conducted with northern blot and a differential expression gene cox3 was identified. Quantitative reverse-transcribed PCR (qRT-PCR) analysis indicated that the expression level of cox3 in the CMS line H276A was 0.39-fold compared to its maintainer line H276B. In addition, the immunoblot analysis revealed that the amount of COX3 protein was decreased to 59.38% in CMS line H276A. The 5` and 3` terminals of the transcript of cox3 in two materials were determined simultaneously with circularized RNA reverse-transcribed PCR (CR-RT-PCR). The data indicated that seven 5` end of transcript of cox3 in H276A (-451/-464/-465/-467/-471/-472/-508 respect to ATG) were identified which were different from that of H276B (-411/-412). A total of 15 single nucleotide polymorphisms (SNPs) was detected by clone sequencing analysis of upstream of cox3. To our knowledge, we are the first to comprehensively analyze the transcripts of the mitochondrial genome in the cotton CMS line and to identify the 5` and 3` terminals of the transcript of cox3 in cotton. Our data will provide a framework for a better understanding of molecular mechanisms of CMS and mitochondrial gene expression in cotton.

BiasAway: command-line and web server to generate nucleotide composition-matched DNA background sequences.

MOTIVATION: Accurate motif enrichment analyses depend on the choice of background DNA sequences used, which should ideally match the sequence composition of the foreground sequences. It is important to avoid false positive enrichment due to sequence biases in the genome, such as GC-bias. Therefore, relying on an appropriate set of background sequences is crucial for enrichment analysis. RESULTS: We developed BiasAway, a command line tool and its dedicated easy-to-use web server to generate synthetic sequences matching any k-mer nucleotide composition or select genomic DNA sequences matching the mononucleotide composition of the foreground sequences through four different models. For genomic sequences, we provide precomputed partitions of genomes from nine species with five different bin sizes to generate appropriate genomic background sequences. AVAILABILITY AND IMPLEMENTATION: BiasAway source code is freely available from Bitbucket (https://bitbucket.org/CBGR/biasaway) and can be easily installed using bioconda or pip. The web server is available at https://biasaway.uio.no and a detailed documentation is available at https://biasaway.readthedocs.io. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

JASPAR 2020: update of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net) is an open-access database of curated, non-redundant transcription factor (TF)-binding profiles stored as position frequency matrices (PFMs) for TFs across multiple species in six taxonomic groups. In this 8th release of JASPAR, the CORE collection has been expanded with 245 new PFMs (169 for vertebrates, 42 for plants, 17 for nematodes, 10 for insects, and 7 for fungi), and 156 PFMs were updated (125 for vertebrates, 28 for plants and 3 for insects). These new profiles represent an 18% expansion compared to the previous release. JASPAR 2020 comes with a novel collection of unvalidated TF-binding profiles for which our curators did not find orthogonal supporting evidence in the literature. This collection has a dedicated web form to engage the community in the curation of unvalidated TF-binding profiles. Moreover, we created a Q&A forum to ease the communication between the user community and JASPAR curators. Finally, we updated the genomic tracks, inference tool, and TF-binding profile similarity clusters. All the data is available through the JASPAR website, its associated RESTful API, and through the JASPAR2020 R/Bioconductor package.

Unmet Need for Further LDL-C Lowering in India despite Statin Therapy: Lipid Association of India Recommendations for the Use of Bempedoic Acid.

Lipid-lowering therapy plays a crucial role in reducing adverse cardiovascular (CV) events in patients with established atherosclerotic cardiovascular disease (ASCVD) and familial hypercholesterolemia. Lifestyle interventions along with high-intensity statin therapy are the first-line management strategy followed by ezetimibe. Only about 20-30% of patients who are on maximally tolerated statins reach recommended low-density lipoprotein cholesterol (LDL-C) goals. Several factors contribute to the problem, including adherence issues, prescription of less than high-intensity statin therapy, and de-escalation of statin dosages, but in patients with very high baseline LDL-C levels, including those with familial hypercholesterolemia and those who are intolerant to statins, it is critical to expand our arsenal of LDL-C-lowering medications. Moreover, in the extreme risk group of patients with an LDL-C goal of ≤30 mg/dL according to the Lipid Association of India (LAI) risk stratification algorithm, there is a significant residual risk requiring the addition of non-statin drugs to achieve LAI recommended targets. This makes bempedoic acid a welcome addition to the existing non-statin therapies such as ezetimibe, bile acid sequestrants, and PCSK9 inhibitors. A low frequency of muscle-related side effects, minimal drug interactions, a significant reduction in high-sensitivity C-reactive protein (hsCRP), and a lower incidence of new-onset or worsening diabetes make it a useful adjunct for LDL-C lowering. However, the CV outcomes trial results are still pending. In this LAI consensus document, we discuss the pharmacology, indications, contraindications, advantages, and evidence-based recommendations for the use of bempedoic acid in clinical practice.

UniBind: maps of high-confidence direct TF-DNA interactions across nine species.

BACKGROUND: Transcription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. It is critical to locate these TF-DNA interactions to understand transcriptional regulation. Efforts to predict bona fide TFBSs benefit from the availability of experimental data mapping DNA binding regions of TFs (chromatin immunoprecipitation followed by sequencing - ChIP-seq). RESULTS: In this study, we processed ~ 10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~ 56 million TFBSs with experimental and computational support for direct TF-DNA interactions for 644 TFs in > 1000 cell lines and tissues. These TFBSs were used to predict > 197,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface ( https://unibind.uio.no/ ), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions. CONCLUSIONS: UniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species.

Analysis of codon usage bias and evolution in the chloroplast genome of Mesona chinensis Benth.

Mesona chinensis Benth (MCB) is one of the main economic crops in tropical and subtropical areas. To understand the codon usage bias (CUB) in M. chinensis Benth, chloroplast genome is essential to study its genetic law, molecular phylogenetic relationships, and exogenous gene expression. Results showed that the GC content of 53 CDS sequences was 37.95%, and GC1, GC2, and GC3 content were 46.02%, 38.26%, and 29.85%, respectively. The general GC content order was GC1>GC2>GC3. Moreover, the majority of genes had an effective number of codon (ENC) value greater than 40, except ndhE, rps8, and rps18. Correlation analysis results revealed that the GC content was significantly correlated with GC1, GC2, GC3, and ENC. Neutrality plot analysis, ENC-plot analysis, and PR2-plot analysis presented that the CUB of M. chinensis Benth chloroplast genome was mainly affected by mutation and selection. In addition, GGG, GCA, and TCC were found to be the optimal codons. Furthermore, results of cluster analysis and evolutionary tree showed that M. chinensis Benth was closely related to Ocimum basilicum, indicating that there was a certain correlation between the CUB of the chloroplast gene and the genetic relationship of plant species. Overall, the study on the CUB of chloroplast genome laid a basis for genetic modification and phylogenetic research of M. chinensis Benth chloroplast genome.

Nucleo-cytoplasmic interactions affect the 5' terminal transcription of mitochondrial genes between the isonuclear CMS line UG93A and its maintainer line UG93B of kenaf (Hibiscus cannabinus).

Gene expression and translation in plant mitochondria remain poorly understood due to the complicated transcription of its mRNA. In this study, we report the 5' and 3' RNA extremities and promoters of five mitochondrial genes, atp1, atp4, atp6, atp9, and cox3. The results reveal that four genes (atp1, atp4, atp6, and cox3) are transcribed from multiple initiation sites but with a uniform transcript at the 3' end, indicating that heterogeneity of the 5' end is a common feature in the transcription of kenaf mitochondrial genes. Furthermore, we found that the transcription initiation sites of these four genes are significantly different in UG93A, UG93B, and the F1 hybrid. These data indicate that nuclear loci and unknown transcription factors within the mitochondria of different cytoplasmic types may be involved in mitochondrial transcription. Promoter architecture analysis showed that the promoter core sequences are conserved in the kenaf mitochondrial genome but are highly divergent, suggesting that these elements are essential for the promoter activity of mitochondrial genes in kenaf. Our results reveal that the heterogeneity of the 5' end and uniformity at the 3' end are common transcriptional features of mitochondrial genes. These data provide essential information for understanding the transcription of mitochondrial genes in kenaf and can be used as a reference for other plants.

Identification of differentially expressed genes and pathways in isonuclear kenaf genotypes under salt stress.

Salinity is a potential abiotic stress and globally threatens crop productivity. However, the molecular mechanisms underlying salt stress tolerance with respect to cytoplasmic effect, gene expression, and metabolism pathway under salt stress have not yet been reported in isonuclear kenaf genotypes. To fill this knowledge gap, growth, physiological, biochemical, transcriptome, and cytoplasm changes in kenaf cytoplasmic male sterile (CMS) line (P3A) and its iso-nuclear maintainer line (P3B) under 200 mM sodium chloride (NaCl) stress and control conditions were analyzed. Salt stress significantly reduced leaf soluble protein, soluble sugars, proline, chlorophyll content, antioxidant enzymatic activity, and induced oxidative damage in terms of higher MDA contents in both genotypes. The reduction of these parameters resulted in a reduced plant growth compared with control. However, P3A was relatively more tolerant to salt stress than P3B. This tolerance of P3A was further confirmed by improved physio-biochemical traits under salt stress conditions. Transcriptome analysis showed that 4256 differentially expressed genes (DEGs) between the two genotypes under salt stress were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that photosynthesis, photosynthesis antenna-protein, and plant hormone signal transduction pathways might be associated with the improved NaCl stress tolerance in P3A. Conclusively, P3A cytoplasmic male sterile could be a potential salt-tolerant material for future breeding program of kenaf and can be used for phytoremediation of salt-affected soils. These data provide a valuable resource on the cytoplasmic effect of salt-responsive genes in kenaf and salt stress tolerance in kenaf.

Degradation of mitochondrial structure and deficiency of complex I were associated with the transgenic CMS of rice.

BACKGROUND: Mitochondria play a significant role in plant cytoplasmic male sterility (CMS). In our previous study, mitochondrial complex I genes, nad4, nad5, and nad7 showed polymorphisms between the transgenic CMS line M2BS and its wild type M2B. The sterility mechanism of the M2BS at cytological, physiological, biochemical, and molecular level is not clear. RESULTS: Cytological observation showed that the anthers were light yellow, fissured, invalid in KI-I2, and full of irregularly typical abortion pollen grains in M2BS. Transmission electron microscopic (TEM) observation revealed no nucleus and degraded mitochondria with obscure cristae in anther cells of M2BS. The results of staining for H2O2 presented a large number of electron dense precipitates (edp) in intercellular space of anther cells of M2BS at anthesis. Moreover, the anther respiration rate and complex I activity of M2BS were significantly lower than those of wild type M2B during pollen development. Furthermore, RNA editing results showed only nad7 presented partially edited at 534th nucleotides. The expression of nad5 and nad7 revealed significant differences between M2B and M2BS. CONCLUSIONS: Our data demonstrated that mitochondrial structural degradation and complex I deficiency might be associated with transgenic CMS of rice.

A map of direct TF-DNA interactions in the human genome.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the most popular assay to identify genomic regions, called ChIP-seq peaks, that are bound in vivo by transcription factors (TFs). These regions are derived from direct TF-DNA interactions, indirect binding of the TF to the DNA (through a co-binding partner), nonspecific binding to the DNA, and noise/bias/artifacts. Delineating the bona fide direct TF-DNA interactions within the ChIP-seq peaks remains challenging. We developed a dedicated software, ChIP-eat, that combines computational TF binding models and ChIP-seq peaks to automatically predict direct TF-DNA interactions. Our work culminated with predicted interactions covering >4% of the human genome, obtained by uniformly processing 1983 ChIP-seq peak data sets from the ReMap database for 232 unique TFs. The predictions were a posteriori assessed using protein binding microarray and ChIP-exo data, and were predominantly found in high quality ChIP-seq peaks. The set of predicted direct TF-DNA interactions suggested that high-occupancy target regions are likely not derived from direct binding of the TFs to the DNA. Our predictions derived co-binding TFs supported by protein-protein interaction data and defined cis-regulatory modules enriched for disease- and trait-associated SNPs. We provide this collection of direct TF-DNA interactions and cis-regulatory modules through the UniBind web-interface (http://unibind.uio.no).

A fractional order Covid-19 epidemic model with Mittag-Leffler kernel.

In this article, we are studying fractional-order COVID-19 model for the analytical and computational aspects. The model consists of five compartments including; ` ` S c ″ which denotes susceptible class, ` ` E c ″ represents exposed population, ` ` I c ″ is the class for infected people who have been developed with COVID-19 and can cause spread in the population. The recovered class is denoted by ` ` R c ″ and ` ` V c ″ is the concentration of COVID-19 virus in the area. The computational study shows us that the spread will be continued for long time and the recovery reduces the infection rate. The numerical scheme is based on the Lagrange's interpolation polynomial and the numerical results for the suggested model are similar to the integer order which gives us the applicability of the numerical scheme and effectiveness of the fractional order derivative.