New methods for next generation sequencing based microRNA expression profiling

<p>Abstract</p> <p>Background</p> <p>MicroRNAs are small non-coding RNA transcripts that regulate post-transcriptional gene expression. The millions of short sequence reads generated by next generation sequencing technologies make this technique explicitly suitable for profiling of known and novel microRNAs. A modification to the small-RNA expression kit (SREK, Ambion) library preparation method for the SOLiD sequencing platform is described to generate microRNA sequencing libraries that are compatible with the Illumina Genome Analyzer.</p> <p>Results</p> <p>High quality sequencing libraries can successfully be prepared from as little as 100 ng small RNA enriched RNA. An easy to use perl-based analysis pipeline called E-miR was developed to handle the sequencing data in several automated steps including data format conversion, 3' adapter removal, genome alignment and annotation to non-coding RNA transcripts. The sample preparation and E-miR pipeline were used to identify 37 cardiac enriched microRNAs in stage 16 chicken embryos. Isomir expression profiles between the heart and embryo were highly correlated for all miRNAs suggesting that tissue or cell specific miRNA modifications do not occur.</p> <p>Conclusions</p> <p>In conclusion, our alternative sample preparation method can successfully be applied to generate high quality miRNA sequencing libraries for the Illumina genome analyzer.</p

Iron accumulation induces oxidative stress, while depressing inflammatory polarization in human iPSC-derived microglia

Iron accumulation in microglia has been observed in Alzheimer’s disease and other neurodegenerative disorders and is thought to contribute to disease progression through various mechanisms, including neuroinflammation. To study this interaction, we treated human induced pluripotent stem cell-derived microglia (iPSC-MG) with iron, in combination with inflammatory stimuli such as interferon gamma (IFN-γ) and amyloid β. Both IFN-γ and iron treatment increased labile iron levels, but only iron treatment led to a consistent increase of ferritin levels, reflecting long-term iron storage. Therefore, in iPSC-MG, ferritin appeared to be regulated by iron revels rather than inflammation. Further investigation showed that while IFN-γ induced pro-inflammatory activation, iron treatment dampened both classic pro- and anti-inflammatory activation on a transcriptomic level. Notably, iron-loaded microglia showed strong upregulation of cellular stress response pathways, the NRF2 pathway, and other oxidative stress pathways. Functionally, iPSC-MG exhibited altered phagocytosis and impaired mitochondrial metabolism following iron treatment. Collectively, these data suggest that in MG, in contrast to current hypotheses, iron treatment does not result in pro-inflammatory activation, but rather dampens it and induces oxidative stress

Determining the quality and complexity of next-generation sequencing data without a reference genome

Genomics, epigenetics, population genetics and bioinformatic

Diagnosis of Fanconi Anemia: Mutation Analysis by Next-Generation Sequencing

Fanconi anemia (FA) is a rare genetic instability syndrome characterized by developmental defects, bone marrow failure, and a high cancer risk. Fifteen genetic subtypes have been distinguished. The majority of patients (≈85%) belong to the subtypes A (≈60%), C (≈15%) or G (≈10%), while a minority (≈15%) is distributed over the remaining 12 subtypes. All subtypes seem to fit within the “classical” FA phenotype, except for D1 and N patients, who have more severe clinical symptoms. Since FA patients need special clinical management, the diagnosis should be firmly established, to exclude conditions with overlapping phenotypes. A valid FA diagnosis requires the detection of pathogenic mutations in a FA gene and/or a positive result from a chromosomal breakage test. Identification of the pathogenic mutations is also important for adequate genetic counselling and to facilitate prenatal or preimplantation genetic diagnosis. Here we describe and validate a comprehensive protocol for the molecular diagnosis of FA, based on massively parallel sequencing. We used this approach to identify BRCA2, FANCD2, FANCI and FANCL mutations in novel unclassified FA patients

Can subtle changes in gene expression be consistently detected with different microarray platforms?

Background: The comparability of gene expression data generated with different microarray platforms is still a matter of concern. Here we address the performance and the overlap in the detection of differentially expressed genes for five different microarray platforms in a challenging biological context where differences in gene expression are few and subtle. Results: Gene expression profiles in the hippocampus of five wild-type and five transgenic δC-doublecortin-like kinase mice were evaluated with five microarray platforms: Applied Biosystems, Affymetrix, Agilent, Illumina, LGTC home-spotted arrays. Using a fixed false discovery rate of 10% we detected surprising differences between the number of differentially expressed genes per platform. Four genes were selected by ABI, 130 by Affymetrix, 3,051 by Agilent, 54 by Illumina, and 13 by LGTC. Two genes were found significantly differentially expressed by all platforms and the four genes identified by the ABI platform were found by at least three other platforms. Quantitative RT-PCR analysis confirmed 20 out of 28 of the genes detected by two or more platforms and 8 out of 15 of the genes detected by Agilent only. We observed improved correlations between platforms when ranking the genes based on the significance level than with a fixed statistical cut-off. We demonstrate significant overlap in the affected gene sets identified by the different platforms, although biological processes were represented by only partially overlapping sets of genes. Aberrances in GABA-ergic signalling in the transgenic mice were consistently found by all platforms. Conclusion: The different microarray platforms give partially complementary views on biological processes affected. Our data indicate that when analyzing samples with only subtle differences in gene expression the use of two different platforms might be more attractive than increasing the number of replicates. Commercial two-color platforms seem to have higher power for finding differentially expressed genes between groups with small differences in expression

Full-length mRNA sequencing uncovers a widespread coupling between transcription initiation and mRNA processing

Background: The multifaceted control of gene expression requires tight coordination of regulatory mechanisms at transcriptional and post-transcriptional level. Here, we studied the interdependence of transcription initiation, splicing and polyadenylation events on single mRNA molecules by full-length mRNA sequencing. Results: In MCF-7 breast cancer cells, we find 2700 genes with interdependent alternative transcription initiation, splicing and polyadenylation events, both in proximal and distant parts of mRNA molecules, including examples of coupling between transcription start sites and polyadenylation sites. The analysis of three human primary tissues (brain, heart and liver) reveals similar patterns of interdependency between transcription initiation and mRNA processing events. We predict thousands of novel open reading frames from full-length mRNA sequences and obtained evidence for their translation by shotgun proteomics. The mapping database rescues 358 previously unassigned peptides and improves the assignment of others. By recognizing sample-specific amino-acid changes and novel splicing patterns, full-length mRNA sequencing improves proteogenomics analysis of MCF-7 cells. Conclusions: Our findings demonstrate that our understanding of transcriptome complexity is far from complete and provides a basis to reveal largely unresolved mechanisms that coordinate transcription initiation and mRNA processing. Electronic supplementary material The online version of this article (10.1186/s13059-018-1418-0) contains supplementary material, which is available to authorized users

Experiences with array-based sequence capture; toward clinical applications

Although sequencing of a human genome gradually becomes an option, zooming in on the region of interest remains attractive and cost saving. We performed array-based sequence capture using 385K Roche NimbleGen, Inc. arrays to zoom in on the protein-coding and immediate intron-flanking sequences of 112 genes, potentially involved in mental retardation and congenital malformation. Captured material was sequenced using Illumina technology. A data analysis pipeline was built that detects sequence variants, positions them in relation to the gene, checks for presence in databases (eg, db single-nucleotide polymorphism (SNP)) and predicts the potential consequences at the level of RNA splicing and protein translation. In the samples analyzed, all known variants were reliably detected, including pathogenic variants from control cases and SNPs derived from array experiments. Although overall coverage varied considerably, it was reproducible per region and facilitated the detection of large deletions and duplications (copy number variations), including a partial deletion in the B3GALTL gene from a patient sample. For ultimate diagnostic application, overall results need to be improved. Future arrays should contain probes from both DNA strands, and to obtain a more even coverage, one could add fewer probes from densely and more probes from sparsely covered regions

Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms

The hippocampal expression profiles of wild-type mice and mice transgenic for δC-doublecortin-like kinase were compared with Solexa/Illumina deep sequencing technology and five different microarray platforms. With Illumina's digital gene expression assay, we obtained ∼2.4 million sequence tags per sample, their abundance spanning four orders of magnitude. Results were highly reproducible, even across laboratories. With a dedicated Bayesian model, we found differential expression of 3179 transcripts with an estimated false-discovery rate of 8.5%. This is a much higher figure than found for microarrays. The overlap in differentially expressed transcripts found with deep sequencing and microarrays was most significant for Affymetrix. The changes in expression observed by deep sequencing were larger than observed by microarrays or quantitative PCR. Relevant processes such as calmodulin-dependent protein kinase activity and vesicle transport along microtubules were found affected by deep sequencing but not by microarrays. While undetectable by microarrays, antisense transcription was found for 51% of all genes and alternative polyadenylation for 47%. We conclude that deep sequencing provides a major advance in robustness, comparability and richness of expression profiling data and is expected to boost collaborative, comparative and integrative genomics studies

Directed adenovirus evolution using engineered mutator viral polymerases

Adenoviruses (Ads) are the most frequently used viruses for oncolytic and gene therapy purposes. Most Ad-based vectors have been generated through rational design. Although this led to significant vector improvements, it is often hampered by an insufficient understanding of Ad’s intricate functions and interactions. Here, to evade this issue, we adopted a novel, mutator Ad polymerase-based, ‘accelerated-evolution’ approach that can serve as general method to generate or optimize adenoviral vectors. First, we site specifically substituted Ad polymerase residues located in either the nucleotide binding pocket or the exonuclease domain. This yielded several polymerase mutants that, while fully supportive of viral replication, increased Ad’s intrinsic mutation rate. Mutator activities of these mutants were revealed by performing deep sequencing on pools of replicated viruses. The strongest identified mutators carried replacements of residues implicated in ssDNA binding at the exonuclease active site. Next, we exploited these mutators to generate the genetic diversity required for directed Ad evolution. Using this new forward genetics approach, we isolated viral mutants with improved cytolytic activity. These mutants revealed a common mutation in a splice acceptor site preceding the gene for the adenovirus death protein (ADP). Accordingly, the isolated viruses showed high and untimely expression of ADP, correlating with a severe deregulation of E3 transcript splicing