OMMA enables population-scale analysis of complex genomic features and phylogenomic relationships from nanochannel-based optical maps.

BackgroundOptical mapping is an emerging technology that complements sequencing-based methods in genome analysis. It is widely used in improving genome assemblies and detecting structural variations by providing information over much longer (up to 1 Mb) reads. Current standards in optical mapping analysis involve assembling optical maps into contigs and aligning them to a reference, which is limited to pairwise comparison and becomes bias-prone when analyzing multiple samples.FindingsWe present a new method, OMMA, that extends optical mapping to the study of complex genomic features by simultaneously interrogating optical maps across many samples in a reference-independent manner. OMMA captures and characterizes complex genomic features, e.g., multiple haplotypes, copy number variations, and subtelomeric structures when applied to 154 human samples across the 26 populations sequenced in the 1000 Genomes Project. For small genomes such as pathogenic bacteria, OMMA accurately reconstructs the phylogenomic relationships and identifies functional elements across 21 Acinetobacter baumannii strains.ConclusionsWith the increasing data throughput of optical mapping system, the use of this technology in comparative genome analysis across many samples will become feasible. OMMA is a timely solution that can address such computational need. The OMMA software is available at https://github.com/TF-Chan-Lab/OMTools

CoRAL: Predicting Non-Coding RNAs from Small RNA-Sequencing Data

The surprising observation that virtually the entire human genome is transcribed means we know little about the function of many emerging classes of RNAs, except their astounding diversities. Traditional RNA function prediction methods rely on sequence or alignment information, which are limited in their abilities to classify the various collections of non-coding RNAs (ncRNAs). To address this, we developed Classification of RNAs by Analysis of Length (CoRAL), a machine learning-based approach for classification of RNA molecules. CoRAL uses biologically interpretable features including fragment length and cleavage specificity to distinguish between different ncRNA populations. We evaluated CoRAL using genome-wide small RNA sequencing data sets from four human tissue types and were able to classify six different types of RNAs with ∼80% cross-validation accuracy. Analysis by CoRAL revealed that microRNAs, small nucleolar and transposon-derived RNAs are highly discernible and consistent across all human tissue types assessed, whereas long intergenic ncRNAs, small cytoplasmic RNAs and small nuclear RNAs show less consistent patterns. The ability to reliably annotate loci across tissue types demonstrates the potential of CoRAL to characterize ncRNAs using small RNA sequencing data in less well-characterized organisms

HAMR: High-Throughput Annotation of Modified Ribonucleotides

RNA is often altered post-transcriptionally by the covalent modification of particular nucleotides; these modifications are known to modulate the structure and activity of their host RNAs. The recent discovery that an RNA methyl-6 adenosine demethylase (FTO) is a risk gene in obesity has brought to light the significance of RNA modifications to human biology. These noncanonical nucleotides, when converted to cDNA in the course of RNA sequencing, can produce sequence patterns that are distinguishable from simple base-calling errors. To determine whether these modifications can be detected in RNA sequencing data, we developed a method that can not only locate these modifications transcriptome-wide with single nucleotide resolution, but can also differentiate between different classes of modifications. Using small RNA-seq data we were able to detect 92% of all known human tRNA modification sites that are predicted to affect RT activity. We also found that different modifications produce distinct patterns of cDNA sequence, allowing us to differentiate between two classes of adenosine and two classes of guanine modifications with 98% and 79% accuracy, respectively. To show the robustness of this method to sample preparation and sequencing methods, as well as to organismal diversity, we applied it to a publicly available yeast data set and achieved similar levels of accuracy. We also experimentally validated two novel and one known 3-methylcytosine (3mC) sites predicted by HAMR in human tRNAs. Researchers can now use our method to identify and characterize RNA modifications using only RNA-seq data, both retrospectively and when asking questions specifically about modified RNA

DASHR: Database of Small Human Noncoding RNAs

Small non-coding RNAs (sncRNAs) are highly abundant RNAs, typically long, that act as key regulators of diverse cellular processes. Although thousands of sncRNA genes are known to exist in the human genome, no single database provides searchable, unified annotation, and expression information for full sncRNA transcripts and mature RNA products derived from these larger RNAs. Here, we present the Database of small human noncoding RNAs (DASHR) . DASHR contains the most comprehensive information to date on human sncRNA genes and mature sncRNA products. DASHR provides a simple user interface for researchers to view sequence and secondary structure, compare expression levels, and evidence of specific processing across all sncRNA genes and mature sncRNA products in various human tissues. DASHR annotation and expression data covers all major classes of sncRNAs including microRNAs (miRNAs), Piwi-interacting (piRNAs), small nuclear, nucleolar, cytoplasmic (sn-, sno-, scRNAs, respectively), transfer (tRNAs), and ribosomal RNAs (rRNAs). Currently, DASHR (v1.0) integrates 187 smRNA high-throughput sequencing (smRNA-seq) datasets with over 2.5 billion reads and annotation data from multiple public sources. DASHR contains annotations for ~48,000 human sncRNA genes and mature sncRNA products, 82% of which are expressed in one of more of the curated tissues. DASHR is available at http://lisanwanglab.org/DASHR

Changes in the Transcriptome of Human Astrocytes Accompanying Oxidative Stress-Induced Senescence

Aging is a major risk factor for many neurodegenerative disorders. A key feature of aging biology that may underlie these diseases is cellular senescence. Senescent cells accumulate in tissues with age, undergo widespread changes in gene expression, and typically demonstrate altered, pro-inflammatory profiles. Astrocyte senescence has been implicated in neurodegenerative disease, and to better understand senescence-associated changes in astrocytes, we investigated changes in their transcriptome using RNA sequencing. Senescence was induced in human fetal astrocytes by transient oxidative stress. Brain-expressed genes, including those involved in neuronal development and differentiation, were downregulated in senescent astrocytes. Remarkably, several genes indicative of astrocytic responses to injury were also downregulated, including glial fibrillary acidic protein and genes involved in the processing and presentation of antigens by major histocompatiblity complex class II proteins, while pro-inflammatory genes were upregulated. Overall, our findings suggest that senescence-related changes in the function of astrocytes may impact the pathogenesis of age-related brain disorders

Phenylthiourea Specifically Reduces Zebrafish Eye Size

Phenylthiourea (PTU) is commonly used for inhibiting melanization of zebrafish embryos. In this study, the standard treatment with 0.2 mM PTU was demonstrated to specifically reduce eye size in larval fish starting at three days post-fertilization. This effect is likely the result of a reduction in retinal and lens size of PTU-treated eyes and is not related to melanization inhibition. This is because the eye size of tyr, a genetic mutant of tyrosinase whose activity is inhibited in PTU treatment, was not reduced. As PTU contains a thiocarbamide group which is presented in many goitrogens, suppressing thyroid hormone production is a possible mechanism by which PTU treatment may reduce eye size. Despite the fact that thyroxine level was found to be reduced in PTU-treated larvae, thyroid hormone supplements did not rescue the eye size reduction. Instead, treating embryos with six goitrogens, including inhibitors of thyroid peroxidase (TPO) and sodium-iodide symporter (NIS), suggested an alternative possibility. Specifically, three TPO inhibitors, including those that do not possess thiocarbamide, specifically reduced eye size; whereas none of the NIS inhibitors could elicit this effect. These observations indicate that TPO inhibition rather than a general suppression of thyroid hormone synthesis is likely the underlying cause of PTU-induced eye size reduction. Furthermore, the tissue-specific effect of PTU treatment might be mediated by an eye-specific TPO expression. Compared with treatment with other tyrosinase inhibitors or bleaching to remove melanization, PTU treatment remains the most effective approach. Thus, one should use caution when interpreting results that are obtained from PTU-treated embryos

Number of People Blind or Visually Impaired by Glaucoma Worldwide and in World Regions 1990 – 2010: A Meta-Analysis

Objective: To assess the number of individuals visually impaired or blind due to glaucoma and to examine regional differences and temporal changes in this parameter for the period from 1990 to 2012. Methods: As part of the Global Burden of Diseases (GBD) Study 2010, we performed a systematic literature review for the period from 1980 to 2012. We primarily identified 14,908 relevant manuscripts, out of which 243 high-quality, population-based studies remained after review by an expert panel that involved application of selection criteria that dwelt on population representativeness and clarity of visual acuity methods used. Sixty-six specified the proportion attributable to glaucoma. The software tool DisMod-MR (Disease Modeling–Metaregression) of the GBD was used to calculate fraction of vision impairment due to glaucoma. Results: In 2010, 2.1 million (95% Uncertainty Interval (UI):1.9,2.6) people were blind, and 4.2 (95% UI:3.7,5.8) million were visually impaired due to glaucoma. Glaucoma caused worldwide 6.6% (95% UI:5.9,7.9) of all blindness in 2010 and 2.2% (95% UI:2.0,2.8) of all moderate and severe visual impairment (MSVI). These figures were lower in regions with younger populations (10%). From 1990 to 2010, the number of blind or visually impaired due to glaucoma increased by 0.8 million (95%UI:0.7, 1.1) or 62% and by 2.3 million (95%UI:2.1,3.5) or 83%, respectively. Percentage of global blindness caused by glaucoma increased between 1990 and 2010 from 4.4% (4.0,5.1) to 6.6%. Age-standardized prevalence of glaucoma related blindness and MSVI did not differ markedly between world regions nor between women. Significance: By 2010, one out of 15 blind people was blind due to glaucoma, and one of 45 visually impaired people was visually impaired, highlighting the increasing global burden of glaucoma

Evolutionary and Transmission Dynamics of Reassortant H5N1 Influenza Virus in Indonesia

H5N1 highly pathogenic avian influenza (HPAI) viruses have seriously affected the Asian poultry industry since their recurrence in 2003. The viruses pose a threat of emergence of a global pandemic influenza through point mutation or reassortment leading to a strain that can effectively transmit among humans. In this study, we present phylogenetic evidences for the interlineage reassortment among H5N1 HPAI viruses isolated from humans, cats, and birds in Indonesia, and identify the potential genetic parents of the reassorted genome segments. Parsimony analyses of viral phylogeography suggest that the reassortant viruses may have originated from greater Jakarta and surroundings, and subsequently spread to other regions in the West Java province. In addition, Bayesian methods were used to elucidate the genetic diversity dynamics of the reassortant strain and one of its genetic parents, which revealed a more rapid initial growth of genetic diversity in the reassortant viruses relative to their genetic parent. These results demonstrate that interlineage exchange of genetic information may play a pivotal role in determining viral genetic diversity in a focal population. Moreover, our study also revealed significantly stronger diversifying selection on the M1 and PB2 genes in the lineages preceding and subsequent to the emergence of the reassortant viruses, respectively. We discuss how the corresponding mutations might drive the adaptation and onward transmission of the newly formed reassortant viruses