Unpaired Regions Research Articles

The evolution of noncoding RNAs in the Drosophila melanogaster genome

Noncoding RNAs (ncRNAs) are important regulatory molecules that function owing to their specific secondary structures. The secondary structure depends on the RNA sequence and can be viewed as a combination of loops, stems, and pseudoknot regions. Each local RNA part, as well as the entire molecule, tends to achieve the minimal free energy by assuming a certain structure. New mutations may change the free energy of the current RNA structure so that another structure becomes optimal, thus altering the RNA function. The chance that a mutation will cause structural alterations depends on the mutation position in the RNA secondary structure. Divergence estimates have been found to differ between paired and unpaired regions, the difference arises because selection is aimed at preserving the ncRNA secondary structure. A more detailed analysis of the evolutionary differences between loops and stems in all of the most important ncRNA classes may be of interest for developing the algorithms to search a genome for ncRNA genes. This study focused on natural selection acting on the loop and stem regions in five ncRNA classes. An attempt was made to associate selection characteristics with changes in free energy. The Drosophila melanogaster genome was used as a model.

On the biased nucleotide composition of the human coronavirus RNA genome

We investigated the nucleotide composition of the RNA genome of the six human coronaviruses. Some general coronavirus characteristics were apparent (e.g. high U, low C count), but we also detected species-specific signatures. Most strikingly, the high U and low C proportions are quite variable and act like communicating vessels, C goes down when U goes up and vice versa. U ranges among virus isolates from 30.7% to 40.3%, and C makes the opposite movement from 20.0% to 12.9%, respectively. The nucleotide biases are more pronounced in the unpaired regions of the structured RNA genome, which may suggest a certain biological function for these distinctive sequence signatures. Coronaviruses have an atypical codon usage that has been linked to mutational events operating on the viral RNA genome on an evolutionary time scale. We suggest that the atypical nucleotide bias may serve a distinct biological function and that it is the direct cause of the characteristic codon usage in these viruses. The relevance for evolution of the novel human pathogens MERS and SARS is discussed.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

A Syntenic Region Conserved from Fish to Mammalian X Chromosome

Sex chromosomes bearing the sex-determining gene initiate development along the male or female pathway, no matter which sex is determined by XY male or ZW female heterogamety. Sex chromosomes originate from ancient autosomes but evolved rapidly after the acquisition of sex-determining factors which are highly divergent between species. In the heterogametic male system (XY system), the X chromosome is relatively evolutionary silent and maintains most of its ancestral genes, in contrast to its Y counterpart that has evolved rapidly and degenerated. Sex in a teleost fish, the Nile tilapia (Oreochromis niloticus), is determined genetically via an XY system, in which an unpaired region is present in the largest chromosome pair. We defined the differences in DNA contents present in this chromosome with a two-color comparative genomic hybridization (CGH) and the random amplified polymorphic DNA (RAPD) approach in XY males. We further identified a syntenic segment within this region that is well conserved in several teleosts. Through comparative genome analysis, this syntenic segment was also shown to be present in mammalian X chromosomes, suggesting a common ancestral origin of vertebrate sex chromosomes.

Direct entry by RNase E is a major pathway for the degradation and processing of RNA in Escherichia coli.

Escherichia coli endoribonuclease E has a major influence on gene expression. It is essential for the maturation of ribosomal and transfer RNA as well as the rapid degradation of messenger RNA. The latter ensures that translation closely follows programming at the level of transcription. Recently, one of the hallmarks of RNase E, i.e. its ability to bind via a 5′-monophosphorylated end, was shown to be unnecessary for the initial cleavage of some polycistronic tRNA precursors. Here we show using RNA-seq analyses of ribonuclease-deficient strains in vivo and a 5′-sensor mutant of RNase E in vitro that, contrary to current models, 5′-monophosphate-independent, ‘direct entry’ cleavage is a major pathway for degrading and processing RNA. Moreover, we present further evidence that direct entry is facilitated by RNase E binding simultaneously to multiple unpaired regions. These simple requirements may maximize the rate of degradation and processing by permitting multiple sites to be surveyed directly without being constrained by 5′-end tethering. Cleavage was detected at a multitude of sites previously undescribed for RNase E, including ones that regulate the activity and specificity of ribosomes. A potentially broad role for RNase G, an RNase E paralogue, in the trimming of 5′-monophosphorylated ends was also revealed.

RNA graph partitioning for the discovery of RNA modularity: a novel application of graph partition algorithm to biology.

Graph representations have been widely used to analyze and design various economic, social, military, political, and biological networks. In systems biology, networks of cells and organs are useful for understanding disease and medical treatments and, in structural biology, structures of molecules can be described, including RNA structures. In our RNA-As-Graphs (RAG) framework, we represent RNA structures as tree graphs by translating unpaired regions into vertices and helices into edges. Here we explore the modularity of RNA structures by applying graph partitioning known in graph theory to divide an RNA graph into subgraphs. To our knowledge, this is the first application of graph partitioning to biology, and the results suggest a systematic approach for modular design in general. The graph partitioning algorithms utilize mathematical properties of the Laplacian eigenvector (µ2) corresponding to the second eigenvalues (λ2) associated with the topology matrix defining the graph: λ2 describes the overall topology, and the sum of µ2′s components is zero. The three types of algorithms, termed median, sign, and gap cuts, divide a graph by determining nodes of cut by median, zero, and largest gap of µ2′s components, respectively. We apply these algorithms to 45 graphs corresponding to all solved RNA structures up through 11 vertices (∼220 nucleotides). While we observe that the median cut divides a graph into two similar-sized subgraphs, the sign and gap cuts partition a graph into two topologically-distinct subgraphs. We find that the gap cut produces the best biologically-relevant partitioning for RNA because it divides RNAs at less stable connections while maintaining junctions intact. The iterative gap cuts suggest basic modules and assembly protocols to design large RNA structures. Our graph substructuring thus suggests a systematic approach to explore the modularity of biological networks. In our applications to RNA structures, subgraphs also suggest design strategies for novel RNA motifs.

Efficient detection of unpaired DNA requires a member of the rad54-like family of homologous recombination proteins.

Meiotic silencing by unpaired DNA (MSUD) is a process that detects unpaired regions between homologous chromosomes and silences them for the duration of sexual development. While the phenomenon of MSUD is well recognized, the process that detects unpaired DNA is poorly understood. In this report, we provide two lines of evidence linking unpaired DNA detection to a physical search for DNA homology. First, we have found that a putative SNF2-family protein (SAD-6) is required for efficient MSUD in Neurospora crassa. SAD-6 is closely related to Rad54, a protein known to facilitate key steps in the repair of double-strand breaks by homologous recombination. Second, we have successfully masked unpaired DNA by placing identical transgenes at slightly different locations on homologous chromosomes. This masking falls apart when the distance between the transgenes is increased. We propose a model where unpaired DNA detection during MSUD is achieved through a spatially constrained search for DNA homology. The identity of SAD-6 as a Rad54 paralog suggests that this process may be similar to the searching mechanism used during homologous recombination.

Structure-specific ribonucleases for MS-based elucidation of higher-order RNA structure.

Supported by high-throughput sequencing technologies, structure-specific nucleases are experiencing a renaissance as biochemical probes for genome-wide mapping of nucleic acid structure. This report explores the benefits and pitfalls of the application of Mung bean (Mb) and V1 nuclease, which attack specifically single- and double-stranded regions of nucleic acids, as possible structural probes to be employed in combination with MS detection. Both enzymes were found capable of operating in ammonium-based solutions that are preferred for high-resolution analysis by direct infusion electrospray ionization (ESI). Sequence analysis by tandem mass spectrometry (MS/MS) was performed to confirm mapping assignments and to resolve possible ambiguities arising from the concomitant formation of isobaric products with identical base composition and different sequences. The observed products grouped together into ladder-type series that facilitated their assignment to unique regions of the substrate, but revealed also a certain level of uncertainty in identifying the boundaries between paired and unpaired regions. Various experimental factors that are known to stabilize nucleic acid structure, such as higher ionic strength, presence of Mg(II), etc., increased the accuracy of cleavage information, but did not completely eliminate deviations from expected results. These observations suggest extreme caution in interpreting the results afforded by these types of reagents. Regardless of the analytical platform of choice, the results highlighted the need to repeat probing experiments under the most diverse possible conditions to recognize potential artifacts and to increase the level of confidence in the observed structural information.

Trans-regulation of RNA-binding protein motifs by microRNA.

The wide array of vital functions that RNA performs is dependent on its ability to dynamically fold into different structures in response to intracellular and extracellular changes. RNA-binding proteins regulate much of this activity by targeting specific RNA structures or motifs. One of these structures, the 3-way RNA junction, is characteristically found in ribosomal RNA and results from the RNA folding in cis, to produce three separate helices that meet around a central unpaired region. Here we demonstrate that 3-way junctions can also form in trans as a result of the binding of microRNAs in an unconventional manner with mRNA by splinting two non-contiguous regions together. This may be used to reinforce the base of a stem-loop motif being targeted by an RNA-binding protein. Trans interactions between non-coding RNA and mRNA may be used to control the post-transcriptional regulatory code and suggests a possible role for some of the recently described transcripts of unknown function expressed from the human genome.

On the nucleotide composition and structure of retroviral RNA genomes

Retroviral RNA genomes display a rich variety in their nucleotide composition. For instance, the single-stranded RNA genome of human T cell leukemia virus (HTLV-1) is C-rich and G-poor and that of the human immunodeficiency virus (HIV-1) is A-rich and C-poor. Animal retroviruses add further variation to this unexplained, but many times remarkable virus-specific property. We previously described that the nucleotide bias is even more extreme in the unpaired regions of the structured HIV-1 RNA genome, which has been probed by SHAPE technology. We now document that the same trend is apparent for the MFold-predicted RNA structure of HIV-1 RNA and subsequently investigated the predicted structures of the RNA genomes of other retroviruses. We conclude that all virus-specific signatures are enhanced for the unpaired nucleotides in the RNA genome. Consequently, the differences in nucleotide count between the diverse human and animal retroviruses are further exposed in the single stranded genome regions. We used a skew analysis to visualize these striking differences in nucleotide usage. Evolutionary events responsible for these nucleotide signatures will be discussed.

Strand-specific Recognition of DNA Damages by XPD Provides Insights into Nucleotide Excision Repair Substrate Versatility

Recognition and removal of DNA damages is essential for cellular and organismal viability. Nucleotide excision repair (NER) is the sole mechanism in humans for the repair of carcinogenic UV irradiation-induced photoproducts in the DNA, such as cyclobutane pyrimidine dimers. The broad substrate versatility of NER further includes, among others, various bulky DNA adducts. It has been proposed that the 5'-3' helicase XPD (xeroderma pigmentosum group D) protein plays a decisive role in damage verification. However, despite recent advances such as the identification of a DNA-binding channel and central pore in the protein, through which the DNA is threaded, as well as a dedicated lesion recognition pocket near the pore, the exact process of target site recognition and verification in eukaryotic NER still remained elusive. Our single molecule analysis by atomic force microscopy reveals for the first time that XPD utilizes different recognition strategies to verify structurally diverse lesions. Bulky fluorescein damage is preferentially detected on the translocated strand, whereas the opposite strand preference is observed for a cyclobutane pyrimidine dimer lesion. Both states, however, lead to similar conformational changes in the resulting specific complexes, indicating a merge to a "final" verification state, which may then trigger the recruitment of further NER proteins.

Adjacent single-stranded regions mediate processing of tRNA precursors by RNase E direct entry.

The RNase E family is renowned for being central to the processing and decay of all types of RNA in many species of bacteria, as well as providing the first examples of endonucleases that can recognize 5′-monophosphorylated ends thereby increasing the efficiency of cleavage. However, there is increasing evidence that some transcripts can be cleaved efficiently by Escherichia coli RNase E via direct entry, i.e. in the absence of the recognition of a 5′-monophosphorylated end. Here, we provide biochemical evidence that direct entry is central to the processing of transfer RNA (tRNA) in E. coli, one of the core functions of RNase E, and show that it is mediated by specific unpaired regions that are adjacent, but not contiguous to segments cleaved by RNase E. In addition, we find that direct entry at a site on the 5′ side of a tRNA precursor triggers a series of 5′-monophosphate-dependent cleavages. Consistent with a major role for direct entry in tRNA processing, we provide additional evidence that a 5′-monophosphate is not required to activate the catalysis step in cleavage. Other examples of tRNA precursors processed via direct entry are also provided. Thus, it appears increasingly that direct entry by RNase E has a major role in bacterial RNA metabolism.

The Core Microprocessor Component DiGeorge Syndrome Critical Region 8 (DGCR8) Is a Nonspecific RNA-binding Protein

MicroRNA (miRNA) biogenesis follows a conserved succession of processing steps, beginning with the recognition and liberation of an miRNA-containing precursor miRNA hairpin from a large primary miRNA transcript (pri-miRNA) by the Microprocessor, which consists of the nuclear RNase III Drosha and the double-stranded RNA-binding domain protein DGCR8 (DiGeorge syndrome critical region protein 8). Current models suggest that specific recognition is driven by DGCR8 detection of single-stranded elements of the pri-miRNA stem-loop followed by Drosha recruitment and pri-miRNA cleavage. Because countless RNA transcripts feature single-stranded-dsRNA junctions and DGCR8 can bind hundreds of mRNAs, we explored correlations between RNA binding properties of DGCR8 and specific pri-miRNA substrate processing. We found that DGCR8 bound single-stranded, double-stranded, and random hairpin transcripts with similar affinity. Further investigation of DGCR8/pri-mir-16 interactions by NMR detected intermediate exchange regimes over a wide range of stoichiometric ratios. Diffusion analysis of DGCR8/pri-mir-16 interactions by pulsed field gradient NMR lent further support to dynamic complex formation involving free components in exchange with complexes of varying stoichiometry, although in vitro processing assays showed exclusive cleavage of pri-mir-16 variants bearing single-stranded flanking regions. Our results indicate that DGCR8 binds RNA nonspecifically. Therefore, a sequential model of DGCR8 recognition followed by Drosha recruitment is unlikely. Known RNA substrate requirements are broad and include 70-nucleotide hairpins with unpaired flanking regions. Thus, specific RNA processing is likely facilitated by preformed DGCR8-Drosha heterodimers that can discriminate between authentic substrates and other hairpins.

Human endonuclease V is a ribonuclease specific for inosine-containing RNA

Deamination of DNA bases can create missense mutations predisposing humans to cancer and also interfere with other basic molecular genetic processes; this deamination generates deoxyinosine from deoxyadenosine. In Escherichia coli, the highly conserved endonuclease V is involved in alternative excision repair that removes deoxyinosine from DNA. However, its exact activities and roles in humans are unknown. Here we characterize the FLJ35220 protein, the human homologue of E. coli endonuclease V, hEndoV as a ribonuclease specific for inosine-containing RNA. hEndoV preferentially binds to RNA and efficiently hydrolyses the second phosphodiester bond located 3′ to the inosine in unpaired inosine-containing ssRNA regions in dsRNA. It localizes to the cytoplasm in cells. The ribonuclease activity is promoted by Tudor staphylococcal nuclease and detected on inosine-containing dsRNA created by the action of adenosine deaminases acting on RNA. These results demonstrate that hEndoV controls the fate of inosine-containing RNA in humans.

SYCP3-like X-linked 2 is expressed in meiotic germ cells and interacts with synaptonemal complex central element protein 2 and histone acetyltransferase TIP60

Meiosis is the process by which diploid germ cells produce haploid gametes. A key event is the formation of the synaptonemal complex. In the pachytene stage, the unpaired regions of X and Y chromosomes form a specialized structure, the XY body, within which gene expression is mostly silenced. In the present study, we showed that SYCP3-like X-linked 2 (SLX2, 1700013H16Rik), a novel member of XLR (X-linked Lymphocyte-Regulated) family, was specifically expressed in meiotic germ cells. In the spermatocyte SLX2 was distributed in the nucleus of germ cells at the preleptotene, leptotene and zygotene stages and is then restricted to the XY body at the pachytene stage. This localization change is coincident with that of phosphorylated histone H2AX (γH2AX), a well-known component of the sex body. Through yeast two-hybrid screening and coimmunoprecipitation assays, we demonstrated that SLX2 interacts with synaptonemal complex central element protein 2 (SYCE2), an important component of synaptonemal complex, and histone acetyltransferase TIP60, which has been implicated in remodeling phospho-H2AX-containing nucleosomes at sites of DNA damage. These results suggest that SLX2 might be involved in DNA recombination, synaptonemal complex formation as well as sex body maintenance during meiosis.

Fitness Costs of Synonymous Mutations in the rpsT Gene Can Be Compensated by Restoring mRNA Base Pairing

We previously reported that the distribution of fitness effects for non-synonymous and synonymous mutations in Salmonella typhimurium ribosomal proteins S20 and L1 are similar, suggesting that fitness constraints are present at the level of mRNA. Here we explore the hypothesis that synonymous mutations confer their fitness-reducing effect by alterating the secondary structure of the mRNA. To this end, we constructed a set of synonymous substitutions in the rpsT gene, encoding ribosomal protein S20, that are located in predicted paired regions in the mRNA and measured their effect on bacterial fitness. Our results show that for 3/9 cases tested, the reduced fitness conferred by a synonymous mutation could be fully or partly restored by introducing a second synonymous substitution that restore base pairing in a mRNA stem. In addition, random mutations in predicted paired regions had larger fitness effects than those in unpaired regions. Finally, we did not observe any correlation between fitness effects of the synonymous mutations and their rarity. These results suggest that for ribosomal protein S20, the deleterious effects of synonymous mutations are not generally due to codon usage effects, but that mRNA secondary structure is a major fitness constraint.

Identification of Small RNAs Associated with Meiotic Silencing by Unpaired DNA

In Neurospora crassa, unpaired genes are silenced by a mechanism called meiotic silencing by unpaired DNA (MSUD). Although some RNA interference proteins are necessary for this process, its requirement of small RNAs has yet to be formally established. Here we report the characterization of small RNAs targeting an unpaired region, using Illumina sequencing.

Reverse Gyrase Transiently Unwinds Double-Stranded DNA in an ATP-Dependent Reaction

Reverse gyrase is a unique DNA topoisomerase that catalyzes the introduction of positive supercoils into DNA in an ATP-dependent reaction. It consists of a helicase domain that functionally cooperates with a topoisomerase domain. Different models for the catalytic mechanism of reverse gyrase that predict a central role of the helicase domain have been put forward. The helicase domain acts as a nucleotide-dependent conformational switch that alternates between open and closed states with different affinities for single- and double-stranded DNA. It has been suggested that the helicase domain can unwind double-stranded regions, but helicase activity has not been demonstrated as yet. Here, we show that the isolated helicase domain and full-length reverse gyrase can transiently unwind double-stranded regions in an ATP-dependent reaction. The latch region of reverse gyrase, an insertion into the helicase domain, is required for DNA supercoiling. Strikingly, the helicase domain lacking the latch cannot unwind DNA, linking unwinding to DNA supercoiling. The unwinding activity may provide and stabilize the single-stranded regions required for strand passage by the topoisomerase domain, either de novo or by expanding already existing unpaired regions that may form at high temperatures.

The Genome Sequence of an Isolate of Indian Citrus Ringspot Virus Infecting the Sweet Orange in India

Whole-genome sequencing of an isolate of Mandarivirus infecting the sweet orange [Citrus sinensis (L) Blanco] in the western part of India (Pune) was done. The single-stranded positive-sense RNA genome of Indian citrus ringspot virus (ICRSV) Pune has 7,560 nucleotides (nt), excluding a poly(A) tail, comprised of 27.98% (2,115 nt) A, 32.12% (2,428 nt) C, 19.68% (1,488 nt) G, and 20.22% (1,529 nt) T residues. The genome, organized into six open reading frames (ORFs), shares 97.7% sequence identity with the complete genome of the ICRSV K1 isolate (AF406744.1) infecting the kinnow (Citrus reticulate Blanco, a hybrid between King and Willow mandarins) in north India. The ICRSV Pune genome formed a complex secondary structure with a large number of unpaired cytosine-rich regions, and recombination analysis highlighted potential recombination in the ICRSV genome.

Monocularaly unpaired regions do not resist suppression in absence of an explicit occluder.

Monocularaly unpaired regions do not resist suppression in absence of an explicit occluder.