Search

BACKGROUND: Mutation of the PRDM16 gene has been associated with human cardiomyopathy. The PRDM16 protein is a transcriptional regulator affecting cardiac development via Tbx5 and Hand1 regulating myocardial structure. Biallelic Prdm16 inactivation induces severe cardiac dysfunction with postnatal lethality and hypertrophy in mice. Early pathological events upon Prdm16 inactivation have not been explored. METHODS: This study performed in depth pathophysiological and molecular analysis of male and female Prdm16csp1/wt mice carrying systemic, monoallelic Prdm16 gene inactivation. We systematically assessed early molecular changes with transcriptomics, proteomics, and metabolomics. Kinetic modelling of the cardiac metabolism was undertaken in silico with CARDIOKIN. RESULTS: Prdm16csp1/wt mice are viable up to 8 months, develop hypoplastic hearts, and diminished systolic performance that is more pronounced in female mice. Prdm16csp1/wt hearts demonstrate moderate alterations of specific transcripts and protein levels with consistent upregulation of pyridine nucleotide-disulphide oxidoreductase domain 2 (Pyroxd2) and the transcriptional regulator pre B-cell leukemia transcription factor interacting protein 1 (Pbxip1). The strongest concordant transcriptional upregulation was detected for Prdm16 itself probably by an autoregulatory mechanism. Prdm16csp1/wt cardiac tissue showed reduction of metabolites associated with amino acid as well as glycerol metabolism, glycolysis, and tricarboxylic acid cycle. Global lipid metabolism was also affected with accumulation of triacylglycerides detected in male Prdm16csp1/wt hearts. In addition, Prdm16csp1/wt cardiac tissue revealed diminished glutathione (GSH) and increased inosine monophosphate (IMP) levels indicating oxidative stress and a dysregulated energetics, respectively. Metabolic modelling in silico suggested lowered fatty acid utilization in male and reduced glucose utilization in female Prdm16csp1/wt cardiac tissue. CONCLUSIONS: Monoallelic Prdm16 mutation restricts cardiac performance in Prdm16csp1/wt mice.Metabolic alterations precede transcriptional dysregulation in Prdm16csp1/wt cardiac tissue. Female Prdm16csp1/wt mice develop a more pronounced phenotype indicating a sexual dimorphism at this early pathological window. This study suggests that metabolic dysregulation is an early event in PRDM16 associated cardiac pathology.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade., The LifeTime initiative is an ambitious, multidisciplinary programme that aims to improve healthcare by tracking individual human cells during disease processes and responses to treatment in order to develop and implement cell-based interceptive medicine in Europe.

Ribosome profiling (Ribo-seq) has catalyzed a paradigm shift in our understanding of the translational ‘vocabulary’ of the human genome, discovering thousands of translated open reading frames (ORFs) within long non-coding RNAs and presumed untranslated regions of protein-coding genes. However, reference gene annotation projects have been circumspect in their incorporation of these ORFs due to uncertainties about their experimental reproducibility and physiological roles. Yet, it is indisputable that certain Ribo-seq ORFs make stable proteins, others mediate gene regulation, and many have medical implications. Ultimately, the absence of standardized ORF annotation has created a circular problem: while Ribo-seq ORFs remain unannotated by reference biological databases, this lack of characterisation will thwart research efforts examining their roles. Here, we outline the initial stages of a community-led effort supported by GENCODE / Ensembl, HGNC and UniProt to produce a consolidated catalog of human Ribo-seq ORFs.

RNA-binding proteins (RBPs) are key regulators of RNA metabolism. Many RBPs possess uncharacterized RNA-binding domains and localize to multiple subcellular compartments, suggesting their involvement in multiple biological processes. We searched for such multifunctionality within a set of 143 RBPs by integrating experimentally validated target genes with the transcriptomes and translatomes of 80 human hearts. This revealed that RBP abundance is predictive of the extent of target regulation in vivo, leading us to newly associate 27 RBPs with translational control. Amongst those were several splicing factors, of which the muscle specific RBM20 modulated target translation rates through switches in isoform production. For 21 RBPs, we newly observed dual regulatory effects impacting both mRNA levels and translation rates, albeit for virtually independent sets of target genes. We highlight a subset, including G3BP1, PUM1, UCHL5, and DDX3X, where dual regulation is achieved by differential affinity for targets of distinct length and functionality. Strikingly, in a manner very similar to DDX3X, the known splicing factors EFTUD2 and PRPF8 selectively influence target translation rates depending on 5’ UTR structure. Our results indicate unanticipated complexity of protein-RNA interactions at consecutive stages of gene expression and implicate multiple core splicing factors as key regulators of translational output.

In this Perspective, owing to an error in the HTML, the surname of author Alejandro López-Tobón of the LifeTime Community Working Groups consortium was indexed as ‘Tobon’ rather than ‘López-Tobón’ and the accents were missing. The HTML version of the original Perspective has been corrected; the PDF and print versions were always correct. © 2021, The Author(s).

Previous genome-wide association studies revealed multiple common variants involved in eczema but the role of rare variants remains to be elucidated. Here, we investigate the role of rare variants in eczema susceptibility. We meta-analyze 21 study populations including 20,016 eczema cases and 380,433 controls. Rare variants are imputed with high accuracy using large population-based reference panels. We identify rare exonic variants in DUSP1, NOTCH4, and SLC9A4 to be associated with eczema. In DUSP1 and NOTCH4 missense variants are predicted to impact conserved functional domains. In addition, five novel common variants at SATB1-AS1/KCNH8, TRIB1/LINC00861, ZBTB1, TBX21/OSBPL7, and CSF2RB are discovered. While genes prioritized based on rare variants are significantly up-regulated in the skin, common variants point to immune cell function. Over 20% of the single nucleotide variant-based heritability is attributable to rare and low-frequency variants. The identified rare/low-frequency variants located in functional protein domains point to promising targets for novel therapeutic approaches to eczema., Genetic studies of eczema to date have mostly explored common genetic variation. Here, the authors perform a large meta-analysis for common and rare variants and discover 8 loci associated with eczema. Over 20% of the heritability of the condition is attributable to rare variants.

The SARS-CoV-2 coronavirus, the etiologic agent responsible for COVID-19 coronavirus disease, is a global threat. To better understand viral tropism, we assessed the RNA expression of the coronavirus receptor, ACE2, as well as the viral S protein priming protease TMPRSS2 thought to govern viral entry in single-cell RNA-sequencing (scRNA-seq) datasets from healthy individuals generated by the Human Cell Atlas consortium. We found that ACE2, as well as the protease TMPRSS2, are differentially expressed in respiratory and gut epithelial cells. In-depth analysis of epithelial cells in the respiratory tree reveals that nasal epithelial cells, specifically goblet/secretory cells and ciliated cells, display the highest ACE2 expression of all the epithelial cells analyzed. The skewed expression of viral receptors/entry-associated proteins towards the upper airway may be correlated with enhanced transmissivity. Finally, we showed that many of the top genes associated with ACE2 airway epithelial expression are innate immune-associated, antiviral genes, highly enriched in the nasal epithelial cells. This association with immune pathways might have clinical implications for the course of infection and viral pathology, and highlights the specific significance of nasal epithelia in viral infection. Our findings underscore the importance of the availability of the Human Cell Atlas as a reference dataset. In this instance, analysis of the compendium of data points to a particularly relevant role for nasal goblet and ciliated cells as early viral targets and potential reservoirs of SARS-CoV-2 infection. This, in turn, serves as a biological framework for dissecting viral transmission and developing clinical strategies for prevention and therapy.

Background: Several long noncoding RNAs (lncRNAs) have been shown to function as central components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, their functions remain largely uncharacterized. This is compounded by the fact that identifying lncRNA loci that have robust and reproducible phenotypes when mutated has been a challenge. Results: We previously generated a cohort of 20 lncRNA loci knockout mice. Here, we extend our initial study and provide a more detailed analysis of the highly conserved lncRNA locus, Taurine Upregulated Gene 1 (Tug1). We report that Tug1 knockout male mice are sterile with complete penetrance due to a low sperm count and abnormal sperm morphology. Having identified a lncRNA loci with a robust phenotype, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region (DNA, RNA, protein, or the act of transcription) have activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors three distinct regulatory activities - two noncoding and one coding: (i) a cis DNA repressor that regulates many neighboring genes, (ii) a lncRNA that can regulate genes by a trans-based function, and finally (iii) Tug1 encodes an evolutionary conserved peptide that when overexpressed impacts mitochondrial membrane potential. Conclusions: Our results reveal an essential role for the Tug1 locus in male fertility and uncover three distinct regulatory activities in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Objective The spontaneously hypertensive rat strain is a frequently used disease model. In a previous study, we measured translational efficiency from this strain and BN-Lx animals. Here, we describe long RNA sequencing reads from ribosomal RNA depleted samples from the same animals. This data can be used to investigate splicing-related events. Results RNA was extracted from rat liver and heart left ventricle from BN-Lx and SHR/Ola rats in biological replicates. Ribosomal RNA was removed and the samples subjected to directional high-throughput RNA-sequencing. Read and alignment statistics indicate high quality of the data. The raw sequencing reads are freely available on the NCBI short read archive and can be used for further research on tissue and strain differences, or analysed together with other published high-throughput data from the same animals.

Abstract Background Fibrosis is a common pathology in many cardiac disorders and is driven by the activation of resident fibroblasts. The global post-transcriptional mechanisms underlying fibroblast-to-myofibroblast conversion in the heart have not been explored. Methods Genome-wide changes of RNA transcription and translation during human cardiac fibroblast activation were monitored with RNA sequencing and ribosome profiling. We then used miRNA-and RNA-binding protein-based analyses to identify translational regulators of fibrogenic genes. To reveal post-transcriptional mechanisms in the human fibrotic heart, we then integrated our findings with cardiac ribosome occupancy levels of 30 dilated cardiomyopathy patients. Results We generated nucleotide-resolution translatome data during the TGFβ1-driven cellular transition of human cardiac fibroblasts to myofibroblasts. This identified dynamic changes of RNA transcription and translation at several time points during the fibrotic response, revealing transient and early-responder genes. Remarkably, about one-third of all changes in gene expression in activated fibroblasts are subject to translational regulation and dynamic variation in ribosome occupancy affects protein abundance independent of RNA levels. Targets of RNA-binding proteins were strongly enriched in post-transcriptionally regulated genes, suggesting genes such as MBNL2 can act as translational activators or repressors. Ribosome occupancy in the hearts of patients with dilated cardiomyopathy suggested an extensive post-transcriptional regulatory network underlying cardiac fibrosis. Key network hubs include RNA-binding proteins such as PUM2 and QKI that work in concert to regulate the translation of target transcripts in human diseased hearts. Conclusions We reveal widespread translational effects of TGFβ1 and define novel post-transcriptional events that control the fibroblast-to-myofibroblast transition. Regulatory networks that affect ribosome occupancy in fibroblasts are paralleled in human heart disease. Our findings show the central importance of translational control in fibrosis and highlight novel pathogenic mechanisms in heart failure.

The inbred Brown Norway (BN) rat develops spontaneous ruptures of the internal elastic lamina (RIEL) of the abdominal aorta (AA) and iliac arteries. Prior studies with crosses of the BN/Orl RJ (susceptible) and LOU/M (resistant) showed the presence of a significant QTL on chromosome 5 and the production of congenic rats proved the involvement of this locus. In this study, we further dissected the above-mentioned QTL by creating a new panel of LOU.BN(chr5) congenic and subcongenic lines and reduced the locus to 5.2 Mb. Then we studied 1,002 heterogeneous stock (HS) rats, whose phenotyping revealed a low prevalence and high variability for RIEL. High-resolution mapping in the HS panel detected the major locus on chromosome 5 (log P > 35) and refined it to 1.4 Mb. Subsequently, RNA-seq analysis on AA of BN, congenics, and LOU revealed expression differences for only protease inhibitor 15 (Pi15) gene and a putative long intergenic noncoding RNA (lincRNA) within the linkage region. The high abundance of lincRNA with respect to reduced Pi15 expression, in conjunction with exertion of longitudinal strain, may be related to RIEL, indicating the potential importance of proteases in biological processes related to defective aortic internal elastic lamina structure. Similar mechanisms may be involved in aneurysm initiation in the human AA.

RATIONALE: Alternative cleavage and polyadenylation (APA) of mRNA represents a layer of gene regulation that to date has remained unexplored in the heart. This phenomenon may be very relevant, as the positioning of the polyA tail in mRNAs influences the length of the 3'UTR, a critical determinant of gene expression. OBJECTIVE: To investigate whether the 3'UTR length is regulated by APA in the human heart and whether this changes in the failing heart. METHODS AND RESULTS: We used 3'end RNA-sequencing (e3'-Seq) to directly measure global patterns of APA in healthy and failing human heart specimens. By monitoring polyadenylation profiles in these hearts, we identified disease-specific APA signatures in numerous genes. Interestingly, many of the genes with shortened 3'UTRs in heart failure were enriched for functional groups such as 'RNA binding', while genes with longer 3'UTRs were enriched for 'cytoskeletal organization' and 'actin binding'. RNA sequencing in a larger series of human hearts revealed that these APA candidates are often differentially expressed in failing hearts, with an inverse correlation between 3'UTR length and the level of gene expression. Protein levels of the APA regulator, Poly(A)-Binding Protein Nuclear 1 were substantially downregulated in failing hearts. CONCLUSIONS: We provide genome-wide, high-resolution polyadenylation maps of the human heart and show that the 3'end formation of mRNA is dynamic in heart failure, suggesting that APA-mediated 3'UTR length modulation represents an additional layer of gene regulation in failing hearts.

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.

Congenital anomalies of the kidney and urogenital tract (CAKUT) occur in approximately 0.5% of live births and represent the most frequent cause of end-stage renal disease in neonates and children. The genetic basis of CAKUT is not well defined. To understand more fully the genetic basis of one type of CAKUT, unilateral renal agenesis (URA), we are studying inbred ACI rats, which spontaneously exhibit URA and associated urogenital anomalies at an incidence of approximately 10%. URA is inherited as an incompletely dominant trait with incomplete penetrance in crosses between ACI and Brown Norway (BN) rats and a single responsible genetic locus, designated Renag1, was previously mapped to rat chromosome 14 (RNO14). The goals of this study were to fine map Renag1, identify the causal genetic variant responsible for URA, confirm that the Renag1 variant is the sole determinant of URA in the ACI rat, and define the embryologic basis of URA in this rat model. Data presented herein localize Renag1 to a 379 kilobase (kb) interval that contains a single protein coding gene, Kit (v-kit Hardy-Zukerman 4 feline sarcoma viral oncogene homolog); identify an endogenous retrovirus-derived long terminal repeat located within Kit intron 1 as the probable causal variant; demonstrate aberrant development of the nephric duct in the anticipated number of ACI rat embryos; and demonstrate expression of Kit and Kit ligand (Kitlg) in the nephric duct. Congenic rats that harbor ACI alleles at Renag1 on the BN genetic background exhibit the same spectrum of urogenital anomalies as ACI rats, indicating that Renag1 is necessary and sufficient to elicit URA and associated urogenital anomalies. These data reveal the first genetic link between Kit and URA and illustrate the value of the ACI rat as a model for defining the mechanisms and cell types in which Kit functions during urogenital development.

Background: ChIP-seq has become a routine method for interrogating the genome-wide distribution of various histone modifications. An important experimental goal is to compare the ChIP-seq profiles between an experimental sample and a reference sample, and to identify regions that show differential enrichment. However, comparative analysis of samples remains challenging for histone modifications with broad domains, such as heterochromatin-associated H3K27me3, as most ChIP-seq algorithms are designed to detect well defined peak-like features. Results: To address this limitation we introduce histoneHMM, a powerful bivariate Hidden Markov Model for the differential analysis of histone modifications with broad genomic footprints. histoneHMM aggregates short-reads over larger regions and takes the resulting bivariate read counts as inputs for an unsupervised classification procedure, requiring no further tuning parameters. histoneHMM outputs probabilistic classifications of genomic regions as being either modified in both samples, unmodified in both samples or differentially modified between samples. We extensively tested histoneHMM in the context of two broad repressive marks, H3K27me3 and H3K9me3, and evaluated region calls with follow up qPCR as well as RNA-seq data. Our results show that histoneHMM outperforms competing methods in detecting functionally relevant differentially modified regions. Conclusion: histoneHMM is a fast algorithm written in C++ and compiled as an R package. It runs in the popular R computing environment and thus seamlessly integrates with the extensive bioinformatic tool sets available through Bioconductor. This makes histoneHMM an attractive choice for the differential analysis of ChIP-seq data. Software is available from http://histonehmm.molgen.mpg.de webcite.

Mutations in several genes cause non-autoimmune diabetes, but numerous patients still have unclear genetic defects, hampering our understanding of the development of the disease and preventing pathogenesis-oriented treatment. We used whole-genome sequencing with linkage analysis to study a consanguineous family with early-onset antibody-negative diabetes and identified a novel deletion in PCBD1 (pterin-4 alpha-carbinolamine dehydratase/dimerization cofactor of hepatocyte nuclear factor-1 alpha), a gene that was recently proposed as a likely cause of diabetes. A subsequent reevaluation of patients with mild neonatal hyperphenylalaninemia due to mutations in PCBD1 from the BIODEF database identified three additional patients who had developed HNF1A-like diabetes in puberty, indicating early β-cell failure. We found that Pcbd1 is expressed in the developing pancreas of both mouse and Xenopus embryos from early specification onward showing colocalization with insulin. Importantly, a morpholino-mediated knockdown in Xenopus revealed that pcbd1 activity is required for the proper establishment of early pancreatic fate within the endoderm. We provide the first genetic evidence that PCBD1 mutations can cause early-onset non-autoimmune diabetes with features similar to dominantly inherited HNF1A-diabetes. This condition responds to and can be treated with oral drugs instead of insulin, which is important clinical information for these patients. Finally, patients at risk can be detected through a newborn screening for phenylketonuria.

Summary A rare germline duplication upstream of the bone morphogenetic protein antagonist GREM1 causes a Mendelian-dominant predisposition to colorectal cancer (CRC). The underlying disease mechanism is strong, ectopic GREM1 overexpression in the intestinal epithelium. Here, we confirm that a common GREM1 polymorphism, rs16969681, is also associated with CRC susceptibility, conferring ∼20% differential risk in the general population. We hypothesized the underlying cause to be moderate differences in GREM1 expression. We showed that rs16969681 lies in a region of active chromatin with allele- and tissue-specific enhancer activity. The CRC high-risk allele was associated with stronger gene expression, and higher Grem1 mRNA levels increased the intestinal tumor burden in ApcMin mice. The intestine-specific transcription factor CDX2 and Wnt effector TCF7L2 bound near rs16969681, with significantly higher affinity for the risk allele, and CDX2 overexpression in CDX2/GREM1-negative cells caused re-expression of GREM1. rs16969681 influences CRC risk through effects on Wnt-driven GREM1 expression in colorectal tumors., Graphical Abstract, Highlights • Polymorphism rs16969681 near GREM1 is associated with bowel cancer risk in humans • rs16969681 is in a region with allele- and tissue-specific enhancer activity • The enhancer is bound by CDX2 and TCF7L2; binding is stronger to the high-risk allele • The resulting small differences in GREM1 expression modulate bowel tumorigenesis, Lewis et al. investigate how a common polymorphism, rs16969681, near the GREM1 gene predisposes to bowel cancer. They show that rs16969681 lies at an enhancer of GREM1 mRNA expression. The rs16969681 alleles differentially bind the intestine-specific transcription factor CDX2 and the Wnt effector TCF7L2, with the cancer risk allele showing stronger binding. This results in differences in GREM1 expression among individuals with the three rs16969681 genotypes. Different levels of GREM1 modulate the growth of colorectal tumors.

Rat hypodactyly (hd) is an autosomal recessive mutation manifesting in homozygotes as reduction or loss of digits II and III. We mapped the hd allele to a short segment of chromosome 10, containing 16 genes. None of these genes has been shown to influence limb development yet. In situ hybridization showed no changes in several important patterning genes (Shh, Fgf8, Bmp2, 4, 7). However, we found that expression of cartilage condensation marker Sox9, and Bmp receptor Bmpr1b (acting as an upstream activator of Sox9 expression) is absent from the subepithelial mesenchyme of the digit condensations II and III. The failure of the chondrogenic condensations to extend towards the subepithelial mesenchyme may reduce the size of digit primordia and underlie the subsequent loss of phalanges and reduction of metacarpals/metatarsals in hd rats.

The pathophysiological basis of Liddle's syndrome, a rare autosomal dominant form of arterial hypertension, has been found to rest on missense mutations or truncations of the beta- and gamma-subunits of the epithelial sodium channel. The hypothesis has been advanced that molecular variants of these genes might also contribute to the common polygenic forms of hypertension. We tested this hypothesis by performing a cosegregation study in a reciprocal cross between the stroke-prone spontaneously hypertensive rat (SHRSPHD) and a Wistar-Kyoto rat (WKY-1HD) reference strain. We carried out genetic mapping and chromosomal assignment of the alpha-, beta-, and gamma-subunits of the epithelial sodium channel using both linkage analysis and fluorescent in situ hybridization techniques. We demonstrate that in the rat, the beta- and gamma-subunits, as in humans, are in close linkage; they map to rat chromosome 1 and cosegregate with systolic pressure after dietary NaCl (logarithm of the odds [LOD] score, 3.7), although the peak LOD score of 5.0 for this quantitative trait locus was detected 4.4 cM away from the beta-/gamma-subunit locus. The alpha-subunit was mapped to chromosome 4 and exhibited no linkage to blood pressure phenotype. Comparative analysis of the complete coding sequences of all three subunits in the SHRSPHD and WKY-1HD strains revealed no biologically relevant mutations. Furthermore, Northern blot comparison of mRNA levels for all three subunits in the kidney showed no differences between SHRSPHD and WKY-1HD. Our results fail to support a material contribution of the epithelial sodium channel genes to blood pressure regulation in this model of polygenic hypertension.