ABSTRACT
Sequential lytic cycles driven by cascading transcriptional waves underlie pathogenesis in the apicomplexan parasite Toxoplasma gondii. This parasite's unique division by internal budding, short cell cycle, and jumbled up classically defined cell cycle stages have restrained in-depth transcriptional program analysis. Here, unbiased transcriptome and chromatin accessibility maps throughout the lytic cell cycle are established at the single-cell level. Correlated pseudo-timeline assemblies of expression and chromatin profiles maps transcriptional versus chromatin level transition points promoting the cell division cycle. Sequential clustering analysis identifies functionally related gene groups promoting cell cycle progression. Promoter DNA motif mapping reveals patterns of combinatorial regulation. Pseudo-time trajectory analysis reveals transcriptional bursts at different cell cycle points. The dominant burst in G1 is driven largely by transcription factor AP2XII-8, which engages a conserved DNA motif, and promotes the expression of 44 ribosomal proteins encoding regulon. Overall, the study provides integrated, multi-level insights into apicomplexan transcriptional regulation.
Subject(s)
Chromatin , Protozoan Proteins , Regulon , Ribosomes , Single-Cell Analysis , Toxoplasma , Toxoplasma/genetics , Toxoplasma/metabolism , Chromatin/metabolism , Chromatin/genetics , Regulon/genetics , Protozoan Proteins/metabolism , Protozoan Proteins/genetics , Ribosomes/metabolism , Ribosomes/genetics , Gene Expression Regulation , Promoter Regions, Genetic/genetics , Cell Cycle/genetics , Humans , Nucleotide Motifs/genetics , Transcriptome , Ribosomal Proteins/metabolism , Ribosomal Proteins/geneticsABSTRACT
Fluorescence in situ hybridization (FISH) with two different probes, the canonical insect telomeric sequence (TTAGG)n and the sequence (TTAGGGATGG)n, was performed on meiotic chromosomes of two members of the true bug family Cimicidae (Cimicomorpha), the common bed bug Cimex lectularius Linnaeus, 1758 and the tropical bed bug C. hemipterus (Fabricius, 1803), whose telomeric motifs were not known. In both species, there were no hybridization signals with the first probe, but strong signals at chromosomal ends were observed with the second probe, indicating the presence of a telomeric motif (TTAGGGATGG)n. This study represents the first FISH confirmation of the presence of a non-canonical telomeric motif not only for the infraorder Cimicomorpha but also for the suborder Heteroptera (Hemiptera) as a whole. The present finding is of key significance for unraveling the evolutionary shifts in the telomeric sequences in this suborder.
Subject(s)
In Situ Hybridization, Fluorescence , Telomere , Animals , Telomere/genetics , Heteroptera/genetics , Nucleotide Motifs/genetics , Chromosomes, Insect/geneticsABSTRACT
Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.
Subject(s)
Gene Expression Regulation , Nucleotide Motifs , Promoter Regions, Genetic , Transcription Factors , Transcription Initiation Site , Transcription Initiation, Genetic , Humans , Binding Sites , Gene Expression Regulation/genetics , Genome, Human/genetics , Nucleotide Motifs/genetics , Promoter Regions, Genetic/genetics , Protein Binding , Transcription Factors/metabolism , Genetic VariationABSTRACT
To combat cancer, a comprehensive understanding of the molecular mechanisms and behaviors involved in carcinogenesis is crucial, as tumorigenesis is a complex process influenced by various genetic events and disease hallmarks. The B-MYB gene encodes a transcription factor involved in cell cycle regulation, survival, and differentiation in normal cells. B-MYB can be transformed into an oncogene through mutations, and abnormal expression of B-MYB has been identified in various cancers, including lung cancer, and is associated with poor prognosis. Targeting this oncogene is a promising approach for anti-cancer drug design. B-MYB has been deemed undruggable in previous reports, necessitating the search for novel therapeutic options. In this study, we found that the B-MYB gene promoter contains several G/C rich motifs compatible with G-quadruplex (G4) formation. We investigated and validated the existence of G4 structures in the promoter region of B-MYB, first in vitro using a combination of bioinformatics, biophysical, and biochemical methods, then in cell with the recently developed G4access method.
Subject(s)
G-Quadruplexes , Promoter Regions, Genetic , Proto-Oncogene Mas , Promoter Regions, Genetic/genetics , Humans , Trans-Activators/genetics , Trans-Activators/metabolism , Cell Cycle Proteins/genetics , Cell Cycle Proteins/metabolism , Nucleotide Motifs/geneticsABSTRACT
BACKGROUND: Trans-acting factors are of special importance in transcription regulation, which is a group of proteins that can directly or indirectly recognize or bind to the 8-12 bp core sequence of cis-acting elements and regulate the transcription efficiency of target genes. The progressive development in high-throughput chromatin capture technology (e.g., Hi-C) enables the identification of chromatin-interacting sequence groups where trans-acting DNA motif groups can be discovered. The problem difficulty lies in the combinatorial nature of DNA sequence pattern matching and its underlying sequence pattern search space. METHOD: Here, we propose to develop MotifHub for trans-acting DNA motif group discovery on grouped sequences. Specifically, the main approach is to develop probabilistic modeling for accommodating the stochastic nature of DNA motif patterns. RESULTS: Based on the modeling, we develop global sampling techniques based on EM and Gibbs sampling to address the global optimization challenge for model fitting with latent variables. The results reflect that our proposed approaches demonstrate promising performance with linear time complexities. CONCLUSION: MotifHub is a novel algorithm considering the identification of both DNA co-binding motif groups and trans-acting TFs. Our study paves the way for identifying hub TFs of stem cell development (OCT4 and SOX2) and determining potential therapeutic targets of prostate cancer (FOXA1 and MYC). To ensure scientific reproducibility and long-term impact, its matrix-algebra-optimized source code is released at http://bioinfo.cs.cityu.edu.hk/MotifHub.
Subject(s)
Algorithms , Software , Nucleotide Motifs/genetics , Reproducibility of Results , Chromatin/geneticsABSTRACT
How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.
Subject(s)
Mutation Rate , Neoplasms , Humans , Nucleotide Motifs/genetics , Mutation/genetics , Epigenesis, Genetic/genetics , Neoplasms/geneticsABSTRACT
RNA is a pivotal molecule that plays critical roles in various cellular processes. Quantifying RNA structures and interactions is essential to understanding RNA function and developing RNA-based therapeutics. Using a unified five-bead model and a non-redundant database, this paper investigates the structural features and interactions of five commonly occurring RNA motifs, i.e., double-stranded helices, hairpin loops, internal/bulge loops, multi-branched junctions, and single-stranded terminal tails. Analyzing detailed distributions of RNA local structural features and base-base interactions reveals a preference for helical structures in both local backbone structures and base orientations. The interactions between adjacent bases exhibit motif-specific and sequence-dependent characteristics, reflecting the distinct topological constraints imposed by different loop-helix connection modes and the varying pairing and stacking interactions among different sequences. These findings shed light on the stability of RNA helices, emphasizing their significance in providing dominant base pairing and stacking interactions for RNA structures and stability. The four non-helix motifs encompass unpaired nucleotide loops and exhibit diverse base-base interactions, contributing to the structural diversity observed in RNA. Overall, the complexity of RNA structure arises from the intricate interplay of base-base interactions.
Subject(s)
RNA , RNA/genetics , RNA/chemistry , Nucleic Acid Conformation , Models, Molecular , Base Pairing , Nucleotide Motifs/geneticsABSTRACT
RNAs that are able to prevent degradation by the 5'-3' exoribonuclease Xrn1 have emerged as crucial structures during infection by an increasing number of RNA viruses. Several plant viruses employ the so-called coremin motif, an Xrn1-resistant RNA that is usually located in 3' untranslated regions. Investigation of its structural and sequence requirements has led to its identification in plant virus families beyond those in which the coremin motif was initially discovered. In this study, we identified coremin-like motifs that deviate from the original in the number of nucleotides present in the loop region of the 5' proximal hairpin. They are present in a number of viral families that previously did not have an Xrn1-resistant RNA identified yet, including the double-stranded RNA virus families Hypoviridae and Chrysoviridae. Through systematic mutational analysis, we demonstrated that a coremin motif carrying a 6-nucleotide loop in the 5' proximal hairpin generally requires a YGNNAD consensus for stalling Xrn1, similar to the previously determined YGAD consensus required for Xrn1 resistance of the original coremin motif. Furthermore, we determined the minimal requirements for the 3' proximal hairpin. Since some putative coremin motifs were found in intergenic regions or coding sequences, we demonstrated their capacity for inhibiting translation through an in vitro ribosomal scanning inhibition assay. Consequently, this study provides a further expansion on the number of viral families with known Xrn1-resistant elements, while adding a novel, potentially regulatory function for this structure.
Subject(s)
Plant Viruses , RNA, Viral , Nucleotide Motifs/genetics , RNA, Viral/metabolism , Exoribonucleases/metabolism , Virome , Ribosomes/metabolism , Nucleotides , Plant Viruses/genetics , Plant Viruses/metabolism , Nucleic Acid Conformation , RNA StabilityABSTRACT
Advances in sequencing technologies have empowered epitranscriptomic profiling at the single-base resolution. Putative RNA modification sites identified from a single high-throughput experiment may contain one type of modification deposited by different writers or different types of modifications, along with false positive results because of the challenge of distinguishing signals from noise. However, current tools are insufficient for subtyping, visualization, and denoising these signals. Here, we present iMVP, which is an interactive framework for epitranscriptomic analysis with a nonlinear dimension reduction technique and density-based partition. As exemplified by the analysis of mRNA m5C and ModTect variant data, we show that iMVP allows the identification of previously unknown RNA modification motifs and writers and the discovery of false positives that are undetectable by traditional methods. Using putative m6A/m6Am sites called from 8 profiling approaches, we illustrate that iMVP enables comprehensive comparison of different approaches and advances our understanding of the difference and pattern of true positives and artifacts in these methods. Finally, we demonstrate the ability of iMVP to analyze an extremely large human A-to-I editing dataset that was previously unmanageable. Our work provides a general framework for the visualization and interpretation of epitranscriptomic data.
Subject(s)
Artifacts , Technology , Humans , Nucleotide Motifs/genetics , RNA, MessengerABSTRACT
The stem-loop II motif (s2m) is an RNA structural element that is found in the 3' untranslated region (UTR) of many RNA viruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Though the motif was discovered over 25 years ago, its functional significance is unknown. In order to understand the importance of s2m, we created viruses with deletions or mutations of the s2m by reverse genetics and also evaluated a clinical isolate harboring a unique s2m deletion. Deletion or mutation of the s2m had no effect on growth in vitro or on growth and viral fitness in Syrian hamsters in vivo. We also compared the secondary structure of the 3' UTR of wild-type and s2m deletion viruses using selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) and dimethyl sulfate mutational profiling and sequencing (DMS-MaPseq). These experiments demonstrate that the s2m forms an independent structure and that its deletion does not alter the overall remaining 3'-UTR RNA structure. Together, these findings suggest that s2m is dispensable for SARS-CoV-2. IMPORTANCE RNA viruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), contain functional structures to support virus replication, translation, and evasion of the host antiviral immune response. The 3' untranslated region of early isolates of SARS-CoV-2 contained a stem-loop II motif (s2m), which is an RNA structural element that is found in many RNA viruses. This motif was discovered over 25 years ago, but its functional significance is unknown. We created SARS-CoV-2 with deletions or mutations of the s2m and determined the effect of these changes on viral growth in tissue culture and in rodent models of infection. Deletion or mutation of the s2m element had no effect on growth in vitro or on growth and viral fitness in Syrian hamsters in vivo. We also observed no impact of the deletion on other known RNA structures in the same region of the genome. These experiments demonstrate that s2m is dispensable for SARS-CoV-2.
Subject(s)
Nucleotide Motifs , SARS-CoV-2 , Animals , Cricetinae , 3' Untranslated Regions/genetics , COVID-19/virology , Mesocricetus , Mutation , SARS-CoV-2/genetics , Nucleotide Motifs/genetics , RNA, Viral/chemistry , RNA, Viral/geneticsABSTRACT
Discovering DNA regulatory sequence motifs and their relative positions is vital to understanding the mechanisms of gene expression regulation. Although deep convolutional neural networks (CNNs) have achieved great success in predicting cis-regulatory elements, the discovery of motifs and their combinatorial patterns from these CNN models has remained difficult. We show that the main difficulty is due to the problem of multifaceted neurons which respond to multiple types of sequence patterns. Since existing interpretation methods were mainly designed to visualize the class of sequences that can activate the neuron, the resulting visualization will correspond to a mixture of patterns. Such a mixture is usually difficult to interpret without resolving the mixed patterns. We propose the NeuronMotif algorithm to interpret such neurons. Given any convolutional neuron (CN) in the network, NeuronMotif first generates a large sample of sequences capable of activating the CN, which typically consists of a mixture of patterns. Then, the sequences are "demixed" in a layer-wise manner by backward clustering of the feature maps of the involved convolutional layers. NeuronMotif can output the sequence motifs, and the syntax rules governing their combinations are depicted by position weight matrices organized in tree structures. Compared to existing methods, the motifs found by NeuronMotif have more matches to known motifs in the JASPAR database. The higher-order patterns uncovered for deep CNs are supported by the literature and ATAC-seq footprinting. Overall, NeuronMotif enables the deciphering of cis-regulatory codes from deep CNs and enhances the utility of CNN in genome interpretation.
Subject(s)
Algorithms , Neural Networks, Computer , Nucleotide Motifs/genetics , Regulatory Sequences, Nucleic Acid/genetics , Databases, FactualABSTRACT
Non-B DNA structures formed by repetitive sequence motifs are known instigators of mutagenesis in experimental systems. Analyzing this phenomenon computationally in the human genome requires careful disentangling of intrinsic confounding factors, including overlapping and interrupted motifs and recurrent sequencing errors. Here, we show that accounting for these factors eliminates all signals of repeat-induced mutagenesis that extend beyond the motif boundary, and eliminates or dramatically shrinks the magnitude of mutagenesis within some motifs, contradicting previous reports. Mutagenesis not attributable to artifacts revealed several biological mechanisms. Polymerase slippage generates frequent indels within every variety of short tandem repeat motif, implicating slipped-strand structures. Interruption-correcting single nucleotide variants within short tandem repeats may originate from error-prone polymerases. Secondary-structure formation promotes single nucleotide variants within palindromic repeats and duplications within direct repeats. G-quadruplex motifs cause recurrent sequencing errors, whereas mutagenesis at Z-DNAs is conspicuously absent.
Subject(s)
DNA , Genome, Human , Humans , Nucleotide Motifs/genetics , Mutagenesis , DNA/genetics , DNA/chemistry , NucleotidesABSTRACT
RNA recognition motifs (RRM) are the most prevalent class of RNA binding domains in eucaryotes. Their RNA binding preferences have been investigated for almost two decades, and even though some RRM domains are now very well described, their RNA recognition code has remained elusive. An increasing number of experimental structures of RRM-RNA complexes has become available in recent years. Here, we perform an in-depth computational analysis to derive an RNA recognition code for canonical RRMs. We present and validate a computational scoring method to estimate the binding between an RRM and a single stranded RNA, based on structural data from a carefully curated multiple sequence alignment, which can predict RRM binding RNA sequence motifs based on the RRM protein sequence. Given the importance and prevalence of RRMs in humans and other species, this tool could help design RNA binding motifs with uses in medical or synthetic biology applications, leading towards the de novo design of RRMs with specific RNA recognition.
Subject(s)
RNA Recognition Motif , RNA , Humans , RNA/chemistry , Amino Acid Sequence , Sequence Alignment , Nucleotide Motifs/genetics , Protein Binding , Binding SitesABSTRACT
A DNA motif is a sequence pattern shared by the DNA sequence segments that bind to a specific protein. Discovering motifs in a given DNA sequence dataset plays a vital role in studying gene expression regulation. As an important attribute of the DNA motif, the motif length directly affects the quality of the discovered motifs. How to determine the motif length more accurately remains a difficult challenge to be solved. We propose a new motif length prediction scheme named MotifLen by using supervised machine learning. First, a method of constructing sample data for predicting the motif length is proposed. Secondly, a deep learning model for motif length prediction is constructed based on the convolutional neural network. Then, the methods of applying the proposed prediction model based on a motif found by an existing motif discovery algorithm are given. The experimental results show that i) the prediction accuracy of MotifLen is more than 90% on the validation set and is significantly higher than that of the compared methods on real datasets, ii) MotifLen can successfully optimize the motifs found by the existing motif discovery algorithms, and iii) it can effectively improve the time performance of some existing motif discovery algorithms.
Subject(s)
Deep Learning , Nucleotide Motifs/genetics , Sequence Analysis, DNA/methods , Algorithms , Neural Networks, ComputerABSTRACT
RNA motif classification is important for understanding structure/function connections and building phylogenetic relationships. Using our coarse-grained RNA-As-Graphs (RAG) representations, we identify recurrent dual graph motifs in experimentally solved RNA structures based on an improved search algorithm that finds and ranks independent RNA substructures. Our expanded list of 183 existing dual graph motifs reveals five common motifs found in transfer RNA, riboswitch, and ribosomal 5S RNA components. Moreover, we identify three motifs for available viral frameshifting RNA elements, suggesting a correlation between viral structural complexity and frameshifting efficiency. We further partition the RNA substructures into 1844 distinct submotifs, with pseudoknots and junctions retained intact. Common modules are internal loops and three-way junctions, and three submotifs are associated with riboswitches that bind nucleotides, ions, and signaling molecules. Together, our library of existing RNA motifs and submotifs adds to the growing universe of RNA modules, and provides a resource of structures and substructures for novel RNA design.
Subject(s)
RNA , Riboswitch , Algorithms , Gene Library , Nucleic Acid Conformation , Nucleotide Motifs/genetics , Phylogeny , RNA/chemistry , RNA/genetics , RNA, Viral/geneticsABSTRACT
Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F1-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.
There are hundreds of different types of cells in the body. Each one performs a unique role, but they all share the same genes. Sequences of the genetic code called enhancers decide which genes each cell uses. Enhancers work like genetic switches: to turn a gene on, proteins called transcription factors assemble on an enhancer. Each transcription factor recognises a short sequence on the enhancer, and several distinct transcription factors work together to promote the activatation of a gene. The relationship between transcription factors, enhancers, and gene activation is complex. The specific genetic sequences of enhancers differ between species, changing the way these genetic switches work. But scientists are not yet able to reliably predict the effects of small changes in the DNA sequence of an enhancer. One way to tackle this problem is to look at different versions of the same enhancers side by side to see how small mutations change their behaviour. Mammalian cells generally carry two copies of each chromosome (the molecules that contain the genetic code), one inherited from each parent. Each of the two copies carries the same genes and enhancers, but there are many small differences in the DNA sequences of enhancers between the chromosomes inherited from each parent, which can potentially alter their function Yang, Ling et al. generated cells from mice that come from different inbred strains, which are similar to purebred dogs. By breeding two distinct inbred mouse strains together that are very different from one another, they generated a panel of hybrid mouse cell lines that have a relatively large number of differences in their DNA sequence between the maternal and paternal chromosomes. Looking at the different versions of each enhancer side-by-side revealed thousands of single letter changes in the DNA sequence of enhancers that changed how they work. Mutations affecting the binding site of one transcription factor within an enhancer can indirectly affect the binding of other types of transcription factors. Yang, Ling et al. found that if a transcription factor could no longer find its place on an enhancer, it stopped others from binding even if their own places had not changed. Sometimes, mutations on either side of the binding sequences also affected transcription factor binding. This suggests a more complex relationship than previously thought may exist between the DNA sequence of an enhancer and the transcription factors that bind to it. Spotting the differences caused by mutations could help further the efforts of scientists to read and write the genetic code. This could have many benefits. It would allow scientists to control natural or artificial genes, and to predict the effects of genetic changes that are identified in humans with genetic diseases. This might improve genetic experiments, medical screening, gene therapy, and our understanding of evolution.
Subject(s)
Enhancer Elements, Genetic , Genetic Variation , Transcription Factor AP-1 , Animals , Humans , Mice , Binding Sites/genetics , Enhancer Elements, Genetic/genetics , Genetic Variation/genetics , Nucleotide Motifs/genetics , Protein Binding/genetics , Transcription Factor AP-1/geneticsABSTRACT
Non-coding RNAs are key regulatory players in bacteria. Many computationally predicted non-coding RNAs, however, lack functional associations. An example is the Bacillaceae-1 RNA motif, whose Rfam model consists of two hairpin loops. We find the motif conserved in nine of 13 non-pathogenic strains of the genus Bacillus but only in one pathogenic strain. To elucidate functional characteristics, we studied 118 hits of the Rfam model in 11 Bacillus spp. and found two distinct classes based on the ensemble diversity of their RNA secondary structure and the genomic context concerning the ribosomal RNA (rRNA) cluster. Forty hits are associated with the rRNA cluster, of which all 19 hits upstream flanking of 16S rRNA have a reverse complementary structure of low structural diversity. Fifty-two hits have large ensemble diversity, of which 38 are located between two coding genes. For eight hits in Bacillus subtilis, we investigated public expression data under various conditions and observed either the forward or the reverse complementary motif expressed. Five hits are associated with the rRNA cluster. Four of them are located upstream of the 16S rRNA and are not transcriptionally active, but instead, their reverse complements with low structural diversity are expressed together with the rRNA cluster. The three other hits are located between two coding genes in non-conserved genomic loci. Two of them are independently expressed from their surrounding genes and are structurally diverse. In summary, we found that Bacillaceae-1 RNA motifs upstream flanking of ribosomal RNA clusters tend to have one stable structure with the reverse complementary motif expressed in B. subtilis. In contrast, a subgroup of intergenic motifs has the thermodynamic potential for structural switches.
Subject(s)
Bacillaceae , Bacillus , Bacillaceae/genetics , Bacillaceae/metabolism , Bacillus/genetics , Bacillus subtilis/genetics , Nucleotide Motifs/genetics , Phylogeny , RNA, Ribosomal/genetics , RNA, Ribosomal, 16S/geneticsABSTRACT
Influenza A virus (IAV) is a member of the single-stranded RNA (ssRNA) family of viruses. The most recent global pandemic caused by the SARS-CoV-2 virus has shown the major threat that RNA viruses can pose to humanity. In comparison, influenza has an even higher pandemic potential as a result of its high rate of mutations within its relatively short (<13 kbp) genome, as well as its capability to undergo genetic reassortment. In light of this threat, and the fact that RNA structure is connected to a broad range of known biological functions, deeper investigation of viral RNA (vRNA) structures is of high interest. Here, for the first time, we propose a secondary structure for segment 8 vRNA (vRNA8) of A/California/04/2009 (H1N1) formed in the presence of cellular and viral components. This structure shows similarities with prior in vitro experiments. Additionally, we determined the location of several well-defined, conserved structural motifs of vRNA8 within IAV strains with possible functionality. These RNA motifs appear to fold independently of regional nucleoprotein (NP)-binding affinity, but a low or uneven distribution of NP in each motif region is noted. This research also highlights several accessible sites for oligonucleotide tools and small molecules in vRNA8 in a cellular environment that might be a target for influenza A virus inhibition on the RNA level.
Subject(s)
Gene Expression Regulation, Viral , Genome, Viral/genetics , Influenza A Virus, H1N1 Subtype/genetics , Nucleic Acid Conformation , RNA, Viral/chemistry , Animals , Base Sequence , Dogs , Humans , Influenza A Virus, H1N1 Subtype/metabolism , Influenza, Human/virology , Madin Darby Canine Kidney Cells , Models, Molecular , Nucleotide Motifs/genetics , RNA Folding , RNA, Viral/genetics , Viral Proteins/genetics , Viral Proteins/metabolismABSTRACT
The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.
Subject(s)
Genome, Plant , Tracheophyta/genetics , Transcription Initiation Site , Base Composition/genetics , Binding Sites , DNA, Plant/genetics , Exons/genetics , Molecular Sequence Annotation , Nucleotide Motifs/genetics , Nucleotides/metabolism , Open Reading Frames/genetics , Promoter Regions, Genetic , Transcription Factors/metabolismABSTRACT
Highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) with enhanced replication capability emerged in China and has become dominant epidemic strain since 2006. Up to now, the replication-regulated genes of PRRSV have not been fully clarified. Here, by swapping the genes or elements between HP-PRRSV and classical PRRSV based on infectious clones, NSP1, NSP2, NSP7, NSP9 and 3'-UTR are found to contribute to the high replication efficiency of HP-PRRSV. Further study revealed that mutations at positions 117th or 119th in the 3'-UTR are significantly related to replication efficiency, and the nucleotide at position 120th is critical for viral rescue. The motif composed by 117-120th nucleotides was quite conservative within each lineage of PRRSV; mutations in the motif of HP-PRRSV and currently epidemic lineage 1 (L1) PRRSV showed higher synthesis ability of viral negative genomic RNA, suggesting that those mutations were beneficial for viral replication. RNA structure analysis revealed that this motif maybe involved into a pseudoknot in the 3'-UTR. The results discovered a novel motif, 117-120th nucleotide in the 3'-UTR, that is critical for replication of PRRSV-2, and mutations in the motif contribute to the enhanced replicative ability of HP-PRRSV or L1 PRRSV. Our findings will help to understand the molecular basis of PRRSV replication and find the potential factors resulting in an epidemic strain of PRRSV.