GTfold: Enabling parallel RNA secondary structure prediction on multi-core desktops

ABSTRACT: BACKGROUND: Accurate and efficient RNA secondary structure prediction remains an important open problem in computational molecular biology. Historically, advances in computing technology have enabled faster and more accurate RNA secondary structure predictions. Previous parallelized prediction programs achieved significant improvements in runtime, but their implementations were not portable from niche high-performance computers or easily accessible to most RNA researchers. With the increasing prevalence of multi-core desktop machines, a new parallel prediction program is needed to take full advantage of today's computing technology. FINDINGS: We present here the first implementation of RNA secondary structure prediction by thermodynamic optimization for modern multi-core computers. We show that GTfold predicts secondary structure in less time than UNAfold and RNAfold, without sacrificing accuracy, on machines with four or more cores. CONCLUSIONS: GTfold supports advances in RNA structural biology by reducing the timescales for secondary structure prediction. The difference will be particularly valuable to researchers working with lengthy RNA sequences, such as RNA viral genomes.     

AreA controls nitrogen source utilisation during both growth programs of the dimorphic fungus Penicillium marneffei

The opportunistic pathogen Penicillium marneffei displays a temperature-dependent dimorphic switching program with saprophytic hyphal growth at 25 °C and yeast growth at 37 °C. The areA gene of P. marneffei has been isolated and found to be required for the utilisation of nonpreferred nitrogen sources during both growth programs of P. marneffei, albeit to differing degrees. Based on this functional characterisation and high degree of sequence conservation with other fungal GATA factors, P. marneffei areA represents an orthologue of Aspergillus nidulans areA and Neurospora crassa NIT2. Based on this study it is proposed that AreA is likely to contribute to the pathogenicity of P. marneffei by facilitating growth in the host environment and regulating the expression of potential virulence factors such as extracellular proteases.     

Full modeling versus summarizing gene-tree uncertainty: Method choice and species-tree accuracy

With the proliferation of species-tree methods, empiricists now have to confront the daunting task of method choice. Such decisions might be made based on theoretical considerations alone. However, the messiness of real data means that theoretical ideals may not hold in practice (e.g., with convergence of complicated MCMC algorithms and computational times that limit analyses to small data sets). On the other hand, simplifying assumptions made by some approaches may compromise the accuracy of species-tree estimates. Here we examine the purported tradeoff between accuracy and computational simplicity for species-tree analysis, focusing on the different ways the approaches treat gene-tree uncertainty. By considering a diversity of species trees, as well as different sampling designs and total sampling efforts, we not only compare the accuracy of species-tree estimates across methods, but we also partition the variation in accuracy across factors to identify their relative importance. This analysis shows that although the method of analysis affects accuracy, other factors - namely, the history of species divergence and aspects of the sampling design - have a larger impact. Despite a full modeling of gene tree uncertainty (e.g., using a Bayesian framework), species-tree estimates may not be accurate, particularly for recent diversification histories. Nevertheless, we demonstrate how factors within the control of the empirical investigator (e.g., decisions about sampling) improve the accuracy of species tree estimates, and more so than the method of analysis. Lastly, with much of the attention on species-tree analyses focused on the discord among loci arising from the coalescent, this work also highlights a previously overlooked key determinant of species-tree accuracy for recent divergences - the level of genetic variation at a locus, which has important implications for improving species-tree estimates in practice.     

Gene Expression in the Spinal Cord in Female Lewis Rats with Experimental Autoimmune Encephalomyelitis Induced with Myelin Basic Protein

Background

Experimental autoimmune encephalomyelitis (EAE), the best available model of multiple sclerosis, can be induced in different animal strains using immunization with central nervous system antigens. EAE is associated with inflammation and demyelination of the nervous system. Micro-array can be used to investigate gene expression and biological pathways that are altered during disease. There are few studies of the changes in gene expression in EAE, and these have mostly been done in a chronic mouse EAE model. EAE induced in the Lewis with myelin basic protein (MBP-EAE) is well characterised, making it an ideal candidate for the analysis of gene expression in this disease model.

Methodology/Principal Findings

MBP-EAE was induced in female Lewis rats by inoculation with MBP and adjuvants. Total RNA was extracted from the spinal cords and used for micro-array analysis using AffimetrixGeneChip Rat Exon 1.0 ST Arrays. Gene expression in the spinal cords was compared between healthy female rats and female rats with MBP-EAE. Gene expression in the spinal cord of rats with MBP-EAE differed from that in the spinal cord of normal rats, and there was regulation of pathways involved with immune function and nervous system function. For selected genes the change in expression was confirmed with real-time PCR.

Conclusions/Significance

EAE leads to modulation of gene expression in the spinal cord. We have identified the genes that are most significantly regulated in MBP-EAE in the Lewis rat and produced a profile of gene expression in the spinal cord at the peak of disease.     

Evaluating the Ribosomal Internal Transcribed Spacer (ITS) as a Candidate Dinoflagellate Barcode Marker

Background

DNA barcoding offers an efficient way to determine species identification and to measure biodiversity. For dinoflagellates, an ancient alveolate group of about 2000 described extant species, DNA barcoding studies have revealed large amounts of unrecognized species diversity, most of which is not represented in culture collections. To date, two mitochondrial gene markers, Cytochrome Oxidase I (COI) and Cytochrome b oxidase (COB), have been used to assess DNA barcoding in dinoflagellates, and both failed to amplify all taxa and suffered from low resolution. Nevertheless, both genes yielded many examples of morphospecies showing cryptic speciation and morphologically distinct named species being genetically similar, highlighting the need for a common marker. For example, a large number of cultured Symbiodinium strains have neither taxonomic identification, nor a common measure of diversity that can be used to compare this genus to other dinoflagellates.

Methodology/Principal Findings

The purpose of this study was to evaluate the Internal Transcribed Spacer units 1 and 2 (ITS) of the rDNA operon, as a high resolution marker for distinguishing species dinoflagellates in culture. In our study, from 78 different species, the ITS barcode clearly differentiated species from genera and could identify 96% of strains to a known species or sub-genus grouping. 8.3% showed evidence of being cryptic species. A quarter of strains identified had no previous species identification. The greatest levels of hidden biodiversity came from Scrippsiella and the Pfiesteriaceae family, whilst Heterocapsa strains showed a high level of mismatch to their given species name.

Conclusions/Significance

The ITS marker was successful in confirming species, revealing hidden diversity in culture collections. This marker, however, may have limited use for environmental barcoding due to paralogues, the potential for unidentifiable chimaeras and priming across taxa. In these cases ITS would serve well in combination with other markers or for specific taxon studies.     

A systematic approach to pair secretory cargo receptors with their cargo suggests a mechanism for cargo selection by Erv14

The endoplasmic reticulum (ER) is the site of synthesis of secreted and membrane proteins. To exit the ER, proteins are packaged into COPII vesicles through direct interaction with the COPII coat or aided by specific cargo receptors. Despite the fundamental role of such cargo receptors in protein traffic, only a few have been identified; their cargo spectrum is unknown and the signals they recognize remain poorly understood. We present here an approach we term "PAIRS" (pairing analysis of cargo receptors), which combines systematic genetic manipulations of yeast with automated microscopy screening, to map the spectrum of cargo for a known receptor or to uncover a novel receptor for a particular cargo. Using PAIRS we followed the fate of ～150 cargos on the background of mutations in nine putative cargo receptors and identified novel cargo for most of these receptors. Deletion of the Erv14 cargo receptor affected the widest range of cargo. Erv14 substrates have a wide array of functions and structures; however, they are all membrane-spanning proteins of the late secretory pathway or plasma membrane. Proteins residing in these organelles have longer transmembrane domains (TMDs). Detailed examination of one cargo supported the hypothesis that Erv14 dependency reflects the length rather than the sequence of the TMD. The PAIRS approach allowed us to uncover new cargo for known cargo receptors and to obtain an unbiased look at specificity in cargo selection. Obtaining the spectrum of cargo for a cargo receptor allows a novel perspective on its mode of action. The rules that appear to guide Erv14 substrate recognition suggest that sorting of membrane proteins at multiple points in the secretory pathway could depend on the physical properties of TMDs. Such a mechanism would allow diverse proteins to utilize a few receptors without the constraints of evolving location-specific sorting motifs.     

High mitochondrial and nuclear genetic diversity in one of the world’s most endangered seabirds,the Chatham Island Taiko (<Emphasis Type="Italic">Pterodroma magentae</Emphasis>)

Interpreting the levels of genetic diversity in organisms with diverse life and population histories can be difficult. The processes and mechanisms regulating this diversity are complex and still poorly understood. However, endangered species typically have low genetic variation as a consequence of the effects of genetic drift in small populations. In this study we examine genetic variation in the critically endangered Chatham Island Taiko (Tchaik, Pterodroma magentae), one of the world’s rarest seabirds. The Taiko has a very small population size of between 120 and 150 individuals, including just 8–15 breeding pairs. We report surprisingly high mitochondrial and nuclear genetic diversity in this critically endangered long-lived species. We hypothesise that the present Taiko population has retained a significant proportion of its past genetic diversity. However, it is also possible that undiscovered birds are breeding in unknown areas, which could increase the population size estimate. Importantly, from a conservation perspective, we show that the high level of variation is unlikely to be maintained in the future since chicks currently being born have only a limited number of the mitochondrial DNA haplotypes found in adults. Reduced genetic variation will mean that our ability to infer past events and the population history of Taiko using genetics could soon be lost and the power to determine, for example, parentage and other close order relationships will be diminished. Therefore, the maintenance of genetic diversity in future generations is an important consideration for conservation management of the Taiko.     

Mutations causing DOK7 congenital myasthenia ablate functional motifs in Dok-7

Dok-7 is a cytoplasmic activator of muscle-specific receptor-tyrosine kinase (MuSK). Both Dok-7 and MuSK are required for neuromuscular synaptogenesis. Mutations in DOK7 underlie a congenital myasthenic syndrome (CMS) associated with small and simplified neuromuscular synapses likely due to impaired Dok-7/MuSK signaling. The overwhelming majority of patients with DOK7 CMS have at least one allele with a frameshift mutation that causes a truncation in the COOH-terminal region of Dok-7 and affects MuSK activation. Dok-7 has pleckstrin homology (PH) and phosphotyrosine binding (PTB) domains in the NH2-terminal moiety, both of which are indispensable for MuSK activation in myotubes, but little is known about additional functional elements. Here, we identify a chromosome region maintenance 1-dependent nuclear export signal (NES) in the COOH-terminal moiety and demonstrate that the NES-mediated cytoplasmic location of Dok-7 is essential for regulating the interaction with MuSK in myotubes. The NH2-terminal PH domain is responsible for the nuclear import of Dok-7. We also show that the Src homology 2 target motifs in the COOH-terminal moiety of Dok-7 are active and crucial for MuSK activation in myotubes. In addition, CMS-associated missense mutations found in the PH or PTB domain inactivate Dok-7. Together, these findings demonstrate that, in addition to the NH2-terminal PH and PTB domains, the COOH-terminal NES and Src homology 2 target motifs play key roles in Dok-7/MuSK signaling for neuromuscular synaptogenesis. Ablation or disruption of these functional elements in Dok-7 probably underlies the neuromuscular junction synaptopathy observed in DOK7 CMS.     

Reproduction and development of Russian wheat aphid biotype 2 on crested wheatgrass, intermediate wheatgrass, and susceptible and resistant wheat

The Russian wheat aphid, Diuraphis noxia (Kurdjumov), is an economically important pest of small grains. Since its introduction into North America in 2003, Russian wheat aphid Biotype 2 has been found to be virulent to all commercially available winter wheat, Triticum aestivum L., cultivars. Our goal was to examine differences in Russian wheat aphid reproduction and development on a variety of plant hosts to gain information about 1) potential alternate host refuges, 2) selective host pressures on Russian wheat aphid genetic variation, and 3) general population dynamics of Russian wheat aphid Biotype 2. We studied host quality of two wheatgrasses (crested wheatgrass, Agropyron cristatum [L.] Gaertn., and intermediate wheatgrass, Agropyron intermedium [Host] Beauvoir) and two types of winter wheat (T. aestivum, one Biotype 2 susceptible wheat, 'Custer' and one biotype 2 resistant wheat, STARS02RWA2414-11). The susceptible wheat had the highest intrinsic rate of increase, greatest longevity and greatest fecundity of the four host studied. Crested wheatgrass and the resistant wheat showed similar growth rates. Intermediate wheatgrass had the lowest intrinsic rate of increase and lowest fecundity of all tested hosts.