Is it time to change the reference genome?

The use of the human reference genome has shaped methods and data across modern genomics. This has offered many benefits while creating a few constraints. In the following opinion, we outline the history, properties, and pitfalls of the current human reference genome. In a few illustrative analyses, we focus on its use for variant-calling, highlighting its nearness to a 'type specimen'. We suggest that switching to a consensus reference would offer important advantages over the continued use of the current reference with few disadvantages

Pan-human consensus genome significantly improves the accuracy of RNA-seq analyses

The Human Reference Genome serves as the foundation for modern genomic analyses. However, in its present form, it does not adequately represent the vast genetic diversity of the human population. In this study, we explored the consensus genome as a potential successor of the current reference genome and assessed its effect on the accuracy of RNA-seq read alignment. In order to find the best haploid genome representation, we constructed consensus genomes at the pan-human, super-population, and population levels, utilizing variant information from the 1000 Genomes Project. Using personal haploid genomes as the ground truth, we compared mapping errors for real RNA-seq reads aligned to the consensus genomes versus the reference genome. For reads overlapping homozygous variants, we found that the mapping error decreased by a factor of ~2-3 when the reference was replaced with the pan-human consensus genome. We also found that using more population-specific consensuses resulted in little to no increase overusing the pan-human consensus, suggesting a limit in the utility of incorporating more specific genomic variation. Replacing reference with consensus genomes impacts functional analyses, such as differential expressions of isoforms, genes, and splice junctions

Improved annotation of 3' untranslated regions and complex loci by combination of strand-specific direct RNA sequencing, RNA-seq and ESTs

The reference annotations made for a genome sequence provide the framework for all subsequent analyses of the genome. Correct annotation is particularly important when interpreting the results of RNA-seq experiments where short sequence reads are mapped against the genome and assigned to genes according to the annotation. Inconsistencies in annotations between the reference and the experimental system can lead to incorrect interpretation of the effect on RNA expression of an experimental treatment or mutation in the system under study. Until recently, the genome-wide annotation of 3-prime untranslated regions received less attention than coding regions and the delineation of intron/exon boundaries. In this paper, data produced for samples in Human, Chicken and A. thaliana by the novel single-molecule, strand-specific, Direct RNA Sequencing technology from Helicos Biosciences which locates 3-prime polyadenylation sites to within +/- 2 nt, were combined with archival EST and RNA-Seq data. Nine examples are illustrated where this combination of data allowed: (1) gene and 3-prime UTR re-annotation (including extension of one 3-prime UTR by 5.9 kb); (2) disentangling of gene expression in complex regions; (3) clearer interpretation of small RNA expression and (4) identification of novel genes. While the specific examples displayed here may become obsolete as genome sequences and their annotations are refined, the principles laid out in this paper will be of general use both to those annotating genomes and those seeking to interpret existing publically available annotations in the context of their own experimental dataComment: 44 pages, 9 figure

A novel class of microRNA-recognition elements that function only within open reading frames.

MicroRNAs (miRNAs) are well known to target 3' untranslated regions (3' UTRs) in mRNAs, thereby silencing gene expression at the post-transcriptional level. Multiple reports have also indicated the ability of miRNAs to target protein-coding sequences (CDS); however, miRNAs have been generally believed to function through similar mechanisms regardless of the locations of their sites of action. Here, we report a class of miRNA-recognition elements (MREs) that function exclusively in CDS regions. Through functional and mechanistic characterization of these 'unusual' MREs, we demonstrate that CDS-targeted miRNAs require extensive base-pairing at the 3' side rather than the 5' seed; cause gene silencing in an Argonaute-dependent but GW182-independent manner; and repress translation by inducing transient ribosome stalling instead of mRNA destabilization. These findings reveal distinct mechanisms and functional consequences of miRNAs that target CDS versus the 3' UTR and suggest that CDS-targeted miRNAs may use a translational quality-control-related mechanism to regulate translation in mammalian cells

MRP3 is a sex determining gene in the diatom Pseudo-nitzschia multistriata

A broad diversity of sex-determining systems has evolved in eukaryotes. However, information on the mechanisms of sex determination for unicellular microalgae is limited, including for diatoms, key-players of ocean food webs. Here we report the identification of a mating type (MT) determining gene for the diatom Pseudo-nitzschia multistriata. By comparing the expression profile of the two MTs, we find five MT-biased genes, of which one, MRP3, is expressed exclusively in MT+ strains in a monoallelic manner. A short tandem repeat of specific length in the region upstream of MRP3 is consistently present in MT+ and absent in MT- strains. MRP3 overexpression in an MT- strain induces sex reversal: the transgenic MT- can mate with another MT- strain and displays altered regulation of the other MT-biased genes, indicating that they lie downstream. Our data show that a relatively simple genetic program is involved in defining the MT in P. multistriata

GI-530159, a novel, selective, mechanosensitive two‐pore‐domain potassium (K2P_{2P}) channel opener, reduces rat dorsal root ganglion neuron excitability

Background and Purpose: TREK two pore domain potassium channels play a critical role in regulating the excitability of somatosensory nociceptive neurons and are important mediators of pain perception. An understanding of the roles of TREK channels in pain perception and, indeed, in other pathophysiological conditions, has been severely hampered by the lack of potent and/or selective activators and inhibitors. In this study we describe a new, selective opener of TREK channels, GI-530159 (4,4?-(Hexafluoroisopropylidene)bis(p-phenyleneoxy)dianiline). Experimental Approach: The effect of GI-530159 on TREK channels was demonstrated using 86Rb efflux assays, whole-cell and single channel patch clamp recordings from recombinant TREK channels. TREK1, TREK2 and TRAAK expression was determined using transcriptome analysis from single dorsal root ganglion (DRG) cells. Current-clamp recordings from cultured rat DRG neurons were used to measure the effect of GI-530159 on neuronal excitability. Key Results: For recombinant human TREK1 channels GI-530159 had an EC50 of 0.8?M from Rb efflux experiments and 0.9?M from electrophysiological recordings. The compound activated TREK2 channels but it had no detectable action on TRAAK channels nor any significant effect on other K channels tested. Current clamp recordings from cultured rat DRG neurones showed that application of GI-530159 at 1?M resulted in a significant reduction in firing frequency and a small hyperpolarisation of resting membrane potential. Conclusions and Implications: This study provides pharmacological evidence for the presence of mechano-sensitive TREK K2P channels in sensory neurones and suggests that development of selective K2P channel openers like GI-530159 could aid to the development of novel analgesic agents