Secondary Structure Model Research Articles

Characterization of putative proteins encoded by variable ORFs in white spot syndrome virus genome

BackgroundWhite Spot Syndrome Virus (WSSV) is an enveloped double-stranded DNA virus which causes mortality of several species of shrimp, being considered one of the main pathogens that affects global shrimp farming. This virus presents a complex genome of ~ 300 kb and viral isolates that present genomes with great identity. Despite this conservation, some variable regions in the WSSV genome occur in coding regions, and these putative proteins may have some relationship with viral adaptation and virulence mechanisms. Until now, the functions of these proteins were little studied. In this work, sequences and putative proteins encoded by WSSV variable regions were characterized in silico.ResultsThe in silico approach enabled determining the variability of some sequences, as well as the identification of some domains resembling the Formin homology 2, RNA recognition motif, Xeroderma pigmentosum group D repair helicase, Hemagglutinin and Ankyrin motif. The information obtained from the sequences and the analysis of secondary and tertiary structure models allow to infer that some of these proteins possibly have functions related to protein modulation/degradation, intracellular transport, recombination and endosome fusion events.ConclusionsThe bioinformatics approaches were efficient in generating three-dimensional models and to identify domains, thereby enabling to propose possible functions for the putative polypeptides produced by the ORFs wsv129, wsv178, wsv249, wsv463a, wsv477, wsv479, wsv492, and wsv497.

A new saxicolous Circinaria species (Megasporaceae) from northeast China

In an investigation of Circinaria species in volcanic mountains of northeast China, one new species, C. laxilobata was identified based on morphological characteristics and DNA sequence data. The new species differs morphologically from other members of the C. contorta group by its greenish olive, thin thallus with dispersed areoles and elongated, diverging, branching marginal areoles distinctly delimited from the inner part, round apothecia with entire thalline margin, moniliform paraphyses, asci with 1–2 ellipsoid spores, and the absence of secondary metabolites. Comparisons of the ITS sequences of the nuclear ribosomal DNA repeat tandem, both in phylogenetic analyses and secondary structure models of ITS2, support the taxonomic distinctness of this species. The new species is described and illustrated in detail. A key to the species of Circinaria from China is also provided.

Length and secondary structure of the 5′ non-coding regions of mouse p53 mRNA transcripts - mouse as a model organism for p53 gene expression studies

ABSTRACTTranscription initiation sites of Trp53 gene in mice were determined using the 5′RACE method. Based on sequence alignment of the 5′-terminal regions of p53 mRNA in mammals, the site for the most abundant transcript turned out to be essentially identical with that determined for human TP53 gene and slightly differed for the longest transcripts, in mice and humans. Secondary structures of the 5′ -terminal regions of the shorter, most abundant and the longest mouse transcripts were determined in vitro and the shorter transcript was also mapped in transfected mouse cells. For the first time, secondary structure models of the 5′ terminus of two mouse p53 mRNAs were proposed. Comparing these models with the conservativeness of the nucleotide sequence of the 5′-terminal region of mRNA in mouse and other mammals, the possible function of the selected structural domains of this region was discussed. To elucidate the translation mechanisms, the two studied mRNAs were translated in the presence of an increasing concentration of the cap analog. For the longest transcript, the data suggested that IRES element(s) was/were involved in translation initiation. Additionally, changes in p53 synthesis under genotoxic and endoplasmic reticulum stress conditions in mouse cells were analyzed.

Signal Transduction in Human Cell Lysate via Dynamic RNA Nanotechnology.

Dynamic RNA nanotechnology with small conditional RNAs (scRNAs) offers a promising conceptual approach to introducing synthetic regulatory links into endogenous biological circuits. Here, we use human cell lysate containing functional Dicer and RNases as a testbed for engineering scRNAs for conditional RNA interference (RNAi). scRNAs perform signal transduction via conditional shape change: detection of a subsequence of mRNA input X triggers formation of a Dicer substrate that is processed to yield small interfering RNA (siRNA) output anti-Y targeting independent mRNA Y for destruction. Automated sequence design is performed using the reaction pathway designer within NUPACK to encode this conditional hybridization cascade into the scRNA sequence subject to the sequence constraints imposed by X and Y. Because it is difficult for secondary structure models to predict which subsequences of mRNA input X will be accessible for detection, here we develop the RNAhyb method to experimentally determine accessible windows within the mRNA that are provided to the designer as sequence constraints. We demonstrate the programmability of scRNA regulators by engineering scRNAs for transducing in both directions between two full-length mRNAs X and Y, corresponding to either the forward molecular logic “if X then not Y” (X Y) or the reverse molecular logic “if Y then not X” (Y X). In human cell lysate, we observe a strong OFF/ON conditional response with low crosstalk, corresponding to a ≈20-fold increase in production of the siRNA output in response to the cognate versus noncognate full-length mRNA input. 2′OMe-RNA chemical modifications protect signal transduction reactants and intermediates against RNase degradation while enabling Dicer processing of signal transduction products. Because diverse biological pathways interact with RNA, scRNAs that transduce between detection of endogenous RNA inputs and production of biologically active RNA outputs hold great promise as a synthetic regulatory paradigm.

Peltigera neodegenii sp. nov. from Central China

Peltigera neodegenii from Central China is described and illustrated as a new species. It is similar in general appearance to P. degenii, P. membranacea, and P. canina, but is distinguished by its shiny upper surface in the central part of lobes, tomenta around the margins of upper surface of lobes, distinct and raised veins on the lower surface of lobes, and simple rhizines. Comparisons of the ITS (ITS1, 5.8S, ITS2) sequences of the nuclear ribosomal DNA repeat tandem, both in phylogenetic analysis and secondary structure models of ITS2, support the taxonomic distinctness of this species.

A consensus secondary structure of ITS2 for the diatom Order Cymatosirales (Mediophyceae, Bacillariophyta) and reappraisal of the order based on DNA, morphology, and reproduction

In 1983, Hasle and colleagues removed cymatosiroid diatoms from the pennates, and erected a new centric diatom family, the Cymatosiraceae, mainly to accommodate for their newly discovered mode of sexual reproduction. The new family consisted of two subfamilies differing in frustule structure. The family was later elevated to the rank of Order Cymatosirales Round and Crawford. We revisited intra-ordinal relationships within Cymatosirales using combined genetic (DNA sequences), morphological (valve and frustule structure), and reproductive (auxospore type) characters. In total, 36 cymatosiroid strains from 19 species representing 13 genera (80% of the total number of extant genera; nine of them represented by their generitypes) were used in this study. Instead of only the commonly used loci (18S rRNA and plastidal genes) to infer diatom phylogeny, we developed a consensus secondary structure model of the Internal Transcribed Spacer 2 (ITS2) for this order and applied it to aid in sequence alignments for ITS2. This improved the alignment and thus the robustness of the phylogenetic framework. The compensatory base changes (CBCs) found in ITS2 secondary structures were mapped onto the multi-gene (18S rRNA + ITS2 + rbcL) phylogenetic tree topology. In all these trees, all species grouped into two morphologically and genetically distinct clades. Each clade was supported by multiple CBCs, as did all the clades representing genera. However, these clades did not correspond to the previously established subfamilies. Consequently, we amend the Order Cymatosirales and family Cymatosiraceae, and propose a new family, the Leyanellaceae. The structure of the auxospore was an additional synapomorphic character for Cymatosirales. Overall, we demonstrate a novel approach to study diatom phylogeny across a broader taxonomic range using ITS2 secondary structural information. Our results suggest that this approach might be useful in establishing higher taxonomic relationships in other groups of diatoms.

Form confers function: Case of the 3'X region of the hepatitis C virus genome.

At the 3’ end of genomic hepatitis C virus (HCV) RNA there is a highly conserved untranslated region, the 3’X-tail, which forms part of the 3’UTR. This region plays key functions in regulation of critical processes of the viral life cycle. The 3’X region is essential for viral replication and infectivity. It is also responsible for regulation of switching between translation and transcription of the viral RNA. There is some evidence indicating the contribution of the 3’X region to the translation efficiency of the viral polyprotein and to the encapsidation process. Several different secondary structure models of the 3’X region, based on computer predictions and experimental structure probing, have been proposed. It is likely that the 3’X region adopts more than one structural form in infected cells and that a specific equilibrium between the various forms regulates several aspects of the viral life cycle. The most intriguing explanations of the structural heterogeneity problem of the 3’X region came with the discovery of its involvement in long-range RNA-RNA interactions and the potential for homodimer formation. This article summarizes current knowledge on the structure and function of the 3’X region of hepatitis C genomic RNA, reviews previous opinions, presents new hypotheses and summarizes the questions that still remain unanswered.

SHAPE probing pictures Mg2+-dependent folding of small self-cleaving ribozymes.

Self-cleaving ribozymes are biologically relevant RNA molecules which catalyze site-specific cleavage of the phosphodiester backbone. Gathering knowledge of their three-dimensional structures is critical toward an in-depth understanding of their function and chemical mechanism. Equally important is collecting information on the folding process and the inherent dynamics of a ribozyme fold. Over the past years, Selective-2′-Hydroxyl Acylation analyzed by Primer Extension (SHAPE) turned out to be a significant tool to probe secondary and tertiary interactions of diverse RNA species at the single nucleotide level under varying environmental conditions. Small self-cleaving ribozymes, however, have not been investigated by this method so far. Here, we describe SHAPE probing of pre-catalytic folds of the recently discovered ribozyme classes twister, twister-sister (TS), pistol and hatchet. The study has implications on Mg2+-dependent folding and reveals potentially dynamic residues of these ribozymes that are otherwise difficult to identify. For twister, TS and pistol ribozymes the new findings are discussed in the light of their crystal structures, and in case of twister also with respect to a smFRET folding analysis. For the hatchet ribozyme where an atomic resolution structure is not yet available, the SHAPE data challenge the proposed secondary structure model and point at selected residues and putative long-distance interactions that appear crucial for structure formation and cleavage activity.

StarSeeker: an automated tool for mature duplex microRNA sequence identification based on secondary structure modeling of precursor molecule

BackgroundMicroRNAs (miRNAs) are small, non-coding RNA molecules that play a key role in gene regulation in both plants and animals. MicroRNA biogenesis involves the enzymatic processing of a primary RNA transcript. The final step is the production of a duplex molecule, often designated as miRNA:miRNA*, that will yield a functional miRNA by separation of the two strands. This miRNA will be incorporated into the RNA-induced silencing complex, which subsequently will bind to its target mRNA in order to suppress its expression. The analysis of miRNAs is still a developing area for computational biology with many open questions regarding the structure and function of this important class of molecules. Here, we present StarSeeker, a simple tool that outputs the putative miRNA* sequence given the precursor and the mature sequences.ResultsWe evaluated StarSeeker using a dataset consisting of all plant sequences available in miRBase (6992 precursor sequences and 8496 mature sequences). The program returned a total of 15,468 predicted miRNA* sequences. Of these, 2650 sequences were matched to annotated miRNAs (~ 90% of the miRBase-annotated sequences). The remaining predictions could not be verified, mainly because they do not comply with the rule requiring the two overhanging nucleotides in the duplex molecule.ConclusionsThe expression pattern of some miRNAs in plants can be altered under various abiotic stress conditions. Potential miRNA* molecules that do not degrade can thus be detected and also discovered in high-throughput sequencing data, helping us to understand their role in gene regulation.

Chloromonas svalbardensis n. sp. with Insights into the Phylogroup Chloromonadinia (Chlorophyceae).

The traditional green algal genus Chloromonas accommodates mesophilic, cold-tolerant and cold-adapted microorganisms. In this paper, we studied a new strain isolated from a wet hummock meadow in the High Arctic. We used morphological, ultrastructural and molecular data to assess the taxonomic position and phylogenetic relationships of the new isolate. The observed morphological features generally corresponded to the cold-tolerant Chloromonas characteristics. However, ellipsoidal or wide ellipsoidal vegetative cells, a massive parietal cup-shaped chloroplast with a number of continuously connected lobes, a thick cell wall, a prominent hemispherical papilla and the anterior position of an oblong or round eyespot distinguished the alga from all previously described Chloromonas species. Analyses of rbcL and 18S rRNA genes showed that the new strain formed an independent lineage within a clade containing mesophilic and psychrotolerant Chloromonas species. Comparisons of secondary structure models of a highly variable ITS2 rDNA marker supported a separate species identity of the new isolate. Considering the morphological and molecular differences from its relatives, a new psychrotolerant species, Chloromonas svalbardensis, is proposed. Further, our results demonstrated the paraphyletic origin of Chloromonas within Chloromonadinia with genetically, morphologically and ecologically well-defined clades. We discuss a scenario of a possible Chloromonas split and revision.

In silico characterization of broad range proteases produced by Serratia marcescens EGD-HP20.

In the present study, Serratia marcescens EGD-HP20 strain was demonstrated to utilize poultry waste comprising of both white non-melanized and dark/brown melanized poultry feathers. The potential of the isolate to hydrolyze diverse keratinous wastes was further corroborated by comparative genomics which indicated the presence of genes for broad substrate specific proteases including metallo-proteases, serine endoprotease, dipeptidase, oligopeptidase, etc. Multiple gene sequence alignments of above genes showed 99-100% sequence identities with that of closely related strains of S. marcescens. The secondary structure, 3D structures and energy models suggested the stable nature of all the observed enzymes. Comparative genomics and hydrolysis of mixed feather waste indicated that the above potential of the isolate was associated with synergistic action of various types of proteases.

Chloromonas arctica sp. nov., a psychrotolerant alga from snow in the High Arctic (Chlamydomonadales, Chlorophyta).

With the advent of molecular phylogenetic methods, it has become possible to assess the bioversity of snow algae more accurately. In this study, we focused on a morphological, ultrastructural and taxonomic description of a new Chloromonas-like alga isolated from snow in the High Arctic (Svalbard). Light and transmission electron microscopy revealed broad ellipsoidal or ellipsoidal-cylindrical, occasionally spherical cells with a chloroplast without a pyrenoid, an inconspicuous eyespot and a papilla. The size difference and the aforementioned morphological traits clearly distinguished the alga from its closest counterparts within the genus Chloromonas. Moreover, we were able to cultivate the alga at both 5 and 20 °C, revealing the psychrotolerant nature of the strain. Phylogenetic analyses of the plastid rbcL and nuclear 18S rRNA gene showed that the alga is nested within a clade containing a number of psychrotolerant strains within the Chloromonadinia phylogroup (Chlorophyceae). In the rbcL phylogeny, the alga formed an independent lineage, sister to the freshwater species Chloromonas paraserbinowii. Comparisons of secondary structure models of a highly variable ITS2 rDNA marker showed support for a distinct species identity for the new strain. The ITS2 secondary structure of the new isolate differed from the closest matches 'Chlamydomonas' gerloffii and Choloromonas reticulata by three and five compensatory base changes, respectively. Considering the morphological and molecular differences from its closest relatives, a new psychrotolerant species from the Arctic, Choromonas arctica sp. nov., is proposed.

Structural analyses of NEAT1 lncRNAs suggest long-range RNA interactions that may contribute to paraspeckle architecture.

Paraspeckles are nuclear bodies that regulate multiple aspects of gene expression. The long non-coding RNA (lncRNA) NEAT1 is essential for paraspeckle formation. NEAT1 has a highly ordered spatial organization within the paraspeckle, such that its 5′ and 3′ ends localize on the periphery of paraspeckle, while central sequences of NEAT1 are found within the paraspeckle core. As such, the structure of NEAT1 RNA may be important as a scaffold for the paraspeckle. In this study, we used SHAPE probing and computational analyses to investigate the secondary structure of human and mouse NEAT1. We propose a secondary structural model of the shorter (3,735 nt) isoform hNEAT1_S, in which the RNA folds into four separate domains. The secondary structures of mouse and human NEAT1 are largely different, with the exception of several short regions that have high structural similarity. Long-range base-pairing interactions between the 5′ and 3′ ends of the long isoform NEAT1 (NEAT1_L) were predicted computationally and verified using an in vitro RNA–RNA interaction assay. These results suggest that the conserved role of NEAT1 as a paraspeckle scaffold does not require extensively conserved RNA secondary structure and that long-range interactions among NEAT1 transcripts may have an important architectural function in paraspeckle formation.

Detailed secondary structure models of invertebrate 7SK RNAs

ABSTRACTThe 7SK RNA is a small nuclear RNA that is involved in the regulation of Pol-II transcription. It is very well conserved in vertebrates, but shows extensive variations in both sequence and structure across invertebrates. A systematic homology search extended the collection of 7SK genes in both Arthropods and Lophotrochozoa making use of the large number of recently published invertebrate genomes. The extended data set made it possible to infer complete consensus structures for invertebrate 7SK RNAs. These show that not only the well-conserved 5′- and 3′- domains but all the interior Stem A domain is universally conserved. In contrast, Stem B region exhibits substantial structural variation and does not adhere to a common structural model beyond phylum level.

Influenza virus segment 5 (+)RNA - secondary structure and new targets for antiviral strategies

Influenza A virus is a threat for humans due to seasonal epidemics and occasional pandemics. This virus can generate new strains that are dangerous through nucleotide/amino acid changes or through segmental recombination of the viral RNA genome. It is important to gain wider knowledge about influenza virus RNA to create new strategies for drugs that will inhibit its spread. Here, we present the experimentally determined secondary structure of the influenza segment 5 (+)RNA. Two RNAs were studied: the full-length segment 5 (+)RNA and a shorter construct containing only the coding region. Chemical mapping data combined with thermodynamic energy minimization were used in secondary structure prediction. Sequence/structure analysis showed that the determined secondary structure of segment 5 (+)RNA is mostly conserved between influenza virus type A strains. Microarray mapping and RNase H cleavage identified accessible sites for oligonucleotides in the revealed secondary structure of segment 5 (+)RNA. Antisense oligonucleotides were designed based on the secondary structure model and tested against influenza virus in cell culture. Inhibition of influenza virus proliferation was noticed, identifying good targets for antisense strategies. Effective target sites fall within two domains, which are conserved in sequence/structure indicating their importance to the virus.

Metal transporter encoding gene families in Fabaceae: II. Cation/H+ exchanger (CAX) encoding genes

The plant CAtion/H+ eXchangers (CAX) proteins belong to Ca2+/cation antiporter (CaCA) superfamily. By using in silico methods, the CAX encoding genes in the genome of six legume species have been identified in this work. In examined legume genomes, the CAX genes belong to a small multigenic family. The number of the CAX genes in these legume species is 17 (soybean), 6 (common bean and C. cajan), 5 (M. truncatula and C. arietinum) and 3 genes (L. japonicus), respectively. The legume CAX genes vary in genomic full-length ranging from 1,213 to 11,561 base pairs. All of the genes exhibit introns (from 4 to 11 introns). Their deduced full-length protein sequences range from 248 to 718 amino acids. Theoretical pI values of most (39/42) of legume CAX proteins were less than 7. The secondary structure modelling of protein exhibit transmembrane helix region (from 3 to 11 regions). Half of all (23/42) included 11 transmembrane helix regions. Based on phylogeny analysis, all of the legume CAX were divided into two groups, A and B, each consisting of two subgroups. The phylogeny suggested an ancient gene duplication in the genome of legumes ancestry. The recent gene duplication even was only detected in the soybean genome after the speciation. The expression analysis showed that all of 3 L. japonicus CAX genes expressed in all examined tissues. However, the expression of C. cajan CAX genes was not detected. For each of 4 remaining legumes, the CAX genes were differed in their expression level depending on studied tissues. The tissue-specific expressions of some CAX genes were observed in 5 out of the 6 legume species, except C. cajan.

Design of Artificial Riboswitches as Biosensors.

RNA aptamers readily recognize small organic molecules, polypeptides, as well as other nucleic acids in a highly specific manner. Many such aptamers have evolved as parts of regulatory systems in nature. Experimental selection techniques such as SELEX have been very successful in finding artificial aptamers for a wide variety of natural and synthetic ligands. Changes in structure and/or stability of aptamers upon ligand binding can propagate through larger RNA constructs and cause specific structural changes at distal positions. In turn, these may affect transcription, translation, splicing, or binding events. The RNA secondary structure model realistically describes both thermodynamic and kinetic aspects of RNA structure formation and refolding at a single, consistent level of modelling. Thus, this framework allows studying the function of natural riboswitches in silico. Moreover, it enables rationally designing artificial switches, combining essentially arbitrary sensors with a broad choice of read-out systems. Eventually, this approach sets the stage for constructing versatile biosensors.

The IRES5'UTR of the dicistrovirus cricket paralysis virus is a type III IRES containing an essential pseudoknot structure.

Cricket paralysis virus (CrPV) is a dicistrovirus. Its positive-sense single-stranded RNA genome contains two internal ribosomal entry sites (IRESs). The 5′ untranslated region (5′UTR) IRES5′UTR mediates translation of non-structural proteins encoded by ORF1 whereas the well-known intergenic region (IGR) IRESIGR is required for translation of structural proteins from open reading frame 2 in the late phase of infection. Concerted action of both IRES is essential for host translation shut-off and viral translation. IRESIGR has been extensively studied, in contrast the IRES5′UTR remains largely unexplored. Here, we define the minimal IRES element required for efficient translation initiation in drosophila S2 cell-free extracts. We show that IRES5′UTR promotes direct recruitment of the ribosome on the cognate viral AUG start codon without any scanning step, using a Hepatitis-C virus-related translation initiation mechanism. Mass spectrometry analysis revealed that IRES5′UTR recruits eukaryotic initiation factor 3, confirming that it belongs to type III class of IRES elements. Using Selective 2′-hydroxyl acylation analyzed by primer extension and DMS probing, we established a secondary structure model of 5′UTR and of the minimal IRES5′UTR. The IRES5′UTR contains a pseudoknot structure that is essential for proper folding and ribosome recruitment. Overall, our results pave the way for studies addressing the synergy and interplay between the two IRES from CrPV.

Visualization of conformational variability in the domains of long single-stranded RNA molecules.

We demonstrate an application of atomic force microscopy (AFM) for the structural analysis of long single-stranded RNA (>1 kb), focusing on 28S ribosomal RNA (rRNA). Generally, optimization of the conditions required to obtain three-dimensional (3D) structures of long RNA molecules is a challenging or nearly impossible process. In this study, we overcome these limitations by developing a method using AFM imaging combined with automated, MATLAB-based image analysis algorithms for extracting information about the domain organization of single RNA molecules. We examined the 5 kb human 28S rRNA since it is the largest RNA molecule for which a 3D structure is available. As a proof of concept, we determined a domain structure that is in accordance with previously described secondary structural models. Importantly, we identified four additional small (200–300 nt), previously unreported domains present in these molecules. Moreover, the single-molecule nature of our method enabled us to report on the relative conformational variability of each domain structure identified, and inter-domain associations within subsets of molecules leading to molecular compaction, which may shed light on the process of how these molecules fold into the final tertiary structure.

A systematic evaluation of nucleotide properties for CRISPR sgRNA design

BackgroundCRISPR is a versatile gene editing tool which has revolutionized genetic research in the past few years. Optimizing sgRNA design to improve the efficiency of target/DNA cleavage is critical to ensure the success of CRISPR screens.ResultsBy borrowing knowledge from oligonucleotide design and nucleosome occupancy models, we systematically evaluated candidate features computed from a number of nucleic acid, thermodynamic and secondary structure models on real CRISPR datasets. Our results showed that taking into account position-dependent dinucleotide features improved the design of effective sgRNAs with area under the receiver operating characteristic curve (AUC) >0.8, and the inclusion of additional features offered marginal improvement (∼2% increase in AUC).ConclusionUsing a machine-learning approach, we proposed an accurate prediction model for sgRNA design efficiency. An R package predictSGRNA implementing the predictive model is available at http://www.ams.sunysb.edu/~pfkuan/softwares.html#predictsgrna.