ScanFold: an approach for genome-wide discovery of local RNA structural elements—applications to Zika virus and HIV

In addition to encoding RNA primary structures, genomes also encode RNA secondary and tertiary structures that play roles in gene regulation and, in the case of RNA viruses, genome replication. Methods for the identification of functional RNA structures in genomes typically rely on scanning analysis windows, where multiple partially-overlapping windows are used to predict RNA structures and folding metrics to deduce regions likely to form functional structure. Separate structural models are produced for each window, where the step size can greatly affect the returned model. This makes deducing unique local structures challenging, as the same nucleotides in each window can be alternatively base paired. We are presenting here a new approach where all base pairs from analysis windows are considered and weighted by favorable folding. This results in unique base pairing throughout the genome and the generation of local regions/structures that can be ranked by their propensity to form unusually thermodynamically stable folds. We applied this approach to the Zika virus (ZIKV) and HIV-1 genomes. ZIKV is linked to a variety of neurological ailments including microcephaly and Guillain–Barré syndrome and its (+)-sense RNA genome encodes two, previously described, functionally essential structured RNA regions. HIV, the cause of AIDS, contains multiple functional RNA motifs in its genome, which have been extensively studied. Our approach is able to successfully identify and model the structures of known functional motifs in both viruses, while also finding additional regions likely to form functional structures. All data have been archived at the RNAStructuromeDB (www.structurome.bb.iastate.edu), a repository of RNA folding data for humans and their pathogens

Constructive and destructive use of compilers in elliptic curve cryptography

Although cryptographic software implementation is often performed by expert programmers, the range of performance and security driven options, as well as more mundane software engineering issues, still make it a challenge. The use of domain specific language and compiler techniques to assist in description and optimisation of cryptographic software is an interesting research challenge. In this paper we investigate two aspects of such techniques, focusing on Elliptic Curve Cryptography (ECC) in particular. Our constructive results show that a suitable language allows description of ECC based software in a manner close to the original mathematics; the corresponding compiler allows automatic production of an executable whose performance is competitive with that of a hand-optimised implementation. In contrast, we study the worrying potential for naïve compiler driven optimisation to render cryptographic software insecure. Both aspects of our work are set within the context of CACE, an ongoing EU funded project on this general topic

Plasmodium falciparum translational machinery condones polyadenosine repeats

Plasmodium falciparum is a causative agent of human malaria. Sixty percent of mRNAs from its extremely AT-rich (81%) genome harbor long polyadenosine (polyA) runs within their ORFs, distinguishing the parasite from its hosts and other sequenced organisms. Recent studies indicate polyA runs cause ribosome stalling and frameshifting, triggering mRNA surveillance pathways and attenuating protein synthesis. Here, we show that P. falciparum is an exception to this rule. We demonstrate that both endogenous genes and reporter sequences containing long polyA runs are efficiently and accurately translated in P. falciparum cells. We show that polyA runs do not elicit any response from No Go Decay (NGD) or result in the production of frameshifted proteins. This is in stark contrast to what we observe in human cells or T. thermophila, an organism with similar AT-content. Finally, using stalling reporters we show that Plasmodium cells evolved not to have a fully functional NGD pathway

EST analysis in Ginkgo biloba: an assessment of conserved developmental regulators and gymnosperm specific genes

BACKGROUND: Ginkgo biloba L. is the only surviving member of one of the oldest living seed plant groups with medicinal, spiritual and horticultural importance worldwide. As an evolutionary relic, it displays many characters found in the early, extinct seed plants and extant cycads. To establish a molecular base to understand the evolution of seeds and pollen, we created a cDNA library and EST dataset from the reproductive structures of male (microsporangiate), female (megasporangiate), and vegetative organs (leaves) of Ginkgo biloba. RESULTS: RNA from newly emerged male and female reproductive organs and immature leaves was used to create three distinct cDNA libraries from which 6,434 ESTs were generated. These 6,434 ESTs from Ginkgo biloba were clustered into 3,830 unigenes. A comparison of our Ginkgo unigene set against the fully annotated genomes of rice and Arabidopsis, and all available ESTs in Genbank revealed that 256 Ginkgo unigenes match only genes among the gymnosperms and non-seed plants – many with multiple matches to genes in non-angiosperm plants. Conversely, another group of unigenes in Gingko had highly significant homology to transcription factors in angiosperms involved in development, including MADS box genes as well as post-transcriptional regulators. Several of the conserved developmental genes found in Ginkgo had top BLAST homology to cycad genes. We also note here the presence of ESTs in G. biloba similar to genes that to date have only been found in gymnosperms and an additional 22 Ginkgo genes common only to genes from cycads. CONCLUSION: Our analysis of an EST dataset from G. biloba revealed genes potentially unique to gymnosperms. Many of these genes showed homology to fully sequenced clones from our cycad EST dataset found in common only with gymnosperms. Other Ginkgo ESTs are similar to developmental regulators in higher plants. This work sets the stage for future studies on Ginkgo to better understand seed and pollen evolution, and to resolve the ambiguous phylogenetic relationship of G. biloba among the gymnosperms

Class I major histocompatibility complexes loaded by a periodate trigger

Class I major histocompatibility complexes (MHCs) present peptide ligands on the cell surface for recognition by appropriate cytotoxic T cells. The unstable nature of unliganded MHC necessitates the production of recombinant class I complexes through in vitro refolding reactions in the presence of an added excess of peptides. This strategy is not amenable to high-throughput production of vast collections of class I complexes. To address this issue, we recently designed photocaged MHC ligands that can be cleaved by a UV light trigger in the MHC bound state under conditions that do not affect the integrity of the MHC structure. The results obtained with photocaged MHC ligands demonstrate that conditional MHC ligands can form a generally applicable concept for the creation of defined peptide−MHCs. However, the use of UV exposure to mediate ligand exchange is unsuited for a number of applications, due to the lack of UV penetration through cell culture systems and due to the transfer of heat upon UV irradiation, which can induce evaporation. To overcome these limitations, here, we provide proof-of-concept for the generation of defined peptide−MHCs by chemical trigger-induced ligand exchange. The crystal structure of the MHC with the novel chemosensitive ligand showcases that the ligand occupies the expected binding site, in a conformation where the hydroxyl groups should be reactive to periodate. We proceed to validate this technology by producing peptide−MHCs that can be used for T cell detection. The methodology that we describe here should allow loading of MHCs with defined peptides in cell culture devices, thereby permitting antigen-specific T cell expansion and purification for cell therapy. In addition, this technology will be useful to develop miniaturized assay systems for performing high-throughput screens for natural and unnatural MHC ligands

Non-gray rotating stellar models and the evolutionary history of the Orion Nebular Cluster

Rotational evolution in the pre-main sequence (PMS) is described with new sets of PMS evolutionary tracks including rotation, non-gray boundary conditions (BCs) and either low (LCE) or high convection efficiency (HCE). Using observational data and our theoretical predictions, we aim at constraining 1) the differences obtained for the rotational evolution of stars within the ONC by means of these different sets of models; 2) the initial angular momentum of low mass stars, by means of their templates in the ONC. We discuss the reliability of current stellar models for the PMS. While the 2D radiation hydrodynamic simulations predict HCE in PMS, semi-empirical calibrations either seem to require that convection is less efficient in PMS than in the following MS phase or are still contradictory. We derive stellar masses and ages for the ONC by using both LCE and HCE. The resulting mass distribution for the bulk of the ONC population is in the range 0.2$-$0.3 {\msun} for our non-gray models and in the range 0.1$-$0.3{\msun} for models having gray BCs. In agreement with Herbst et al. (2002) we find that a large percentage ($\sim$70%) of low-mass stars (M\simlt 0.5{\msun} for LCE; M\simlt0.35{\msun} for HCE) in the ONC appears to be fast rotators (P$<$4days). Three possibilities are open: 1) $\sim$70% of the ONC low mass stars lose their disk at early evolutionary phases; 2)their locking period is shorter; 3) the period evolution is linked to a different morphology of the magnetic fields of the two groups of stars. We also estimate the range of initial angular momentum consistent with the observed periods. The comparisons made indicate that a second parameter is needed to describe convection in the PMS, possibly related to the structural effect of a dynamo magnetic field.Comment: 17 pages, 11 figure

The effect of time since measles vaccination and age at first dose on measles vaccine effectiveness - A systematic review.

BACKGROUND: In settings where measles has been eliminated, vaccine-derived immunity may in theory wane more rapidly due to a lack of immune boosting by circulating measles virus. We aimed to assess whether measles vaccine effectiveness (VE) waned over time, and if so, whether differentially in measles-eliminated and measles-endemic settings. METHODS: We performed a systematic literature review of studies that reported VE and time since vaccination with measles-containing vaccine (MCV). We extracted information on case definition (clinical symptoms and/or laboratory diagnosis), method of vaccination status ascertainment (medical record or vaccine registry), as well as any biases which may have arisen from cold chain issues and a lack of an age at first dose of MCV. We then used linear regression to evaluate VE as a function of age at first dose of MCV and time since MCV. RESULTS: After screening 14,782 citations, we identified three full-text articles from measles-eliminated settings and 33 articles from measles-endemic settings. In elimination settings, two-dose VE estimates increased as age at first dose of MCV increased and decreased as time since MCV increased; however, the small number of studies available limited interpretation. In measles-endemic settings, one-dose VE increased by 1.5% (95% CI 0.5, 2.5) for every month increase in age at first dose of MCV. We found no evidence of waning VE in endemic settings. CONCLUSIONS: The paucity of data from measles-eliminated settings indicates that additional studies and approaches (such as studies using proxies including laboratory correlates of protection) are needed to answer the question of whether VE in measles-eliminated settings wanes. Age at first dose of MCV was the most important factor in determining VE. More VE studies need to be conducted in elimination settings, and standards should be developed for information collected and reported in such studies

Influenza A Virus Coding Regions Exhibit Host-Specific Global Ordered RNA Structure

Influenza A is a significant public health threat, partially because of its capacity to readily exchange gene segments between different host species to form novel pandemic strains. An understanding of the fundamental factors providing species barriers between different influenza hosts would facilitate identification of strains capable of leading to pandemic outbreaks and could also inform vaccine development. Here, we describe the difference in predicted RNA secondary structure stability that exists between avian, swine and human coding regions. The results predict that global ordered RNA structure exists in influenza A segments 1, 5, 7 and 8, and that ranges of free energies for secondary structure formation differ between host strains. The predicted free energy distributions for strains from avian, swine, and human species suggest criteria for segment reassortment and strains that might be ideal candidates for viral attenuation and vaccine development

Isoenergetic penta- and hexanucleotide microarray probing and chemical mapping provide a secondary structure model for an RNA element orchestrating R2 retrotransposon protein function

LNA (locked nucleic acids, i.e. oligonucleotides with a methyl bridge between the 2′ oxygen and 4′ carbon of ribose) and 2,6-diaminopurine were incorporated into 2′-O-methyl RNA pentamer and hexamer probes to make a microarray that binds unpaired RNA approximately isoenergetically. That is, binding is roughly independent of target sequence if target is unfolded. The isoenergetic binding and short probe length simplify interpretation of binding to a structured RNA to provide insight into target RNA secondary structure. Microarray binding and chemical mapping were used to probe the secondary structure of a 323 nt segment of the 5′ coding region of the R2 retrotransposon from Bombyx mori (R2Bm 5′ RNA). This R2Bm 5′ RNA orchestrates functioning of the R2 protein responsible for cleaving the second strand of DNA during insertion of the R2 sequence into the genome. The experimental results were used as constraints in a free energy minimization algorithm to provide an initial model for the secondary structure of the R2Bm 5′ RNA