Dynamic interplay between RPL3- and RPL3L-containing ribosomes modulates mitochondrial activity in the mammalian heart

The existence of naturally occurring ribosome heterogeneity is now a well-acknowledged phenomenon. However, whether this heterogeneity leads to functionally diverse 'specialized ribosomes' is still a controversial topic. Here, we explore the biological function of RPL3L (uL3L), a ribosomal protein (RP) paralogue of RPL3 (uL3) that is exclusively expressed in skeletal muscle and heart tissues, by generating a viable homozygous Rpl3l knockout mouse strain. We identify a rescue mechanism in which, upon RPL3L depletion, RPL3 becomes up-regulated, yielding RPL3-containing ribosomes instead of RPL3L-containing ribosomes that are typically found in cardiomyocytes. Using both ribosome profiling (Ribo-seq) and a novel orthogonal approach consisting of ribosome pulldown coupled to nanopore sequencing (Nano-TRAP), we find that RPL3L modulates neither translational efficiency nor ribosome affinity towards a specific subset of transcripts. In contrast, we show that depletion of RPL3L leads to increased ribosome-mitochondria interactions in cardiomyocytes, which is accompanied by a significant increase in ATP levels, potentially as a result of fine-tuning of mitochondrial activity. Our results demonstrate that the existence of tissue-specific RP paralogues does not necessarily lead to enhanced translation of specific transcripts or modulation of translational output. Instead, we reveal a complex cellular scenario in which RPL3L modulates the expression of RPL3, which in turn affects ribosomal subcellular localization and, ultimately, mitochondrial activity.

Identification of novel genes whose expression in adipose tissue affects body fat mass and distribution: an RNA-Seq and Mendelian Randomization study

Many studies have shown that abdominal adiposity is more strongly related to health risks than peripheral adiposity. However, the underlying pathways are still poorly understood. In this cross-sectional study using data from RNA-sequencing experiments and whole-body MRI scans of 200 participants in the EPIC-Potsdam cohort, our aim was to identify novel genes whose gene expression in subcutaneous adipose tissue has an effect on body fat mass (BFM) and body fat distribution (BFD). The analysis identified 625 genes associated with adiposity, of which 531 encode a known protein and 487 are novel candidate genes for obesity. Enrichment analyses indicated that BFM-associated genes were characterized by their higher than expected involvement in cellular, regulatory and immune system processes, and BFD-associated genes by their involvement in cellular, metabolic, and regulatory processes. Mendelian Randomization analyses suggested that the gene expression of 69 genes was causally related to BFM and BFD. Six genes were replicated in UK Biobank. In this study, we identified novel genes for BFM and BFD that are BFM- and BFD-specific, involved in different molecular processes, and whose up-/downregulated gene expression may causally contribute to obesity.

Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies

Pathogenic variants in genes that cause dilated cardiomyopathy (DCM) and arrhythmogenic cardiomyopathy (ACM) convey high risks for the development of heart failure through unknown mechanisms. Using single-nucleus RNA sequencing, we characterized the transcriptome of 880,000 nuclei from 18 control and 61 failing, nonischemic human hearts with pathogenic variants in DCM and ACM genes or idiopathic disease. We performed genotype-stratified analyses of the ventricular cell lineages and transcriptional states. The resultant DCM and ACM ventricular cell atlas demonstrated distinct right and left ventricular responses, highlighting genotype-associated pathways, intercellular interactions, and differential gene expression at single-cell resolution. Together, these data illuminate both shared and distinct cellular and molecular architectures of human heart failure and suggest candidate therapeutic targets.

Ablation of lysophosphatidic acid receptor 1 attenuates hypertrophic cardiomyopathy in a mouse model

Myocardial fibrosis is a key pathologic feature of hypertrophic cardiomyopathy (HCM). However, the fibrotic pathways activated by HCM-causing sarcomere protein gene mutations are poorly defined. Because lysophosphatidic acid is a mediator of fibrosis in multiple organs and diseases, we tested the role of the lysophosphatidic acid pathway in HCM. Lysphosphatidic acid receptor 1 (LPAR1), a cell surface receptor, is required for lysophosphatidic acid mediation of fibrosis. We bred HCM mice carrying a pathogenic myosin heavy-chain variant (403+/-) with Lpar1-ablated mice to create mice carrying both genetic changes (403+/- LPAR1 -/-) and assessed development of cardiac hypertrophy and fibrosis. Compared with 403+/- LPAR1WT, 403+/- LPAR1 -/- mice developed significantly less hypertrophy and fibrosis. Single-nucleus RNA sequencing of left ventricular tissue demonstrated that Lpar1 was predominantly expressed by lymphatic endothelial cells (LECs) and cardiac fibroblasts. Lpar1 ablation reduced the population of LECs, confirmed by immunofluorescence staining of the LEC markers Lyve1 and Ccl21a and, by in situ hybridization, for Reln and Ccl21a. Lpar1 ablation also altered the distribution of fibroblast cell states. FB1 and FB2 fibroblasts decreased while FB0 and FB3 fibroblasts increased. Our findings indicate that Lpar1 is expressed predominantly by LECs and fibroblasts in the heart and is required for development of hypertrophy and fibrosis in an HCM mouse model. LPAR1 antagonism, including agents in clinical trials for other fibrotic diseases, may be beneficial for HCM.

A trans locus causes a ribosomopathy in hypertrophic hearts that affects mRNA translation in a protein length-dependent fashion

BACKGROUND

Little is known about the impact of trans-acting genetic variation on the rates with which proteins are synthesized by ribosomes. Here, we investigate the influence of such distant genetic loci on the efficiency of mRNA translation and define their contribution to the development of complex disease phenotypes within a panel of rat recombinant inbred lines.

RESULTS

We identify several tissue-specific master regulatory hotspots that each control the translation rates of multiple proteins. One of these loci is restricted to hypertrophic hearts, where it drives a translatome-wide and protein length-dependent change in translational efficiency, altering the stoichiometric translation rates of sarcomere proteins. Mechanistic dissection of this locus across multiple congenic lines points to a translation machinery defect, characterized by marked differences in polysome profiles and misregulation of the small nucleolar RNA SNORA48. Strikingly, from yeast to humans, we observe reproducible protein length-dependent shifts in translational efficiency as a conserved hallmark of translation machinery mutants, including those that cause ribosomopathies. Depending on the factor mutated, a pre-existing negative correlation between protein length and translation rates could either be enhanced or reduced, which we propose to result from mRNA-specific imbalances in canonical translation initiation and reinitiation rates.

CONCLUSIONS

We show that distant genetic control of mRNA translation is abundant in mammalian tissues, exemplified by a single genomic locus that triggers a translation-driven molecular mechanism. Our work illustrates the complexity through which genetic variation can drive phenotypic variability between individuals and thereby contribute to complex disease.

Impairment of the ER/mitochondria compartment in human cardiomyocytes with PLN p.Arg14del mutation

The phospholamban (PLN) p.Arg14del mutation causes dilated cardiomyopathy, with the molecular disease mechanisms incompletely understood. Patient dermal fibroblasts were reprogrammed to hiPSC, isogenic controls were established by CRISPR/Cas9, and cardiomyocytes were differentiated. Mutant cardiomyocytes revealed significantly prolonged Ca2+ transient decay time, Ca2+ -load dependent irregular beating pattern, and lower force. Proteomic analysis revealed less endoplasmic reticulum (ER) and ribosomal and mitochondrial proteins. Electron microscopy showed dilation of the ER and large lipid droplets in close association with mitochondria. Follow-up experiments confirmed impairment of the ER/mitochondria compartment. PLN p.Arg14del end-stage heart failure samples revealed perinuclear aggregates positive for ER marker proteins and oxidative stress in comparison with ischemic heart failure and non-failing donor heart samples. Transduction of PLN p.Arg14del EHTs with the Ca2+ -binding proteins GCaMP6f or parvalbumin improved the disease phenotype. This study identified impairment of the ER/mitochondria compartment without SR dysfunction as a novel disease mechanism underlying PLN p.Arg14del cardiomyopathy. The pathology was improved by Ca2+ -scavenging, suggesting impaired local Ca2+ cycling as an important disease culprit.

Functional analysis of a gene-edited mouse model to gain insights into the disease mechanisms of a titin missense variant

Titin truncating variants are a well-established cause of cardiomyopathy; however, the role of titin missense variants is less well understood. Here we describe the generation of a mouse model to investigate the underlying disease mechanism of a previously reported titin A178D missense variant identified in a family with non-compaction and dilated cardiomyopathy. Heterozygous and homozygous mice carrying the titin A178D missense variant were characterised in vivo by echocardiography. Heterozygous mice had no detectable phenotype at any time point investigated (up to 1 year). By contrast, homozygous mice developed dilated cardiomyopathy from 3 months. Chronic adrenergic stimulation aggravated the phenotype. Targeted transcript profiling revealed induction of the foetal gene programme and hypertrophic signalling pathways in homozygous mice, and these were confirmed at the protein level. Unsupervised proteomics identified downregulation of telethonin and four-and-a-half LIM domain 2, as well as the upregulation of heat shock proteins and myeloid leukaemia factor 1. Loss of telethonin from the cardiac Z-disc was accompanied by proteasomal degradation; however, unfolded telethonin accumulated in the cytoplasm, leading to a proteo-toxic response in the mice.We show that the titin A178D missense variant is pathogenic in homozygous mice, resulting in cardiomyopathy. We also provide evidence of the disease mechanism: because the titin A178D variant abolishes binding of telethonin, this leads to its abnormal cytoplasmic accumulation. Subsequent degradation of telethonin by the proteasome results in proteasomal overload, and activation of a proteo-toxic response. The latter appears to be a driving factor for the cardiomyopathy observed in the mouse model.

Phosphatidylinositol 4-kinase β mutations cause nonsyndromic sensorineural deafness and inner ear malformation

Congenital hearing loss is a common disorder worldwide. Heterogeneous gene variation accounts for approximately 20-25% of such patients. We investigated a five-generation Chinese family with autosomal-dominant nonsyndromic sensorineural hearing loss (SNHL). No wave was detected in the pure-tone audiometry, and the auditory brainstem response was absent in all patients. Computed tomography of the patients, as well as of two sporadic SNHL cases, showed bilateral inner ear anomaly, cochlear maldevelopment, absence of the osseous spiral lamina, and an enlarged vestibular aqueduct. Such findings were absent in nonaffected persons. We used linkage analysis and exome sequencing and uncovered a heterozygous missense mutation in the PI4KB gene (p.Gln121Arg) encoding phosphatidylinositol 4-kinase β (PI4KB) from the patients in this family. In addition, 3 missense PI4KB (p.Val434Gly, p.Glu667Lys, and p.Met739Arg) mutations were identified in five patients with nonsyndromic SNHL from 57 sporadic cases. No such mutations were present within 600 Chinese controls, the 1000 genome project, gnomAD, or similar databases. Depleting pi4kb mRNA expression in zebrafish caused inner ear abnormalities and audiosensory impairment, mimicking the patient phenotypes. Moreover, overexpression of 4 human missense PI4KB mutant mRNAs in zebrafish embryos resulted in impaired hearing function, suggesting dominant-negative effects. Taken together, our results reveal that PI4KB mutations can cause SNHL and inner ear malformation. PI4KB should be included in neonatal deafness screening.

Cells of the adult human heart

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

Statin-induced myopathic changes in primary human muscle cells and reversal by a prostaglandin F2 alpha analogue

Statin-related muscle side effects are a constant healthcare problem since patient compliance is dependent on side effects. Statins reduce plasma cholesterol levels and can prevent secondary cardiovascular diseases. Although statin-induced muscle damage has been studied, preventive or curative therapies are yet to be reported. We exposed primary human muscle cell populations (n = 22) to a lipophilic (simvastatin) and a hydrophilic (rosuvastatin) statin and analyzed their expressome. Data and pathway analyses included GOrilla, Reactome and DAVID. We measured mevalonate intracellularly and analyzed eicosanoid profiles secreted by human muscle cells. Functional assays included proliferation and differentiation quantification. More than 1800 transcripts and 900 proteins were differentially expressed after exposure to statins. Simvastatin had a stronger effect on the expressome than rosuvastatin, but both statins influenced cholesterol biosynthesis, fatty acid metabolism, eicosanoid synthesis, proliferation, and differentiation of human muscle cells. Cultured human muscle cells secreted ω-3 and ω-6 derived eicosanoids and prostaglandins. The ω-6 derived metabolites were found at higher levels secreted from simvastatin-treated primary human muscle cells. Eicosanoids rescued muscle cell differentiation. Our data suggest a new aspect on the role of skeletal muscle in cholesterol metabolism. For clinical practice, the addition of omega-n fatty acids might be suitable to prevent or treat statin-myopathy.

Widespread Translational Control of Fibrosis in the Human Heart by RNA-Binding Proteins

Supplemental Digital Content is available in the text.

Fibrosis is a common pathology in many cardiac disorders and is driven by the activation of resident fibroblasts. The global posttranscriptional mechanisms underlying fibroblast-to-myofibroblast conversion in the heart have not been explored.

IL-11 is a crucial determinant of cardiovascular fibrosis

Fibrosis is a common pathology in cardiovascular disease. In the heart, fibrosis causes mechanical and electrical dysfunction and in the kidney, it predicts the onset of renal failure. Transforming growth factor β1 (TGFβ1) is the principal pro-fibrotic factor, but its inhibition is associated with side effects due to its pleiotropic roles. We hypothesized that downstream effectors of TGFβ1 in fibroblasts could be attractive therapeutic targets and lack upstream toxicity. Here we show, using integrated imaging-genomics analyses of primary human fibroblasts, that upregulation of interleukin-11 (IL-11) is the dominant transcriptional response to TGFβ1 exposure and required for its pro-fibrotic effect. IL-11 and its receptor (IL11RA) are expressed specifically in fibroblasts, in which they drive non-canonical, ERK-dependent autocrine signalling that is required for fibrogenic protein synthesis. In mice, fibroblast-specific Il11 transgene expression or Il-11 injection causes heart and kidney fibrosis and organ failure, whereas genetic deletion of Il11ra1 protects against disease. Therefore, inhibition of IL-11 prevents fibroblast activation across organs and species in response to a range of important pro-fibrotic stimuli. These results reveal a central role of IL-11 in fibrosis and we propose that inhibition of IL-11 is a potential therapeutic strategy to treat fibrotic diseases.

Tumour ischaemia by interferon-γ resembles physiological blood vessel regression

The relative contribution of the effector molecules produced by T cells to tumour rejection is unclear, but interferon-γ (IFNγ) is critical in most of the analysed models. Although IFNγ can impede tumour growth by acting directly on cancer cells, it must also act on the tumour stroma for effective rejection of large, established tumours. However, which stroma cells respond to IFNγ and by which mechanism IFNγ contributes to tumour rejection through stromal targeting have remained unknown. Here we use a model of IFNγ induction and an IFNγ-GFP fusion protein in large, vascularized tumours growing in mice that express the IFNγ receptor exclusively in defined cell types. Responsiveness to IFNγ by myeloid cells and other haematopoietic cells, including T cells or fibroblasts, was not sufficient for IFNγ-induced tumour regression, whereas responsiveness of endothelial cells to IFNγ was necessary and sufficient. Intravital microscopy revealed IFNγ-induced regression of the tumour vasculature, resulting in arrest of blood flow and subsequent collapse of tumours, similar to non-haemorrhagic necrosis in ischaemia and unlike haemorrhagic necrosis induced by tumour necrosis factor. The early events of IFNγ-induced tumour ischaemia resemble non-apoptotic blood vessel regression during development, wound healing or IFNγ-mediated, pregnancy-induced remodelling of uterine arteries. A better mechanistic understanding of how solid tumours are rejected may aid the design of more effective protocols for adoptive T-cell therapy.

A roadmap of constitutive NF-κB activity in Hodgkin lymphoma: Dominant roles of p50 and p52 revealed by genome-wide analyses

Background

NF-κB is widely involved in lymphoid malignancies; however, the functional roles and specific transcriptomes of NF-κB dimers with distinct subunit compositions have been unclear.

Methods

Using combined ChIP-sequencing and microarray analyses, we determined the cistromes and target gene signatures of canonical and non-canonical NF-κB species in Hodgkin lymphoma (HL) cells.

Results

We found that the various NF-κB subunits are recruited to regions with redundant κB motifs in a large number of genes. Yet canonical and non-canonical NF-κB dimers up- and downregulate gene sets that are both distinct and overlapping, and are associated with diverse biological functions. p50 and p52 are formed through NIK-dependent p105 and p100 precursor processing in HL cells and are the predominant DNA binding subunits. Logistic regression analyses of combinations of the p50, p52, RelA, and RelB subunits in binding regions that have been assigned to genes they regulate reveal a cross-contribution of p52 and p50 to canonical and non-canonical transcriptomes. These analyses also indicate that the subunit occupancy pattern of NF-κB binding regions and their distance from the genes they regulate are determinants of gene activation versus repression. The pathway-specific signatures of activated and repressed genes distinguish HL from other NF-κB-associated lymphoid malignancies and inversely correlate with gene expression patterns in normal germinal center B cells, which are presumed to be the precursors of HL cells.

Conclusions

We provide insights that are relevant for lymphomas with constitutive NF-κB activation and generally for the decoding of the mechanisms of differential gene regulation through canonical and non-canonical NF-κB signaling.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0280-5) contains supplementary material, which is available to authorized users.

Natural variation of histone modification and its impact on gene expression in the rat genome

Histone modifications are epigenetic marks that play fundamental roles in many biological processes including the control of chromatin-mediated regulation of gene expression. Little is known about interindividual variability of histone modification levels across the genome and to what extent they are influenced by genetic variation. We annotated the rat genome with histone modification maps, identified differences in histone trimethyl-lysine levels among strains, and described their underlying genetic basis at the genome-wide scale using ChIP-seq in heart and liver tissues in a panel of rat recombinant inbred and their progenitor strains. We identified extensive variation of histone methylation levels among individuals and mapped hundreds of underlying cis- and trans-acting loci throughout the genome that regulate histone methylation levels in an allele-specific manner. Interestingly, most histone methylation level variation was trans-linked and the most prominent QTL identified influenced H3K4me3 levels at 899 putative promoters throughout the genome in the heart. Cis- acting variation was enriched in binding sites of distinct transcription factors in heart and liver. The integrated analysis of DNA variation together with histone methylation and gene expression levels showed that histoneQTLs are an important predictor of gene expression and that a joint analysis significantly enhanced the prediction of gene expression traits (eQTLs). Our data suggest that genetic variation has a widespread impact on histone trimethylation marks that may help to uncover novel genotype-phenotype relationships.

Increased proliferative cells in the medullary thick ascending limb of the loop of Henle in the Dahl salt-sensitive rat

Studies of transcriptome profiles have provided new insights into mechanisms underlying the development of hypertension. Cell type heterogeneity in tissue samples, however, has been a significant hindrance in these studies. We performed a transcriptome analysis in medullary thick ascending limbs of the loop of Henle isolated from Dahl salt-sensitive rats. Genes differentially expressed between Dahl salt-sensitive rats and salt-insensitive consomic SS.13(BN) rats on either 0.4% or 7 days of 8.0% NaCl diet (n=4) were highly enriched for genes located on chromosome 13, the chromosome substituted in the SS.13(BN) rat. A pathway involving cell proliferation and cell cycle regulation was identified as one of the most highly ranked pathways based on differentially expressed genes and by a Bayesian model analysis. Immunofluorescent analysis indicated that just 1 week of a high-salt diet resulted in a severalfold increase in proliferative medullary thick ascending limb cells in both rat strains, and that Dahl salt-sensitive rats exhibited a significantly greater proportion of medullary thick ascending limb cells in a proliferative state than in SS.13(BN) rats (15.0±1.4% versus 10.1±0.6%; n=7-9; P<0.05). The total number of cells per medullary thick ascending limb section analyzed was not different between the 2 strains. The study revealed alterations in regulatory pathways in Dahl salt-sensitive rats in tissues highly enriched for a single cell type, leading to the unexpected finding of a greater increase in the number of proliferative medullary thick ascending limb cells in Dahl salt-sensitive rats on a high-salt diet.

The genome sequence of the spontaneously hypertensive rat: Analysis and functional significance

The spontaneously hypertensive rat (SHR) is the most widely studied animal model of hypertension. Scores of SHR quantitative loci (QTLs) have been mapped for hypertension and other phenotypes. We have sequenced the SHR/OlaIpcv genome at 10.7-fold coverage by paired-end sequencing on the Illumina platform. We identified 3.6 million high-quality single nucleotide polymorphisms (SNPs) between the SHR/OlaIpcv and Brown Norway (BN) reference genome, with a high rate of validation (sensitivity 96.3%-98.0% and specificity 99%-100%). We also identified 343,243 short indels between the SHR/OlaIpcv and reference genomes. These SNPs and indels resulted in 161 gain or loss of stop codons and 629 frameshifts compared with the BN reference sequence. We also identified 13,438 larger deletions that result in complete or partial absence of 107 genes in the SHR/OlaIpcv genome compared with the BN reference and 588 copy number variants (CNVs) that overlap with the gene regions of 688 genes. Genomic regions containing genes whose expression had been previously mapped as cis-regulated expression quantitative trait loci (eQTLs) were significantly enriched with SNPs, short indels, and larger deletions, suggesting that some of these variants have functional effects on gene expression. Genes that were affected by major alterations in their coding sequence were highly enriched for genes related to ion transport, transport, and plasma membrane localization, providing insights into the likely molecular and cellular basis of hypertension and other phenotypes specific to the SHR strain. This near complete catalog of genomic differences between two extensively studied rat strains provides the starting point for complete elucidation, at the molecular level, of the physiological and pathophysiological phenotypic differences between individuals from these strains.

Air transmissivity of feathers

The flight feathers and their coverts of a European kestrel Falco tinnunculus have been tested for their air transmissivity. The transmissivity was measured in both directions, i.e. from ventral to dorsal and vice versa; the mean difference between the two directions was less than 10 %. However, the transmissivity of the inner versus outer vanes of the remiges and coverts differed significantly, with the outer vanes being more transmissive than the inner vanes. A functional interpretation of the different transmissivities of the inner and outer vanes is given, and we propose that its significance lies in the formation of a smooth, continuous wing surface. The individual feathers are pushed firmly towards one another as a result of the different transmissivities, which cause a pressure gradient to build up from the less-transmissive inner vane towards the overlying, more-transmissive outer vane of the adjacent feather. In another test series, the transmissivity of the flight and covert feathers of 27 species was measured; the differences found between species were small.