Search

Age-associated chronic inflammation (inflammageing) is a central hallmark of ageing1, but its influence on specific cells remains largely unknown. Fibroblasts are present in most tissues and contribute to wound healing2,3. They are also the most widely used cell type for reprogramming to induced pluripotent stem (iPS) cells, a process that has implications for regenerative medicine and rejuvenation strategies4. Here we show that fibroblast cultures from old mice secrete inflammatory cytokines and exhibit increased variability in the efficiency of iPS cell reprogramming between mice. Variability between individuals is emerging as a feature of old age5–8, but the underlying mechanisms remain unknown. To identify drivers of this variability, we performed multi-omics profiling of fibroblast cultures from young and old mice that have different reprogramming efficiencies. This approach revealed that fibroblast cultures from old mice contain ‘activated fibroblasts’ that secrete inflammatory cytokines, and that the proportion of activated fibroblasts in a culture correlates with the reprogramming efficiency of that culture. Experiments in which conditioned medium was swapped between cultures showed that extrinsic factors secreted by activated fibroblasts underlie part of the variability between mice in reprogramming efficiency, and we have identified inflammatory cytokines, including TNF, as key contributors. Notably, old mice also exhibited variability in wound healing rate in vivo. Single-cell RNA-sequencing analysis identified distinct subpopulations of fibroblasts with different cytokine expression and signalling in the wounds of old mice with slow versus fast healing rates. Hence, a shift in fibroblast composition, and the ratio of inflammatory cytokines that they secrete, may drive the variability between mice in reprogramming in vitro and influence wound healing rate in vivo. This variability may reflect distinct stochastic ageing trajectories between individuals, and could help in developing personalized strategies to improve iPS cell generation and wound healing in elderly individuals. Fibroblasts from old mice are heterogeneous, which affects the ability of these fibroblasts to reprogram into induced pluripotent stem cells in vitro and influences wound healing rate in vivo.

Rationale: Regulatory DNA elements in the human genome play important roles in determining the transcriptional abundance and spatiotemporal gene expression during embryonic heart development and somatic cell reprogramming. It is not well known how chromatin marks in regulatory DNA elements are modulated to establish cell type–specific gene expression in the human heart. Objective: We aimed to decipher the cell type–specific epigenetic signatures in regulatory DNA elements and how they modulate heart-specific gene expression. Methods and Results: We profiled genome-wide transcriptional activity and a variety of epigenetic marks in the regulatory DNA elements using massive RNA-seq (n=12) and ChIP-seq (chromatin immunoprecipitation combined with high-throughput sequencing; n=84) in human endothelial cells (CD31 + CD144 + ), cardiac progenitor cells (Sca-1 + ), fibroblasts (DDR2 + ), and their respective induced pluripotent stem cells. We uncovered 2 classes of regulatory DNA elements: class I was identified with ubiquitous enhancer (H3K4me1) and promoter (H3K4me3) marks in all cell types, whereas class II was enriched with H3K4me1 and H3K4me3 in a cell type–specific manner. Both class I and class II regulatory elements exhibited stimulatory roles in nearby gene expression in a given cell type. However, class I promoters displayed more dominant regulatory effects on transcriptional abundance regardless of distal enhancers. Transcription factor network analysis indicated that human induced pluripotent stem cells and somatic cells from the heart selected their preferential regulatory elements to maintain cell type–specific gene expression. In addition, we validated the function of these enhancer elements in transgenic mouse embryos and human cells and identified a few enhancers that could possibly regulate the cardiac-specific gene expression. Conclusions: Given that a large number of genetic variants associated with human diseases are located in regulatory DNA elements, our study provides valuable resources for deciphering the epigenetic modulation of regulatory DNA elements that fine-tune spatiotemporal gene expression in human cardiac development and diseases.

Essential tremor (ET) is the most common adult-onset movement disorder. In the present study, we performed whole exome sequencing of a large ET-affected family (10 affected and 6 un-affected family members) and identified a TUB p.V431I variant (rs75594955) segregating in a manner consistent with autosomal-dominant inheritance. Subsequent targeted re-sequencing of TUB in 820 unrelated individuals with sporadic ET and 630 controls revealed significant enrichment of rare nonsynonymous TUB variants (e.g. rs75594955: p.V431I, rs1241709665: p.Ile20Phe, rs55648406: p.Arg49Gln) in the ET cohort (SKAT-O test p-value = 6.20e-08). TUB encodes a transcription factor predominantly expressed in neuronal cells and has been previously implicated in obesity. ChIP-seq analyses of the TUB transcription factor across different regions of the mouse brain revealed that TUB regulates the pathways responsible for neurotransmitter production as well thyroid hormone signaling. Together, these results support the association of rare variants in TUB with ET., Author summary Essential tremor (ET) is the most common adult-onset movement disorder and in most affected families it appears to be inherited in an autosomal dominant pattern. The causes of essential tremor are unknown. Although many genetic studies in affected families and sporadic cases of ET have shown that genes may play a role, it has proven quite challenging to identify the specific genetic variants involved. Here, we use state-of-the-art technologies to identify the role of genetic variants on ET through exome sequencing of a large affected ET family and subsequent validation in a large population of cases and controls. We show that rare nonsynonymous variants of the TUB gene are significantly enriched in ET cases versus healthy controls. Further studies of biological pathways regulated by TUB in the mouse brain reveal key pathways related to ET. Our work expands our knowledge of the genetic basis of ET.

Left ventricular non-compaction (LVNC) is the third most prevalent cardiomyopathy in children and its pathogenesis has been associated with the developmental defect of the embryonic myocardium. We show that patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) generated from LVNC patients carrying a mutation in the cardiac transcription factor TBX20 recapitulate a key aspect of the pathological phenotype at the single-cell level and this was associated with perturbed transforming growth factor beta (TGF-β) signalling. LVNC iPSC-CMs have decreased proliferative capacity due to abnormal activation of TGF-β signalling. TBX20 regulates the expression of TGF-β signalling modifiers including one known to be a genetic cause of LVNC, PRDM16, and genome editing of PRDM16 caused proliferation defects in iPSC-CMs. Inhibition of TGF-β signalling and genome correction of the TBX20 mutation were sufficient to reverse the disease phenotype. Our study demonstrates that iPSC-CMs are a useful tool for the exploration of pathological mechanisms underlying poorly understood cardiomyopathies including LVNC.

Age-associated chronic inflammation (inflammaging) has emerged as a central hallmark of aging1-3, but its impact on specific cells is still largely unknown. Fibroblasts are present in all tissues and contribute to wound healing4-6. They are also the cell type that is mostly used for induced pluripotent stem cell (iPSC) reprogramming7 – a process that has implications for regenerative medicine and rejuvenation strategies8-17. Here we show that primary fibroblasts from old mice secrete inflammatory cytokines and that there is an increased variability in reprogramming efficiency between fibroblast cultures from old individuals. Individual-to-individual variability is emerging as a key feature of old age18-21, which could reflect distinct aging trajectories, but the underlying causes remain unknown. To identify drivers of this variability, we perform a multi-omic assessment of young and old fibroblast cultures with different reprogramming efficiency. This approach, coupled with single cell transcriptomics, reveals that old fibroblast cultures are heterogeneous and show a greater proportion of ‘activated fibroblasts’ that secrete inflammatory cytokines, which correlates with reprogramming efficiency. We experimentally validate that activated fibroblasts express inflammatory cytokines in vivo and that their presence is linked to enhanced reprogramming efficiency in culture. Conditioned-media swapping experiments show that extrinsic factors secreted by activated fibroblasts are more critical than intrinsic factors for the individual-to-individual variability in reprogramming efficiency, and we identify TNFα as a key inflammatory cytokine underlying this variability. Interestingly, old mice also exhibit variability in wound healing efficiency in vivo and old wounds show an increased subpopulation of activated fibroblasts with a unique TNFα signature. Our study shows that a switch in fibroblast composition, and the ratio of inflammatory cytokines they secrete, drives variability in reprogramming in vitro and may influence wound healing in vivo. These findings could help identify personalized strategies to improve iPSC generation and wound healing in older individuals.

Identifying bacterial strains in metagenome and microbiome samples using computational analyses of short-read sequences remains a difficult problem. Here, we present an analysis of a human gut microbiome using TruSeq synthetic long reads combined with computational tools for metagenomic long-read assembly, variant calling and haplotyping (Nanoscope and Lens). Our analysis identifies 178 bacterial species, of which 51 were not found using shotgun reads alone. We recover bacterial contigs that comprise multiple operons, including 22 contigs of >1 Mbp. Furthermore, we observe extensive intraspecies variation within microbial strains in the form of haplotypes that span up to hundreds of Kbp. Incorporation of synthetic long-read sequencing technology with standard short-read approaches enables more precise and comprehensive analyses of metagenomic samples.

Understanding transcriptome complexity is crucial for understanding human biology and disease. Technologies such as Synthetic long-read RNA sequencing (SLR-RNA-seq) delivered 5 million isoforms and allowed assessing splicing coordination. Pacific Biosciences and Oxford Nanopore increase throughput also but require high input amounts or amplification. Our new droplet-based method, sparse isoform sequencing (spISO-seq), sequences 100k–200k partitions of 10–200 molecules at a time, enabling analysis of 10–100 million RNA molecules. SpISO-seq requires less than 1 ng of input cDNA, limiting or removing the need for prior amplification with its associated biases. Adjusting the number of reads devoted to each molecule reduces sequencing lanes and cost, with little loss in detection power. The increased number of molecules expands our understanding of isoform complexity. In addition to confirming our previously published cases of splicing coordination (e.g., BIN1), the greater depth reveals many new cases, such as MAPT. Coordination of internal exons is found to be extensive among protein coding genes: 23.5%–59.3% (95% confidence interval) of highly expressed genes with distant alternative exons exhibit coordination, showcasing the need for long-read transcriptomics. However, coordination is less frequent for noncoding sequences, suggesting a larger role of splicing coordination in shaping proteins. Groups of genes with coordination are involved in protein–protein interactions with each other, raising the possibility that coordination facilitates complex formation and/or function. We also find new splicing coordination types, involving initial and terminal exons. Our results provide a more comprehensive understanding of the human transcriptome and a general, cost-effective method to analyze it.

Alopecia with mental retardation syndrome (APMR) is a very rare autosomal recessive condition that is associated with total or partial absence of hair from the scalp and other parts of the body as well as variable intellectual disability. Here we present whole-exome sequencing results of a large consanguineous family segregating APMR syndrome with seven affected family members. Our study revealed a novel predicted pathogenic, homozygous missense mutation in the AHSG (OMIM 138680) gene (AHSG: NM_001622:exon7:c.950G>A:p.Arg317His). The variant is predicted to affect a region of the protein required for protein processing and disrupts a phosphorylation motif. In addition, the altered protein migrates with an aberrant size relative to healthy individuals. Consistent with the phenotype, AHSG maps within APMR linkage region 1 (APMR 1) as reported before, and falls within runs of homozygosity (ROH). Previous families with APMR syndrome have been studied through linkage analyses and the linkage resolution did not allow pointing out to a single gene candidate. Our study is the first report to identify a homozygous missense mutation for APMR syndrome through whole-exome sequencing.

Human induced pluripotent stem cells (iPSCs) can be derived from various types of somatic cells by transient overexpression of 4 Yamanaka factors (OCT4, SOX2, C-MYC, and KLF4). Patient-specific iPSC derivatives (e.g., neuronal, cardiac, hepatic, muscular, and endothelial cells [ECs]) hold great promise in drug discovery and regenerative medicine. In this study, we aimed to evaluate whether the cellular origin can affect the differentiation, in vivo behavior, and single-cell gene expression signatures of human iPSC–derived ECs. We derived human iPSCs from 3 types of somatic cells of the same individuals: fibroblasts (FB-iPSCs), ECs (EC-iPSCs), and cardiac progenitor cells (CPC-iPSCs). We then differentiated them into ECs by sequential administration of Activin, BMP4, bFGF, and VEGF. EC-iPSCs at early passage (10 < P < 20) showed higher EC differentiation propensity and gene expression of EC-specific markers (PECAM1 and NOS3) than FB-iPSCs and CPC-iPSCs. In vivo transplanted EC-iPSC–ECs were recovered with a higher percentage of CD31+ population and expressed higher EC-specific gene expression markers (PECAM1, KDR, and ICAM) as revealed by microfluidic single-cell quantitative PCR (qPCR). In vitro EC-iPSC–ECs maintained a higher CD31+ population than FB-iPSC–ECs and CPC-iPSC–ECs with long-term culturing and passaging. These results indicate that cellular origin may influence lineage differentiation propensity of human iPSCs; hence, the somatic memory carried by early passage iPSCs should be carefully considered before clinical translation.

© 2016 The Author(s). Insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) expression correlates with malignancy, but its role(s) in pathogenesis remains enigmatic. We interrogated the IGF2BP3-RNA interaction network in pancreatic ductal adenocarcinoma (PDAC) cells. Using a combination of genome-wide approaches, we have identified 164 direct mRNA targets of IGF2BP3. These transcripts encode proteins enriched for functions such as cell migration, proliferation, and adhesion. Loss of IGF2BP3 reduced PDAC cell invasiveness and remodeled focal adhesion junctions. Individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) revealed significant overlap of IGF2BP3 and microRNA (miRNA) binding sites. IGF2BP3 promotes association of the RNA-induced silencing complex (RISC) with specific transcripts. Our results show that IGF2BP3 influences a malignancy-associated RNA regulon by modulating miRNA-mRNA interactions.

Age-associated chronic inflammation (inflammageing) is a central hallmark of ageing

Left ventricular non-compaction (LVNC) is the third most prevalent cardiomyopathy in children and has a unique phenotype with characteristically extensive hypertrabeculation of the left ventricle, similar to the embryonic left ventricle, suggesting a developmental defect of the embryonic myocardium. However, studying this disease has been challenging due to the lack of an animal model that can faithfully recapitulate the clinical phenotype of LVNC. To address this, we showed that patient-specific induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) generated from a family with LVNC history recapitulated a developmental defect consistent with the LVNC phenotype at the single-cell level. We then utilized hiPSC-CMs to show that increased transforming growth factor beta (TGFβ) signaling is one of the central mechanisms underlying the pathogenesis of LVNC. LVNC hiPSC-CMs demonstrated decreased proliferative capacity due to abnormal activation of TGFβ signaling (Figs A-B). Exome sequencing demonstrated a mutation in TBX20, which regulates TGFβ signaling through upregulation of ITGAV, contributing to the LVNC phenotype. Inhibition of abnormal TGFβ signaling or genetic correction of the TBX20 mutation (Figs C-D) using TALEN reversed the proliferation defects seen in LVNC hiPSC-CMs. Our results demonstrate that hiPSC-CMs are a useful tool for the exploration of novel mechanisms underlying poorly understood cardiomyopathies such as LVNC. Here we provide the first evidence of activation of TGFβ signaling as playing a role in the pathogenesis of LVNC.

Background The genetic diversity of loci and mutations underlying hereditary hearing loss is an active area of investigation. To identify loci associated with predominantly non-syndromic sensorineural hearing loss, we performed exome sequencing of families and of single probands, as well as copy number variation (CNV) mapping in a case–control cohort. Results Analysis of three distinct families revealed several candidate loci in two families and a single strong candidate gene, MYH7B, for hearing loss in one family. MYH7B encodes a Type II myosin, consistent with a role for cytoskeletal proteins in hearing. High-resolution genome-wide CNV analysis of 150 cases and 157 controls revealed deletions in genes known to be involved in hearing (e.g. GJB6, OTOA, and STRC, encoding connexin 30, otoancorin, and stereocilin, respectively), supporting CNV contributions to hearing loss phenotypes. Additionally, a novel region on chromosome 16 containing part of the PDXDC1 gene was found to be frequently deleted in hearing loss patients (OR = 3.91, 95% CI: 1.62-9.40, p = 1.45 × 10-7). Conclusions We conclude that many known as well as novel loci and distinct types of mutations not typically tested in clinical settings can contribute to the etiology of hearing loss. Our study also demonstrates the challenges of exome sequencing and genome-wide CNV mapping for direct clinical application, and illustrates the need for functional and clinical follow-up as well as curated open-access databases. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-1155) contains supplementary material, which is available to authorized users.

Alternative splicing shapes mammalian transcriptomes, with many RNA molecules undergoing multiple distant alternative splicing events. Comprehensive transcriptome analysis, including analysis of exon co-association in the same molecule, requires deep, long-read sequencing. Here we introduce an RNA sequencing method, synthetic long-read RNA sequencing (SLR-RNA-seq), in which small pools (≤1,000 molecules/pool, ≤1 molecule/gene for most genes) of full-length cDNAs are amplified, fragmented and short-read-sequenced. We demonstrate that these RNA sequences reconstructed from the short reads from each of the pools are mostly close to full length and contain few insertion and deletion errors. We report many previously undescribed isoforms (human brain: ∼13,800 affected genes, 14.5% of molecules; mouse brain ∼8,600 genes, 18% of molecules) and up to 165 human distant molecularly associated exon pairs (dMAPs) and distant molecularly and mutually exclusive pairs (dMEPs). Of 16 associated pairs detected in the mouse brain, 9 are conserved in human. Our results indicate conserved mechanisms that can produce distant but phased features on transcript and proteome isoforms.

Increasing evidence suggests that interactions between regulatory genomic elements play an important role in regulating gene expression. We generated a genome-wide interaction map of regulatory elements in human cells (ENCODE tier 1 cells, K562, GM12878) using Chromatin Interaction Analysis by Paired-End Tag sequencing (ChIA-PET) experiments targeting six broadly distributed factors. Bound regions covered 80% of DNase I hypersensitive sites including 99.7% of TSS and 98% of enhancers. Correlating this map with ChIP-seq and RNA-seq data sets revealed cohesin, CTCF, and ZNF143 as key components of three-dimensional chromatin structure and revealed how the distal chromatin state affects gene transcription. Comparison of interactions between cell types revealed that enhancer–promoter interactions were highly cell-type-specific. Construction and comparison of distal and proximal regulatory networks revealed stark differences in structure and biological function. Proximal binding events are enriched at genes with housekeeping functions, while distal binding events interact with genes involved in dynamic biological processes including response to stimulus. This study reveals new mechanistic and functional insights into regulatory region organization in the nucleus.