A model of DNA sequence evolution

Statistical studies of gene populations on the purine/pyrimidine alphabet have shown that the mean occurrence probability of thei-motif YRY(N) _i YRY (R=purine, Y=pyrimidine, N=R or Y) is not uniform by varyingi in the range [1,99], but presents a maximum ati=6 in the following populations: protein coding genes of eukaryotes, prokaryotes, chloroplasts and mitrochondria, and also viral introns, ribosomal RNA genes and transfer RNA genes (Arquès and Michel, 1987b,J. theor. Biol. 128, 457–461). From the “universality” of this observation, we suggested that the oligonucleotide YRY(N)₆ is a primitive one and that it has a central function in DNA sequence evolution (Arquès and Michel, 1987b,J. theor. Biol. 128, 457–461). Following this idea, we introduce a concept of a model of DNA sequence evolution which will be validated according to a shema presented in three parts. In the first part, using the last version of the gene database, the YRY(N)₆YRY preferential occurrence (maximum ati=6) is confirmed for the populations mentioned above and is extended to some newly analysed populations: chloroplast introns, chloroplast 5′ regions, mitochondrial 5′ regions and small nuclear RNA genes. On the other hand, the YRY(N)₆YRY preferential occurrence and periodicities are used in order to classify 18 gene populations. In the second part, we will demonstrate that several statistical features characterizing different gene populations (in particular the YRY(N)₆YRY preferential occurrence and the periodicities) can be retrieved from a simple Markov model based on the mixing of the two oligonucleotides YRY(N)₆ and YRY(N)₃ and based on the percentages of RYR and YRY in the unspecified trinucleotides (N)₃ of YRY(N)₆ and YRY(N)₃. Several properties are identified and prove in particular that the oligonucleotide mixing is an independent process and that several different features are functions of a unique parameter. In the third part, the return of the model to the reality shows a strong correlation between reality and simulation concerning the presence of large alternating purine/pyrimidine stretches and of periodicities. It also contributes to a greater understanding of biological reality, e.g. the presence or the absence of large alternating purine/pyrimidine stretches can be explained as being a simple consequence of the mixing of two particular oligonucleotides. Finally, we believe that such an approach is the first step toward a unified model of DNA sequence evolution allowing the molecular understanding of both the origin of life and the actual biological reality.     

A simulation of the genetic periodicities modulo 2 and 3 with processes of nucleotide insertions and deletions.

Recently, a new genetic process termed RNA editing has been identified showing insertions and deletions of nucleotides in particular RNA molecules. On the other hand, there are a few non-random statistical properties in genes: in particular, the periodicity modulo 3 (P3) associated with an open reading frame, the periodicity modulo 2 (P2) associated with alternating purine/pyrimidine stretches, the YRY(N)6YRY preferential occurrence (R = purine = adenine or guanine, Y = pyrimidine = cytosine or thymine, N = R or Y) representing a "code" of the DNA helix pitch, etc. The problem investigated here is whether a process of the type RNA editing can lead to the non-random statistical properties commonly observed in genes. This paper will show in particular that: The process of insertions and deletions of mononucleotides in the initial sequence [YRY(N)3]* [series of YRY(N)3] can lead to the periodicity modulo 2 (P2). The process of insertions and deletions of trinucleotides in the initial sequence [YRY(N)6]* [series of YRY(N)6] can lead to the periodicity modulo 3 (P3) and the YRY(N)6YRY preferential occurrence. Furthermore, these two processes lead to a strong correlation with the reality, namely the mononucleotide insertion/deletion process, with the 5' eukaryotic regions and the trinucleotide insertion/deletion process, with the eukaryotic protein coding genes.     

Monocyte chemoattractant protein-induced protein 1 directly degrades viral miRNAs with a specific motif and inhibits KSHV infection

Kaposi''s sarcoma-associated herpesvirus (KSHV) expresses miRNAs during latency. However, regulation of viral miRNAs remains largely unknown. Our prior studies demonstrated that MCPIP1 regulates KSHV miRNA biogenesis by degrading most KSHV pre-miRNAs through its RNase activity. Some viral pre-miRNAs are partially resistant to degradation by MCPIP1. Here, we further characterized MCPIP1 substrate specificity and its antiviral potential against KSHV infection. In vitro cleavage assays and binding assays showed that MCPIP1 cleavage efficiency is related to binding affinity. Motif-based sequence analysis identified that KSHV pre-miRNAs that are well degraded by MCPIP1 have a 5-base motif (M5 base motif) within their terminal loops and this motif region consists of multiple pyrimidine-purine-pyrimidine (YRY) motifs. We further demonstrated that mutation of this M5 base motif within terminal loop of pre-miRNAs inhibited MCPIP1-mediated RNA degradation. We also revealed that MCPIP1 has an antiviral effect against KSHV infection. MCPIP1 can reduce the expression of Dicer, which in turn restricts KSHV infection. Conclusively, our findings demonstrated that MCPIP1 inhibited KSHV infection and suppressed viral miRNA biogenesis by directly degrading KSHV pre-miRNAs and altering the expression of miRNA biogenesis factors.     

Structural architecture of an RNA that competitively inhibits RNase L

Activation of RNase L endonuclease activity is part of the mammalian innate immune response to viral infection. The poliovirus RNA genome contains a sequence in its protein-coding region that can act as a competitive inhibitor of RNase L. Mutation, sequence, and functional analysis of this competitive inhibitor RNA (ciRNA) revealed that its activity depends on specific sequences, showed that a loop-loop hairpin interaction forms in the ciRNA, and suggested the presence of a loop E motif. These features lead to the hypothesis that the ciRNA's function is conferred in part by a specific three-dimensional folded RNA architecture. By using a combination of biophysical, mutational, and functional studies, we have mapped features of the three-dimensional architecture of the ciRNA in its unbound form. We show that the loop-loop interaction forms in the free ciRNA and affects the overall structure, perhaps forming long-range tertiary interactions with the loop E motif. Local tight RNA-RNA backbone packing occurs in parts of the structure, but the fold appears to be less stable than many other tightly packed RNAs. This feature may allow the ciRNA to accommodate the translocation of ribosomes and polymerase across this multifunctional region of the viral RNA but also to function as an RNase L inhibitor.     

Transduplication resulted in the incorporation of two protein-coding sequences into the <Emphasis Type="Italic">Turmoil</Emphasis>-1 transposable element of <Emphasis Type="Italic">C. elegans</Emphasis>

   

Transposable elements may acquire unrelated gene fragments into their sequences in a process called transduplication. Transduplication of protein-coding genes is common in plants, but is unknown of in animals. Here, we report that the Turmoil-1 transposable element in C. elegans has incorporated two protein-coding sequences into its inverted terminal repeat (ITR) sequences. The ITRs of Turmoil-1 contain a conserved RNA recognition motif (RRM) that originated from the rsp-2 gene and a fragment from the protein-coding region of the cpg-3 gene. We further report that an open reading frame specific to C. elegans may have been created as a result of a Turmoil-1 insertion. Mutations at the 5' splice site of this open reading frame may have reactivated the transduplicated RRM motif.     

How to find the real one (at the level of pre-mRNA splicing)

The mature mRNA always carries nucleotide sequences that faithfully mirror the protein product according to the niles of the genetic code. However, in the chromosome, the nucleotide sequence that represents a certain protein is interrupted by additional sequences. Therefore, most eukaryotic genes are longer than their final mRNA products. The human genome project revealed that only a tiny portion of sequences serves as protein-coding region and almost one quarter of the genome is occupied by non-coding intervening sequences. The elimination of these non-coding regions from the precursor RNA in a process termed splicing must be extremely precise, because even a single nucleotide mistake may cause a fatal error. At present, two types of intervening sequences have been identified in protein-coding genes. One of them, the U2-dependent or major-class is prevalent and represents 99% of known sequences. The other one, the so-called U12-dependent or minor-class of introns, occurs in much lesser amounts in the genome. The basic problem of nuclear splicing concerns i/ the molecular mechanisms, which ensure that the coding regions are correctly recognized and spliced together: ii/ the principles and mechanisms that guarantee the high fidelity of the splicing system; iii/ the differences in the excision mechanisms of the two classes of introns. We are going to present models explaining how intervening sequences are accurately removed and the coding regions correctly juxtaposed. The two splicing mechanisms will also be compared.     

Complete mitochondrial DNA sequence of the Japanese spiny lobster,Panulirus japonicus (Crustacea: Decapoda)

We determined the complete nucleotide sequence of the mitochondrial genome for a Japanese spiny lobster, Panulirus japonicus (Crustacea: Decapoda). The entire genome was amplified using long polymerase chain reaction, and the products were subsequently used as templates for direct sequencing using a primer-walking strategy. The genome (15,717 base pairs) contained the same 37 genes (two ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes) plus the putative control region as found in other arthropods, with the gene order identical to that of typical arthropods. Preliminary phylogenetic analyses of selected arthropods using concatenated amino acid sequences of the 13 protein-coding genes strongly supported monophyly of Decapoda species and confidently rejected "Macroura", a conventional taxon that shares an elongated abdominal body.     

Complete mitochondrial DNA sequence of Aulopus japonicus (Teleostei: Aulopiformes), a basal Eurypterygii: longer DNA sequences and higher-level relationships

For the first step toward resolution of the higher-level relationships of the order Aulopiformes (Teleostei: Eurypterygii) using longer DNA sequences, we determined the complete mitochondrial DNA sequence for Aulopus japonicus (Aulopodidae). The entire genome was purified by gene amplification using a long PCR technique, and the products were subsequently used as templates for PCR with 63 fish-versatile and 3 species-specific primers that amplify contiguous, overlapping segments of the entire genome. Direct sequencing of the PCR products demonstrated that the genome (16 653 base pairs [bp]) contained the same 37 mitochondrial genes (2 ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes) as found in other vertebrates, with the gene order identical to that in typical vertebrates. Maximum-parsimony analysis using nucleotide sequences from the concatenated 12 protein-coding genes (no third codon positions and excluding the ND6 gene) plus 22 tRNA genes (stem regions only) from eight teleosts placed A. japonicus in a reasonable phylogenetic position; those from individual protein-coding genes and the concatenated 22 tRNA genes alone, however, did not reproduce the expected phylogeny with few exceptions, probably owing to insufficient phylogenetic information in these smaller data sets. This result suggests that further taxonomic sampling and sequencing efforts may clarify limits and intra- and interrelationships of this morphologically and ecologically diverse group of fishes using mitochondrial genomic (mitogenomic) data. Received: August 31, 2000 / Revised: December 20, 2000 / Accepted: January 23, 2001     

Intermolecular interaction between a branching ribozyme and associated homing endonuclease mRNA

RNA tertiary interactions involving docking of GNRA (N; any base; R; purine) hairpin loops into helical stem structures on other regions of the same RNA are one of the most common RNA tertiary interactions. In this study, we investigated a tertiary association between a GAAA hairpin tetraloop in a small branching ribozyme (DiGIR1) and a receptor motif (HEG P1 motif) present in a hairpin structure on a separate mRNA molecule. DiGIR1 generates a 2', 5' lariat cap at the 5' end of its downstream homing endonuclease mRNA by catalysing a self-cleavage branching reaction at an internal processing site. Upon release, the 5' end of the mRNA forms a distinct hairpin structure termed HEG P1. Our biochemical data, in concert with molecular 3D modelling, provide experimental support for an intermolecular tetraloop receptor interaction between the L9 GAAA in DiGIR1 and a GNRA tetraloop receptor-like motif (UCUAAG-CAAGA) found within the HEG P1. The biological role of this interaction appears to be linked to the homing endonuclease expression by promoting post-cleavage release of the lariat capped mRNA. These findings add to our understanding of how protein-coding genes embedded in nuclear ribosomal DNA are expressed in eukaryotes and controlled by ribozymes.     

The complete mitochondrial DNA sequence of the Atlantic salmon, Salmo salar.

The complete sequence of the Atlantic salmon (Salmo salar) mitochondrial genome has been determined. The entire sequence is 16665 base pairs (bp) in length, with a gene content (13 protein-coding, two ribosomal RNA [rRNA] and 22 transfer RNA [tRNA] genes) and order conforming to that observed in most other vertebrates. Base composition and codon usage have been detailed. Nucleotide and derived amino acid sequences of the 13 protein-coding genes from Atlantic salmon have been compared with their counterparts in rainbow trout. A putative structure for the origin of L-strand replication (O(L)) is proposed, and sequence features of the control region (D-loop) are described.     

本氏烟Ⅰ型启动子的克隆及其转录起始位点分析

启动子是转录水平上一个重要的调控元件,其决定着基因的表达模式和表达强度。Ⅰ型启动子具有高转录活性和种属间特异性等特点。如将其应用于植物RNA病毒载体表达系统,有利于提高表达系统的表达效率和生物安全性。本氏烟(Nicotiana benthaminana)是一种被广泛地应用于植物生物反应器和植物病理学的模式生物,但是现有核酸数据库中尚没有其Ⅰ型启动子的相关信息。因此,克隆本氏烟Ⅰ型启动子并分析其转录起始位点就具有重要的应用价值。通过半巢式PCR获得了514 bp的本氏烟Ⅰ型启动子序列(KC352713);生物信息学分析初步预测其转录起始位点位于其核心序列TATA(G)TA(N)GGGGG中的第3位A处;通过植物RNA病毒表达载体和5'RACE技术在体内验证本氏烟Ⅰ型启动子转录起始位点与生物信息学预测结果一致。研究结果为深入研究Ⅰ型启动子和构建Ⅰ型启动子介导转录的植物RNA病毒载体表达系统奠定了基础。     

A Group I Intron in the Nuclear Small Subunit rRNA Gene of Cryptendoxyla hypophloia, an Ascomycetous Fungus: Evidence for a New Major Class of Group I Introns

The ascomycetous fungus Cryptendoxyla hypophloia contains an insertion of 433 base pairs in the genes encoding nuclear small subunit ribosomal RNA. Secondary structure analyses of the insert reveal characteristics indicative of a Group I intron, including elements P, Q, R, and S; however, the sequences of these conserved regions deviate significantly from recognized consensus sequences for Group I introns. Principal-components analysis, based on 79 nucleotide positions from the conserved core sequences of 93 Group I introns, identified 17 introns similar to that of C. hypophloia. This grouping, which includes inserts from phylogenetically diverse organisms, cannot readily be classified in any previously recognized major group of Group I introns. We propose the creation of a new group, IE, to accommodate these sequences, and discuss the evolutionary relationships between group IE and other major groups of Group I introns. Received: 11 January 1998 / Accepted: 12 October 1998