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. To find the best haploid genome representation, we constructed consensus genomes at the pan-human, superpopulation, and population levels, using variant information from The 1000 Genomes Project Consortium. 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 approximately two to three 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 over using the pan-human consensus, suggesting a limit in the utility of incorporating a more specific genomic variation. Replacing the reference with consensus genomes impacts functional analyses, such as differential expressions of isoforms, genes, and splice junctions

High-fidelity promoter profiling reveals widespread alternative promoter usage and transposon-driven developmental gene expression

Many Eukaryotic genes possess multiple alternative promoters with distinct expression specificities. Therefore, comprehensively annotating promoters and deciphering their individual regulatory dynamics is critical for gene expression profiling applications, and for our understanding of regulatory complexity. We introduce RAMPAGE, a novel promoter activity profiling approach that combines extremely specific 5'-complete cDNA sequencing with an integrated data analysis workflow to address the limitations of current techniques. RAMPAGE features a streamlined protocol for fast and easy generation of highly multiplexed sequencing libraries, offers very high transcription start site specificity, generates accurate and reproducible promoter expression measurements, and yields extensive transcript connectivity information through paired-end cDNA sequencing. We used RAMPAGE in a genome-wide study of promoter activity throughout 36 stages of the life cycle of Drosophila melanogaster, and describe here a comprehensive dataset that represents the first available developmental timecourse of promoter usage. We found that over 40% of developmentally expressed genes have at least 2 promoters, and that alternative promoters generally implement distinct regulatory programs. Transposable elements, long proposed to play a central role in the evolution of their host genomes through their ability to regulate gene expression, contribute at least 1,300 promoters shaping the developmental transcriptome of D. melanogaster. Hundreds of these promoters drive the expression of annotated genes, and transposons often impart their own expression specificity upon the genes they regulate. These observations provide support for the theory that transposons may drive regulatory innovation through the distribution of stereotyped cis-regulatory modules throughout their host genomes

The fractured landscape of RNA-seq alignment: the default in our STARs

Many tools are available for RNA-seq alignment and expression quantification, with comparative value being hard to establish. Benchmarking assessments often highlight methods' good performance, but are focused on either model data or fail to explain variation in performance. This leaves us to ask, what is the most meaningful way to assess different alignment choices? And importantly, where is there room for progress? In this work, we explore the answers to these two questions by performing an exhaustive assessment of the STAR aligner. We assess STAR's performance across a range of alignment parameters using common metrics, and then on biologically focused tasks. We find technical metrics such as fraction mapping or expression profile correlation to be uninformative, capturing properties unlikely to have any role in biological discovery. Surprisingly, we find that changes in alignment parameters within a wide range have little impact on both technical and biological performance. Yet, when performance finally does break, it happens in difficult regions, such as X-Y paralogs and MHC genes. We believe improved reporting by developers will help establish where results are likely to be robust or fragile, providing a better baseline to establish where methodological progress can still occur

Exchange Field Induced Magnetoresistance in Colossal Magnetoresistance Manganites

The effect of an exchange field on electrical transport in thin films of metallic ferromagnetic manganites has been investigated. The exchange field was induced both by direct exchange coupling in a ferromagnet/antiferromagnet multilayer and by indirect exchange interaction in a ferromagnet/paramagnet superlattice. The electrical resistance of the manganite layers was found to be determined by the absolute value of the vector sum of the effective exchange field and the external magnetic field.Comment: 5 pages, 4 figure

Comparison of the transcriptional landscapes between human and mouse tissues

Although the similarities between humans and mice are typically highlighted, morphologically and genetically, there are many differences. To better understand these two species on a molecular level, we performed a comparison of the expression profiles of 15 tissues by deep RNA sequencing and examined the similarities and differences in the transcriptome for both protein-coding and -noncoding transcripts. Although commonalities are evident in the expression of tissue-specific genes between the two species, the expression for many sets of genes was found to be more similar in different tissues within the same species than between species. These findings were further corroborated by associated epigenetic histone mark analyses. We also find that many noncoding transcripts are expressed at a low level and are not detectable at appreciable levels across individuals. Moreover, the majority lack obvious sequence homologs between species, even when we restrict our attention to those which are most highly reproducible across biological replicates. Overall, our results indicate that there is considerable RNA expression diversity between humans and mice, well beyond what was described previously, likely reflecting the fundamental physiological differences between these two organisms

Accuracy assessment of fusion transcript detection via read-mapping and de novo fusion transcript assembly-based methods

BACKGROUND: Accurate fusion transcript detection is essential for comprehensive characterization of cancer transcriptomes. Over the last decade, multiple bioinformatic tools have been developed to predict fusions from RNA-seq, based on either read mapping or de novo fusion transcript assembly. RESULTS: We benchmark 23 different methods including applications we develop, STAR-Fusion and TrinityFusion, leveraging both simulated and real RNA-seq. Overall, STAR-Fusion, Arriba, and STAR-SEQR are the most accurate and fastest for fusion detection on cancer transcriptomes. CONCLUSION: The lower accuracy of de novo assembly-based methods notwithstanding, they are useful for reconstructing fusion isoforms and tumor viruses, both of which are important in cancer research

Oscillatory Exchange Coupling and Positive Magnetoresistance in Epitaxial Oxide Heterostructures

Oscillations in the exchange coupling between ferromagnetic $La_{2/3}Ba_{1/3}MnO_3$ layers with paramagnetic $LaNiO_3$ spacer layer thickness has been observed in epitaxial heterostructures of the two oxides. This behavior is explained within the RKKY model employing an {\it ab initio} calculated band structure of $LaNiO_3$, taking into account strong electron scattering in the spacer. Antiferromagnetically coupled superlattices exhibit a positive current-in-plane magnetoresistance.Comment: 4 pages (RevTeX), 5 figures (EPS

Diabetes causes marked inhibition of mitochondrial metabolism in pancreatic β-cells

Diabetes is a global health problem caused primarily by the inability of pancreatic β-cells to secrete adequate levels of insulin. The molecular mechanisms underlying the progressive failure of β-cells to respond to glucose in type-2 diabetes remain unresolved. Using a combination of transcriptomics and proteomics, we find significant dysregulation of major metabolic pathways in islets of diabetic βV59M mice, a non-obese, eulipidaemic diabetes model. Multiple genes/proteins involved in glycolysis/gluconeogenesis are upregulated, whereas those involved in oxidative phosphorylation are downregulated. In isolated islets, glucose-induced increases in NADH and ATP are impaired and both oxidative and glycolytic glucose metabolism are reduced. INS-1 β-cells cultured chronically at high glucose show similar changes in protein expression and reduced glucose-stimulated oxygen consumption: targeted metabolomics reveals impaired metabolism. These data indicate hyperglycaemia induces metabolic changes in β-cells that markedly reduce mitochondrial metabolism and ATP synthesis. We propose this underlies the progressive failure of β-cells in diabetes.Peer reviewe

Genome-wide polyadenylation site mapping datasets in the rice blast fungus Magnaporthe oryzae

Polyadenylation plays an important role in gene regulation, thus affecting a wide variety of biological processes. In the rice blast fungus Magnaporthe oryzae the cleavage factor I protein Rpb35 is required for pre-mRNA polyadenylation and fungal virulence. Here we present the bioinformatic approach and output data related to a global survey of polyadenylation site usage in M. oryzae wild-type and Delta rbp35 strains under a variety of nutrient conditions, some of which simulate the conditions experienced by the fungus during part of its infection cycle