Genefinding.Pdf

Total Page:16

File Type:pdf, Size:1020Kb

Genefinding.Pdf Gene Finding GenBank Growth GenBank Growth • In 2003 • ~ 31 million sequences • ~ 37 billion base pairs GenBank: Exponential Growth Growth of GenBank in billions of base pairs from release 3 in April of 1994 to the current release, 142. Number of base pairs of sequence in GenBank release 142 for selected organisms Human Genome • Human genome = 3.4*109 bp • Number of genes = 30,000 - 100,000 • Gene percentage ~= 1% Gene Finding Genome Annotation kb 60 80 100 120 140 Single-exon gene Initial exon Internal exon Terminal exon Annotation amounts to finding: - content segments (exons - aa coding) - signal segments (Promoter, UTR, donor/acceptor site, Poly-Adenylation, etc.) The Central Dogma RNA polymerase performs transcription Ribosomes performs translation Transcription • Transcription starts a few bps before the coding region and ends a few bps after the coding region • Segments before and after the coding region are called untranslated regions (UTRs) • RNA polymerase starts transcription by binding to promoter regions before the transcription start sites • Promoters have signals which help cells control the expression of different genes Translation • Ribosomes produce amino acid sequences by translating information in the mRNA • Translation is based on the genetic code • One start codon • Three stop codons The Genetic Code Start codon Prokaryotic vs. Eukaryotic Genes Prokaryotes Eukaryotes small genomes large genomes high gene density low gene density no introns (or splicing) introns (splicing) no RNA processing RNA processing Eukaryotic and Prokaryotic Gene Structure Gene Finding in Prokaryotes • Most of the DNA sequence codes for proteins • 70% of the H.influenzea genome is coding • Each gene is one continuous stretch of base pairs (no introns) Prokaryotic DNA 5’ UTR 3’ UTR Coding Region Promoter Start codon Stop codon Transcription start site Open Reading Frames (ORFs) • Any sequence of nucleotides can be translated in three possible ways, depending on where the coding starts ACCUUAGCGUA Thr-Leu-Ala Pro-Stop-Arg Leu-Ser-Val • An open reading frame (ORF) is a reading frame with no stop codons Finding Long ORFs • Examine frequencies of stop codons to distinguish between coding and non-coding regions • Assuming a distribution on the codons, one can calculate the expected number a stop codon would be observed Finding Long ORFs • Scan the DNA sequence, looking for long ORFs in all three reading frames • Once a stop codon is detected, scan backwards to the start codon Shortcomings of the Algorithm • Fails to detect very short genes • Fails to detect overlapping long ORFs on opposite strands • #ORFs >> #genes (6500 ORFs in E.Coli, compared to 1100 genes only) Detecting Coding Regions • Consider the frequencies of each of the codons in coding regions • Amino acids Leu, Ala, and Try are coded by 6,4, and 1 different codons, respectively • Expected frequency ratio: 6:4:1 • However, in a protein, they appear with ratio 6.9:6.5:1 • Coding DNA is not random Using Markov Chains • Two Markov chains • G: “inside” a coding ORF • R: “inside” a non-coding ORF (NORF) • Given a sequence X, compute the ratio AG P (X|G) XiXi log = log +1 P (X|R) ! AR i XiXi+1 ORFs as Markov Chains • Given the set of all ORFs in a sequence • translate each ORF into a codon sequence • get 64-state Markov chain • Transition probabilities: probability of codon D being followed by codon E • Using this Markov chain, one can calculate the probability that a given ORF is a coding region Using Codon Frequencies Coding sequence a1, b1, c1, a2, b2, c2, . , an, bn, cn = × · · · × p1 fa1b1c1 fanbncn Prob. of three possible = × · · · × p2 fb c a fbncnan reading frames 1 1 2 +1 = × · · · × { p3 fc1a2b2 fcnan+1bn+1 Prob. reading frame #i being coding reading frame: pi Pi = p1 + p2 + p3 • Algorithm: • slide a window and compute the probabilities • Problems: • Dependence on codon frequency of already found genes • Detecting HGT and other heterogeneous regions Detecting Promoter Regions • Special sequences (with high frequency) in the promoter region • In E.Coli • TTGACA appears 35 bases before the transcription start point • TATAAT appears 12 bases before the transcription start point • Both strings appear with high frequency Another Approach • Define two models for “inside” and “outside” promoter regions • Calculate log-likelihood ratio P (X|Inside Promoter) log P (X|outside Promoter) Prokaryotic Gene Finding Tools Glimmer http://www.tigr.org/~salzberg/glimmer.html GeneMark http://opal.biology.gatech.edu/GeneMark/gmhmm2_prok.cgi Critica http://www.ttaxus.com/index.php?pagename=Software ORNL Annotation Pipeline http://compbio.ornl.gov/GP3/pro.shtml Gene Finding in Eukaryotes Eukaryotic DNA 5’ UTR 3’ UTR Introns Promoter Start codon Exons Stop codon Transcription PolyA Donor Acceptor start site Site Site Typical Figures: Vertebrate • Gene: ~30Kb • Coding region: ~1-2Kb • ~6 exons, each ~150bp • Promoter is ~6bp, and ~30bp before the transcription start site • Deviation from those figures is huge (dystrophin gene is ~2.4Mb, etc.) Distribution of Introns/Exons in Human Genome Sakharkar et al., In Silico Biology, 2004 Evolution of Gene Finding Tools 1982 extrinsic intrinsic hybrid Ab-initio Alignment-based Genie 1996 Genscan 1997 Comparative Genomics cDNA, Protein DNA Protein ExoFish Procrustes GenieEST GenieESTHOM 2000 1996 Informant HMM-based Rosetta Twinscan 2000 2001 Pair-HMM Phylo-HMM Slam Siepel-Haussler DoubleScan Jojic-Haussler 2002 2004 Performance Evaluation . Benchmark training set Burset-Guigo: • 570 vertebrate sequences • 6 accuracy measures • Sensitivity, specificity most common actual class positive negative At different levels: positive true positives false positives • Nucleotide (TP) (FP) predicted • Exon negative false negatives true negatives (FN) (TN) • Gene (unreliable) TP TP sensitivity = = actual positives TP + FN TP TP specificity = = predicted positives TP + FP A Solution: Markov Models • Usually a window of width 6, and a table of 4 6 observed frequencies for coding regions, and another for non-coding regions • (fifth other Markov model) • Homogeneous model: does not take into account any ORF information • Non-homogeneous model: different tables for the three possible reading frames • Problem when exons are too short Splicing and Splice Junctions • Splicing is the removal of introns by enzymes called spliceosomes • Spliceosomes contain both proteins and snRNA • snRNA recognizes the splice sites through RNA- RNA base-pairing • Splice sites recognition must be precise (a shift in the reading frame affects its message) • Many genes have alternative splicing: the exons used differ among different variations • Alternative splicing happens in more than 50% of the genes (on average, a gene has more than two variations) Consensus Eukaryotic Gene Sequence 100% frequency C-T rich AGGUAAGU.............CTGAC.........CAGG........ 15bp Donor Branch Acceptor site point site • This typical structure leads naturally to an algorithm based on position-specific weight matrices • However, this does not exploit all the information (reading frames, intron/exon states, etc.) • Not suitable for short genes An HMM Solution k P (exon of length k) = p (1 − p) Geometric dist. However, exon length does not seem to have a geometric distribution Solution • Generalized Hidden Markov Model (GHMM) GHMMs • The output of a state is a string of finite length • For a given state, the output string and its length are randomly chosen according to some probability distribution • Different states may have different probability distributions GHMMs • A finite set Q of states • Initial state probability distribution π T • Transition probabilities i j , for i, j ∈ Q • Length distribution f of the states ( f q is the length distribution of state q ) • Probabilistic models for each of the states, according to which output strings are generated upon visiting a state GenScan Model Burge&Karlin, 97 Gene Prediction Using the GenScan Model • A parse Φ of a sequence S of length L is a sequence of states ( q 1 , . , q t ) with an associated duration d i to each state, where t L = ! di i=1 • Intuitively, a parse is an annotation of a sequence, matching each region with a functional unit of a gene • Let Φ be a parse of sequence S • P ( S i |d i ) : probability of generating S i by the sequence generation model of state q i with length di The probability of generating S based on Φ is • t P (Φ, S) = πq1 fq1 (d1)P (S1|d1) ! Tqk−1qk fqk (dk)P (Sk|dk) • k=2 We have • P (Φ, S) P (Φ, S) P (Φ|S) = = P (S) !Φi P (Φi, S) Φi is a parse of length L GENSCAN • A computer program for gene identification • Uses the GHMM model described earlier • Uses a set of completely sequenced genes from GenBank for training C+G Content • The training set is divided into four categories depending on the C+G content of the sequence • For each of the categories, separate initial state probabilities, transition probabilities, and state length distributions are computed Initial State Probabilities • Initial state probabilities should be close to the frequencies with which various functional units occur in the actual genomic data • For example, if the estimated frequency of the non-coding intergenic region in 80%, then the initial probability for the state N (in the GenScan model) must be around 0.8 Transition Probabilities • Vary quite a bit with the C+G content • Computed separately for each of the categories • Transitions must be biologically permissible • E.g., transition from a P+ to an F+ state must have probability 1 State Length Distributions • Different functional units on a gene have different lengths • E.g., average exon length = 150bp; introns of length 1Kbp are common • Separate distributions for intron states in each category • Different distributions for initial/internal/terminal exons • Geometric distributions are used for UTR states Signal Models • Different signal models for different functional units • Weighted Matrix Model (WMM), where every position has its own independent distribution: PolyA, translation initiation and termination, promoters • Weighted Array Model (WAM), which allows for dependencies between adjacent positions: splice sites • Etc.
Recommended publications
  • Insights Into Comparative Genomics, Codon Usage Bias, And
    plants Article Insights into Comparative Genomics, Codon Usage Bias, and Phylogenetic Relationship of Species from Biebersteiniaceae and Nitrariaceae Based on Complete Chloroplast Genomes Xiaofeng Chi 1,2 , Faqi Zhang 1,2 , Qi Dong 1,* and Shilong Chen 1,2,* 1 Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810008, China; [email protected] (X.C.); [email protected] (F.Z.) 2 Qinghai Provincial Key Laboratory of Crop Molecular Breeding, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810008, China * Correspondence: [email protected] (Q.D.); [email protected] (S.C.) Received: 29 October 2020; Accepted: 17 November 2020; Published: 18 November 2020 Abstract: Biebersteiniaceae and Nitrariaceae, two small families, were classified in Sapindales recently. Taxonomic and phylogenetic relationships within Sapindales are still poorly resolved and controversial. In current study, we compared the chloroplast genomes of five species (Biebersteinia heterostemon, Peganum harmala, Nitraria roborowskii, Nitraria sibirica, and Nitraria tangutorum) from Biebersteiniaceae and Nitrariaceae. High similarity was detected in the gene order, content and orientation of the five chloroplast genomes; 13 highly variable regions were identified among the five species. An accelerated substitution rate was found in the protein-coding genes, especially clpP. The effective number of codons (ENC), parity rule 2 (PR2), and neutrality plots together revealed that the codon usage bias is affected by mutation and selection. The phylogenetic analysis strongly supported (Nitrariaceae (Biebersteiniaceae + The Rest)) relationships in Sapindales. Our findings can provide useful information for analyzing phylogeny and molecular evolution within Biebersteiniaceae and Nitrariaceae.
    [Show full text]
  • Mutation Bias Shapes Gene Evolution in Arabidopsis Thaliana ​
    bioRxiv preprint doi: https://doi.org/10.1101/2020.06.17.156752; this version posted June 18, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. Mutation bias shapes gene evolution in Arabidopsis thaliana ​ 1,2† 1 1 3,4 Monroe, J. Grey ,​ Srikant, Thanvi ,​ Carbonell-Bejerano, Pablo ,​ Exposito-Alonso, Moises ,​ 5​ ​ 6 7 ​ 1† ​ Weng, Mao-Lun ,​ Rutter, Matthew T. ,​ Fenster, Charles B. ,​ Weigel, Detlef ​ ​ ​ ​ 1 Department​ of Molecular Biology, Max Planck Institute for Developmental Biology, 72076 Tübingen, Germany 2 Department​ of Plant Sciences, University of California Davis, Davis, CA 95616, USA 3 Department​ of Plant Biology, Carnegie Institution for Science, Stanford, CA 94305, USA 4 Department​ of Biology, Stanford University, Stanford, CA 94305, USA 5 Department​ of Biology, Westfield State University, Westfield, MA 01086, USA 6 Department​ of Biology, College of Charleston, SC 29401, USA 7 Department​ of Biology and Microbiology, South Dakota State University, Brookings, SD 57007, USA † corresponding​ authors: [email protected], [email protected] ​ ​ ​ Classical evolutionary theory maintains that mutation rate variation between genes should be random with respect to fitness 1–4 and evolutionary optimization of genic 3,5 ​ mutation rates remains controversial .​ However, it has now become known that ​ cytogenetic (DNA sequence + epigenomic) features influence local mutation probabilities 6 ,​ which is predicted by more recent theory to be a prerequisite for beneficial mutation 7 rates between different classes of genes to readily evolve .​ To test this possibility, we ​ used de novo mutations in Arabidopsis thaliana to create a high resolution predictive ​ ​ ​ model of mutation rates as a function of cytogenetic features across the genome.
    [Show full text]
  • Chapter 3. the Beginnings of Genomic Biology – Molecular
    Chapter 3. The Beginnings of Genomic Biology – Molecular Genetics Contents 3. The beginnings of Genomic Biology – molecular genetics 3.1. DNA is the Genetic Material 3.6.5. Translation initiation, elongation, and termnation 3.2. Watson & Crick – The structure of DNA 3.6.6. Protein Sorting in Eukaryotes 3.3. Chromosome structure 3.7. Regulation of Eukaryotic Gene Expression 3.3.1. Prokaryotic chromosome structure 3.7.1. Transcriptional Control 3.3.2. Eukaryotic chromosome structure 3.7.2. Pre-mRNA Processing Control 3.3.3. Heterochromatin & Euchromatin 3.4. DNA Replication 3.7.3. mRNA Transport from the Nucleus 3.4.1. DNA replication is semiconservative 3.7.4. Translational Control 3.4.2. DNA polymerases 3.7.5. Protein Processing Control 3.4.3. Initiation of replication 3.7.6. Degradation of mRNA Control 3.4.4. DNA replication is semidiscontinuous 3.7.7. Protein Degradation Control 3.4.5. DNA replication in Eukaryotes. 3.8. Signaling and Signal Transduction 3.4.6. Replicating ends of chromosomes 3.8.1. Types of Cellular Signals 3.5. Transcription 3.8.2. Signal Recognition – Sensing the Environment 3.5.1. Cellular RNAs are transcribed from DNA 3.8.3. Signal transduction – Responding to the Environment 3.5.2. RNA polymerases catalyze transcription 3.5.3. Transcription in Prokaryotes 3.5.4. Transcription in Prokaryotes - Polycistronic mRNAs are produced from operons 3.5.5. Beyond Operons – Modification of expression in Prokaryotes 3.5.6. Transcriptions in Eukaryotes 3.5.7. Processing primary transcripts into mature mRNA 3.6. Translation 3.6.1.
    [Show full text]
  • "The" Genetic Code?
    Evolutionary Anthropology 14:6–11 (2005) CROTCHETS & QUIDDITIES “The” Genetic Code? KENNETH M. WEISS AND ANNE V. BUCHANAN The DNA-based code for protein through messenger and transfer RNA is widely themselves, that carry the informa- regarded as the code of life. But genomes are littered with other kinds of coding tion. elements as well, and all of them probably came after a supercode for the tRNA Your life depends on the fidelity of system itself. these many codes. Aberrant codes re- lated to cell behavior can lead to dys- genesis or various metabolic diseases. Evolution and the diversification of Everyone knows of “the” genetic Anomalous cell-surface proteins can organisms are made possible by code, by which nucleotide triplets in cause autoimmune destruction, and vi- codes, or arbitrary assignments of DNA in the nucleus of cells specify the ruses are the Alan Turings of life that “meaning,” in multiple ways. Many amino acid (aa) sequence of proteins. evolve ways to break their receptor are not widely appreciated. Codes al- This is the code described in text- codes to gain illicit entry into cells (Fig. low the same system of components books as the heart of the genetic the- 1). to be used for multiple purposes. ory of life and its evolution. Discover- But there is an additional code, a These can be open-ended, the way the ies in recent years have made things code of codes, that makes all of this alphabet and vocabulary make this more complicated by showing that ge- possible, including “the” genetic code column possible, but the flexibility of nomes are littered with all sorts of itself, and may be the oldest and most a code can become constrained once a other kinds of coding elements.
    [Show full text]
  • Designing Lentiviral Vectors for Gene Therapy of Genetic Diseases
    viruses Review Designing Lentiviral Vectors for Gene Therapy of Genetic Diseases Valentina Poletti 1,2,3,* and Fulvio Mavilio 4 1 Department of Woman and Child Health, University of Padua, 35128 Padua, Italy 2 Harvard Medical School, Harvard University, Boston, MA 02115, USA 3 Pediatric Research Institute City of Hope, 35128 Padua, Italy 4 Department of Life Sciences, University of Modena and Reggio Emilia, 41125 Modena, Italy; [email protected] * Correspondence: [email protected] Abstract: Lentiviral vectors are the most frequently used tool to stably transfer and express genes in the context of gene therapy for monogenic diseases. The vast majority of clinical applications involves an ex vivo modality whereby lentiviral vectors are used to transduce autologous somatic cells, ob- tained from patients and re-delivered to patients after transduction. Examples are hematopoietic stem cells used in gene therapy for hematological or neurometabolic diseases or T cells for immunotherapy of cancer. We review the design and use of lentiviral vectors in gene therapy of monogenic diseases, with a focus on controlling gene expression by transcriptional or post-transcriptional mechanisms in the context of vectors that have already entered a clinical development phase. Keywords: lentiviral vectors; transcriptional regulation; post-transcriptional regulation; miRNA; promoters; retroviral integration; ex vivo gene therapy Citation: Poletti, V.; Mavilio, F. 1. Introduction Designing Lentiviral Vectors for Gene Therapy of Genetic Diseases.
    [Show full text]
  • Analysis of Codon Usage Patterns in Giardia Duodenalis Based on Transcriptome Data from Giardiadb
    G C A T T A C G G C A T genes Article Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB Xin Li, Xiaocen Wang, Pengtao Gong, Nan Zhang, Xichen Zhang and Jianhua Li * Key Laboratory of Zoonosis Research, Ministry of Education, College of Veterinary Medicine, Jilin University, Changchun 130062, China; [email protected] (X.L.); [email protected] (X.W.); [email protected] (P.G.); [email protected] (N.Z.); [email protected] (X.Z.) * Correspondence: [email protected]; Tel.: +86-431-8783-6172; Fax: +86-431-8798-1351 Abstract: Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite- induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duo- denalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G.
    [Show full text]
  • Transcription and Open Reading Frame
    Transcription The expression of genetic information stored in the DNA sequence starts with synthesis of the RNA copy of the gene in a process called transcription. The RNA copy of the gene is called messenger RNA (mRNA). A special enzyme, RNA polymerase, recognizes a sequence, called promoter, on the DNA double helix upstream from the protein coding sequence. RNA polymerase binds to the promoter and opens it up: separates the complementary strands at about 12-nt-long region of the promoter. Then the enzyme starts the mRNA synthesis using all 4 NTPs. The DNA strand, which is used as a template by RNA polymerase, is called the template or antisense strand. The opposite strand of the gene, which sequence is identical to the sequence of mRNA (except the substitution T U, of course), is called the coding or sense strand. There is also a special sequence after the end of the gene, which signals to RNA polymerase to terminate the mRNA synthesis. Thus synthesized mRNA molecule, which includes the protein coding region flanked by short unrelated sequences from both sides, is either translated by the ribosome into the protein molecule at the spot (in case of prokaryotes) or is transported from the nucleus to the cytoplasm (in case of eukaryotes) and there it is translated by the ribosome. Open Reading Frame (ORF) Since the genetic code is triplet, three reading frames are possible for the same mRNA molecule. For instance, the sequence: …..AUUGCCUAACCCUUAGGG…. can be separated into triplets by three possible ways: ….AUUGCCUAACCCUUAGGG…. ….AUUGCCUAACCCUUAGGG…. ….AUUGCCUAACCCUUAGGG…. The frame, in which no stop codons are encountered, is called the Open Reading Frame (ORF).
    [Show full text]
  • Codon Usage Biases Co-Evolve with Transcription Termination Machinery
    RESEARCH ARTICLE Codon usage biases co-evolve with transcription termination machinery to suppress premature cleavage and polyadenylation Zhipeng Zhou1†, Yunkun Dang2,3†*, Mian Zhou4, Haiyan Yuan1, Yi Liu1* 1Department of Physiology, The University of Texas Southwestern Medical Center, Dallas, United States; 2State Key Laboratory for Conservation and Utilization of Bio- Resources in Yunnan, Yunnan University, Kunming, China; 3Center for Life Science, School of Life Sciences, Yunnan University, Kunming, China; 4State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, China Abstract Codon usage biases are found in all genomes and influence protein expression levels. The codon usage effect on protein expression was thought to be mainly due to its impact on translation. Here, we show that transcription termination is an important driving force for codon usage bias in eukaryotes. Using Neurospora crassa as a model organism, we demonstrated that introduction of rare codons results in premature transcription termination (PTT) within open reading frames and abolishment of full-length mRNA. PTT is a wide-spread phenomenon in Neurospora, and there is a strong negative correlation between codon usage bias and PTT events. Rare codons lead to the formation of putative poly(A) signals and PTT. A similar role for codon *For correspondence: usage bias was also observed in mouse cells. Together, these results suggest that codon usage [email protected] (YD); biases co-evolve with the transcription termination machinery to suppress premature termination of [email protected] (YL) transcription and thus allow for optimal gene expression. †These authors contributed DOI: https://doi.org/10.7554/eLife.33569.001 equally to this work Competing interests: The authors declare that no Introduction competing interests exist.
    [Show full text]
  • Small Open Reading Frames Tiny Treasures of the Non-Coding Genomic Regions
    GENERAL ARTICLE Small Open Reading Frames Tiny Treasures of the Non-coding Genomic Regions A Yazhini Open Reading Frames (ORFs) are the DNA sequences in the genome that has the potential to be translated. Generally, only long ORFs (≥ 300 nucleotides or nt) are thought to be protein coding regions and are considered as genes in the genome annotation pipeline. Until recent years, small ORFs (smORFs) of less than 100 codons (< 300 nt) were regarded as non-functional on the basis of empirical observations. How- A Yazhini is a research ever, recent work on ribosome profiling and mass spectrom- student at Molecular etry have led to the discovery of many translating functional Biophysics Unit, Indian Institute of Science. She small ORFs and presence of their stable peptide products. works mainly on protein Further, examples of biologically active peptides with vital evolution, protein structure regulatory functions underline the importance of smORFs prediction and structure of in cell functions. Genome-wide analysis shows that smORFs macromolecule assemblies under the guidance of Prof. N are conserved across diverse species, and the functional char- Srinivasan. Overall, her acterization of their peptides reveals their critical role in a research revolves around broad spectrum of regulatory mechanisms. Further analy- structural and mechanistic sis of small ORFs is likely throw light on many exciting, un- understanding of proteins. explored regulatory mechanisms in different developmental stages and tissue types. 1. Introduction Gene is the protein encoding functional unit in the genome. It consists of promoter, protein coding (exon) and terminal regions. The protein coding regions of the gene are called the ‘open read- ing frames’ (ORFs).
    [Show full text]
  • Unconventional Viral Gene Expression Mechanisms As Therapeutic Targets
    Review Unconventional viral gene expression mechanisms as therapeutic targets https://doi.org/10.1038/s41586-021-03511-5 Jessica Sook Yuin Ho1,3, Zeyu Zhu1,3 & Ivan Marazzi1,2 ✉ Received: 8 June 2020 Accepted: 22 March 2021 Unlike the human genome that comprises mostly noncoding and regulatory sequences, Published online: 19 May 2021 viruses have evolved under the constraints of maintaining a small genome size while expanding the efciency of their coding and regulatory sequences. As a result, viruses Check for updates use strategies of transcription and translation in which one or more of the steps in the conventional gene–protein production line are altered. These alternative strategies of viral gene expression (also known as gene recoding) can be uniquely brought about by dedicated viral enzymes or by co-opting host factors (known as host dependencies). Targeting these unique enzymatic activities and host factors exposes vulnerabilities of a virus and provides a paradigm for the design of novel antiviral therapies. In this Review, we describe the types and mechanisms of unconventional gene and protein expression in viruses, and provide a perspective on how future basic mechanistic work could inform translational eforts that are aimed at viral eradication. Expression of a gene in the human genome is a multistep and heavily (for example, alternative splicing) or use unique strategies. Here we regulated process that resembles a production line. Protein-coding describe the diverse ways by which viral genomes give rise to genes and genes are transcribed almost exclusively by RNA polymerase II (RNAPII). proteins that deviate from the canonical framework of human genes, During transcription, quality-control checkpoints are implemented to restricting our analyses to eukaryotes and their viruses.
    [Show full text]
  • Recitation 8 Solutions (PDF)
    SOLUTION KEY- SECTION 8 Questions: 1. Assume that the sequence of DNA below is a short protein-encoding gene; the sequences in between the transcription start and stop sites are shown. The entire DNA sequence of the very short gene is: 5' CGCTTATAGAACCCAATCTCTCATAGGC 3' 3' GCGAATATCTTGGGTTAGAGAGTATCCG 5' iv. What would the resulting mRNA be if the top strand of this DNA molecule were used as a template in transcription? Label the 5’ and 3’ ends of your molecule. 5’GCCUAUGAGAGAUUGGGUUCUAUAAGCG3’ v. What is the full sequence of the protein that would be translated from this RNA? Label the N and C termini of your molecule. N-met-arg-asp-trp-val-leu-C vi. What would happen to the encoded protein if the underlined nucleotide C were mutated to a T? N-met-arg-asp-C (This creates a premature stop codon resulting in a truncated protein). 2. Drawn below is part of a wild-type gene. The DNA sequence shown encodes the last amino acids of a protein that is normally 380 amino acids long. The bold & underlined codon indicates the correct reading frame of this gene. The lower strand of the gene is used as the template during the transcription of mRNA from this gene. …GCTAAGTATTGCTCAAGATTAGGATGATAAATAACTGG–3’ …CGATTCATAACGAGTTCTAATCCTACTATTTATTGACC–5’ iii. In the copy of the sequence drawn below, circle one base pair that you could change to make a mutant form of the gene that produces a protein that is now 381 amino acids long. Indicate the identity of one new base pair that could take its place.
    [Show full text]
  • Mutations in Noncoding Regions of GJB1 Are a Major Cause of X-Linked CMT
    Mutations in noncoding regions of GJB1 are a major cause of X-linked CMT Pedro J. Tomaselli, MD, ABSTRACT * MSc Objective: To determine the prevalence and clinical and genetic characteristics of patients with X- Alexander M. Rossor, linked Charcot-Marie-Tooth disease (CMT) due to mutations in noncoding regions of the gap junc- * MRCP, PhD tion b-1 gene (GJB1). Alejandro Horga, MD, Methods: Mutations were identified by bidirectional Sanger sequence analysis of the 595 bases PhD of the upstream promoter region, and 25 bases of the 39 untranslated region (UTR) sequence in Zane Jaunmuktane, patients in whom mutations in the coding region had been excluded. Clinical and neurophysiologic FRCPath data were retrospectively collected. Aisling Carr, MRCP, PhD Results: Five mutations were detected in 25 individuals from 10 kindreds representing 11.4% of Paola Saveri, BSc all cases of CMTX1 diagnosed in our neurogenetics laboratory between 1996 and 2016. Four . 1 . Giuseppe Piscosquito, pathogenic mutations, c.-17G A, c.-17 1G T, c.-103C T, and c.-146-90_146-89insT were 9 * . 9 GJB1 MD detected in the 5 UTR. A novel mutation, c. 15C T, was detected in the 3 UTR of in 2 9 Davide Pareyson, MD unrelated families with CMTX1 and is the first pathogenic mutation in the 3 UTR of any myelin- Matilde Laura, MD, PhD associated CMT gene. Mutations segregated with the phenotype, were at sites predicted to be Julian C. Blake, FRCP pathogenic, and were not present in the normal population. Roy Poh, PhD Conclusions: Mutations in noncoding DNA are a major cause of CMTX1 and highlight the impor- James Polke, PhD tance of mutations in noncoding DNA in human disease.
    [Show full text]