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.

Integrating Whole-Genome Sequencing in Clinical Genetics: A Novel Disruptive Structural Rearrangement Identified in the Dystrophin Gene (DMD).

While in most patients the identification of genetic alterations causing dystrophinopathies is a relatively straightforward task, a significant number require genomic and transcriptomic approaches that go beyond a routine diagnostic set-up. In this work, we present a Becker Muscular Dystrophy patient with elevated creatinine kinase levels, progressive muscle weakness, mild intellectual disability and a muscle biopsy showing dystrophic features and irregular dystrophin labelling. Routine molecular techniques (Southern-blot analysis, multiplex PCR, MLPA and genomic DNA sequencing) failed to detect a defect in the DMD gene. Muscle DMD transcript analysis (RT-PCR and cDNA-MLPA) showed the absence of exons 75 to 79, seen to be present at the genomic level. These results prompted the application of low-coverage linked-read whole-genome sequencing (WGS), revealing a possible rearrangement involving DMD intron 74 and a region located upstream of the PRDX4 gene. Breakpoint PCR and Sanger sequencing confirmed the presence of a ~8 Mb genomic inversion. Aberrant DMD transcripts were subsequently identified, some of which contained segments from the region upstream of PRDX4. Besides expanding the mutational spectrum of the disorder, this study reinforces the importance of transcript analysis in the diagnosis of dystrophinopathies and shows how WGS has a legitimate role in clinical laboratory genetics.

Immune stimuli shape the small non-coding transcriptome of extracellular vesicles released by dendritic cells.

The release and uptake of nano-sized extracellular vesicles (EV) is a highly conserved means of intercellular communication. The molecular composition of EV, and thereby their signaling function to target cells, is regulated by cellular activation and differentiation stimuli. EV are regarded as snapshots of cells and are, therefore, in the limelight as biomarkers for disease. Although research on EV-associated RNA has predominantly focused on microRNAs, the transcriptome of EV consists of multiple classes of small non-coding RNAs with potential gene-regulatory functions. It is not known whether environmental cues imposed on cells induce specific changes in a broad range of EV-associated RNA classes. Here, we investigated whether immune-activating or -suppressing stimuli imposed on primary dendritic cells affected the release of various small non-coding RNAs via EV. The small RNA transcriptomes of highly pure EV populations free from ribonucleoprotein particles were analyzed by RNA sequencing and RT-qPCR. Immune stimulus-specific changes were found in the miRNA, snoRNA, and Y-RNA content of EV from dendritic cells, whereas tRNA and snRNA levels were much less affected. Only part of the changes in EV-RNA content reflected changes in cellular RNA, which urges caution in interpreting EV as snapshots of cells. By comprehensive analysis of RNA obtained from highly purified EV, we demonstrate that multiple RNA classes contribute to genetic messages conveyed via EV. The identification of multiple RNA classes that display cell stimulation-dependent association with EV is the prelude to unraveling the function and biomarker potential of these EV-RNAs.

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.

Huntington's disease biomarker progression profile identified by transcriptome sequencing in peripheral blood.

With several therapeutic approaches in development for Huntington's disease, there is a need for easily accessible biomarkers to monitor disease progression and therapy response. We performed next-generation sequencing-based transcriptome analysis of total RNA from peripheral blood of 91 mutation carriers (27 presymptomatic and, 64 symptomatic) and 33 controls. Transcriptome analysis by DeepSAGE identified 167 genes significantly associated with clinical total motor score in Huntington's disease patients. Relative to previous studies, this yielded novel genes and confirmed previously identified genes, such as H2AFY, an overlap in results that has proven difficult in the past. Pathway analysis showed enrichment of genes of the immune system and target genes of miRNAs, which are downregulated in Huntington's disease models. Using a highly parallelized microfluidics array chip (Fluidigm), we validated 12 of the top 20 significant genes in our discovery cohort and 7 in a second independent cohort. The five genes (PROK2, ZNF238, AQP9, CYSTM1 and ANXA3) that were validated independently in both cohorts present a candidate biomarker panel for stage determination and therapeutic readout in Huntington's disease. Finally we suggest a first empiric formula predicting total motor score from the expression levels of our biomarker panel. Our data support the view that peripheral blood is a useful source to identify biomarkers for Huntington's disease and monitor disease progression in future clinical trials.

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.

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.

High-resolution whole-genome sequencing reveals that specific chromatin domains from most human chromosomes associate with nucleoli.

The nuclear space is mostly occupied by chromosome territories and nuclear bodies. Although this organization of chromosomes affects gene function, relatively little is known about the role of nuclear bodies in the organization of chromosomal regions. The nucleolus is the best-studied subnuclear structure and forms around the rRNA repeat gene clusters on the acrocentric chromosomes. In addition to rDNA, other chromatin sequences also surround the nucleolar surface and may even loop into the nucleolus. These additional nucleolar-associated domains (NADs) have not been well characterized. We present here a whole-genome, high-resolution analysis of chromatin endogenously associated with nucleoli. We have used a combination of three complementary approaches, namely fluorescence comparative genome hybridization, high-throughput deep DNA sequencing and photoactivation combined with time-lapse fluorescence microscopy. The data show that specific sequences from most human chromosomes, in addition to the rDNA repeat units, associate with nucleoli in a reproducible and heritable manner. NADs have in common a high density of AT-rich sequence elements, low gene density and a statistically significant enrichment in transcriptionally repressed genes. Unexpectedly, both the direct DNA sequencing and fluorescence photoactivation data show that certain chromatin loci can specifically associate with either the nucleolus, or the nuclear envelope.

Terminal osseous dysplasia is caused by a single recurrent mutation in the FLNA gene.

Terminal osseous dysplasia (TOD) is an X-linked dominant male-lethal disease characterized by skeletal dysplasia of the limbs, pigmentary defects of the skin, and recurrent digital fibroma with onset in female infancy. After performing X-exome capture and sequencing, we identified a mutation at the last nucleotide of exon 31 of the FLNA gene as the most likely cause of the disease. The variant c.5217G>A was found in six unrelated cases (three families and three sporadic cases) and was not found in 400 control X chromosomes, pilot data from the 1000 Genomes Project, or the FLNA gene variant database. In the families, the variant segregated with the disease, and it was transmitted four times from a mildly affected mother to a more seriously affected daughter. We show that, because of nonrandom X chromosome inactivation, the mutant allele was not expressed in patient fibroblasts. RNA expression of the mutant allele was detected only in cultured fibroma cells obtained from 15-year-old surgically removed material. The variant activates a cryptic splice site, removing the last 48 nucleotides from exon 31. At the protein level, this results in a loss of 16 amino acids (p.Val1724_Thr1739del), predicted to remove a sequence at the surface of filamin repeat 15. Our data show that TOD is caused by this single recurrent mutation in the FLNA gene.

Genome-wide assessment of differential roles for p300 and CBP in transcription regulation.

Despite high levels of homology, transcription coactivators p300 and CREB binding protein (CBP) are both indispensable during embryogenesis. They are largely known to regulate the same genes. To identify genes preferentially regulated by p300 or CBP, we performed an extensive genome-wide survey using the ChIP-seq on cell-cycle synchronized cells. We found that 57% of the tags were within genes or proximal promoters, with an overall preference for binding to transcription start and end sites. The heterogeneous binding patterns possibly reflect the divergent roles of CBP and p300 in transcriptional regulation. Most of the 16 103 genes were bound by both CBP and p300. However, after stimulation 89 and 1944 genes were preferentially bound by CBP or p300, respectively. Target genes were found to be primarily involved in the regulation of metabolic and developmental processes, and transcription, with CBP showing a stronger preference than p300 for genes active in negative regulation of transcription. Analysis of transcription factor binding sites suggest that CBP and p300 have many partners in common, but AP-1 and Serum Response Factor (SRF) appear to be more prominent in CBP-specific sequences, whereas AP-2 and SP1 are enriched in p300-specific targets. Taken together, our findings further elucidate the distinct roles of coactivators p300 and CBP in transcriptional regulation.

Deep sequencing to reveal new variants in pooled DNA samples.

We evaluated massive parallel sequencing and long-range PCR (LRP) for rare variant detection and allele frequency estimation in pooled DNA samples. Exons 2 to 16 of the MUTYH gene were analyzed in breast cancer patients with Illumina's (Solexa) technology. From a pool of 287 genomic DNA samples we generated a single LRP product, while the same LRP was performed on 88 individual samples and the resulting products then pooled. Concentrations of constituent samples were measured with fluorimetry for genomic DNA and high-resolution melting curve analysis (HR-MCA) for LRP products. Illumina sequencing results were compared to Sanger sequencing data of individual samples. Correlation between allele frequencies detected by both methods was poor in the first pool, presumably because the genomic samples amplified unequally in the LRP, due to DNA quality variability. In contrast, allele frequencies correlated well in the second pool, in which all expected alleles at a frequency of 1% and higher were reliably detected, plus the majority of singletons (0.6% allele frequency). We describe custom bioinformatics and statistics to optimize detection of rare variants and to estimate required sequencing depth. Our results provide directions for designing high-throughput analyses of candidate genes.

Large deletions in the polycystic kidney disease 1 (PKD1) gene.

Since identification of the genes mutated in patients with Autosomal Dominant Polycystic Kidney Disease, PKD1 and PKD2, a large number of different germ line mutations in both genes have been found by conventional PCR-based mutation detection methods. Nevertheless, in approximately 40% of the PKD1 families the disease-causing mutation remains to be elucidated. Complex germ line rearrangements are often not detectable by these standard diagnostic techniques. To detect large deletions in the PKD1 gene we performed Field Inversion Gel Electrophoresis (FIGE) followed by Southern blot analysis with probes selected in the unique and in the reiterated region of this gene. Our analysis revealed 4 deletions in 125 patients, indicating that large deletions in PKD1 are rare. Likely, patients with a deletion that also affects the neighbouring Tuberous Sclerosis Complex 2 (TSC2) gene will be diagnosed as patients with tuberous sclerosis. It was speculated that the exceptional polypyrimidine tract located in intron 21 and the small tract in intron 22, might play a role in the pathogenesis of ADPKD. Since this region is extremely difficult to amplify by PCR, we analysed the 5.8 kb BamHI fragment that contains the polypyrimidine tracts. We did not observe a disease-linked alteration although we detected two different rare variants either in PKD1 or in one of its homologues.

Loss of CBP acetyltransferase activity by PHD finger mutations in Rubinstein-Taybi syndrome.

Disruption of one copy of the human CREB binding protein (CBP or CREBBP) gene leads to the Rubinstein-Taybi syndrome (RTS), a developmental disorder characterized by retarded growth and mental functions, broad thumbs, broad big toes and typical facial abnormalities. The CREB binding protein (CBP) is an essential transcriptional coactivator for many different transcription factors. CBP has the intrinsic ability to acetylate histones and other proteins, which is regarded as an important step in transcription activation. In vitro studies have shown that this enzymatic activity critically depends on the integrity of a plant homeodomain (PHD)-type zinc finger in the HAT domain of CBP. We therefore investigated whether PHD finger mutations are present in RTS patients. Mutational analysis of 39 patients revealed eight novel heterozygous mutations in the HAT domain of CBP, one of which alters a conserved PHD finger amino acid (E1278K), while a second mutation deletes exon 22, which encodes the central region of the PHD finger. Functional analysis of these RTS-associated PHD finger mutants showed that they lacked in vitro acetyltransferase activity towards histones and CBP itself and displayed reduced coactivator function for the transcription factor CREB. Importantly, in EBV-transformed lymphoblastoid cells from the exon 22 deletion patient we found approximately 50% less endogenous CBP HAT activity. These findings therefore underscore the functional importance of the PHD finger in vivo and imply that reduction of CBP HAT activity, as exemplified here by disruption of the PHD finger, is sufficient to cause RTS.