Ian M. Silverman

Loss of NAD Homeostasis Leads to Progressive and Reversible Degeneration of Skeletal Muscle

NAD is an obligate co-factor for the catabolism of metabolic fuels in all cell types. However, the availability of NAD in several tissues can become limited during genotoxic stress and the course of natural aging. The point at which NAD restriction imposes functional limitations on tissue physiology remains unknown. We examined this question in murine skeletal muscle by specifically depleting Nampt, an essential enzyme in the NAD salvage pathway. Knockout mice exhibited a dramatic 85% decline in intramuscular NAD content, accompanied by fiber degeneration and progressive loss of both muscle strength and treadmill endurance. Administration of the NAD precursor nicotinamide riboside rapidly ameliorated functional deficits and restored muscle mass despite having only a modest effect on the intramuscular NAD pool. Additionally, lifelong overexpression of Nampt preserved muscle NAD levels and exercise capacity in aged mice, supporting a critical role for tissue-autonomous NAD homeostasis in maintaining muscle mass and function.

HAMR: high-throughput annotation of modified ribonucleotides

RNA is often altered post-transcriptionally by the covalent modification of particular nucleotides; these modifications are known to modulate the structure and activity of their host RNAs. The recent discovery that an RNA methyl-6 adenosine demethylase (FTO) is a risk gene in obesity has brought to light the significance of RNA modifications to human biology. These noncanonical nucleotides, when converted to cDNA in the course of RNA sequencing, can produce sequence patterns that are distinguishable from simple base-calling errors. To determine whether these modifications can be detected in RNA sequencing data, we developed a method that can not only locate these modifications transcriptome-wide with single nucleotide resolution, but can also differentiate between different classes of modifications. Using small RNA-seq data we were able to detect 92% of all known human tRNA modification sites that are predicted to affect RT activity. We also found that different modifications produce distinct patterns of cDNA sequence, allowing us to differentiate between two classes of adenosine and two classes of guanine modifications with 98% and 79% accuracy, respectively. To show the robustness of this method to sample preparation and sequencing methods, as well as to organismal diversity, we applied it to a publicly available yeast data set and achieved similar levels of accuracy. We also experimentally validated two novel and one known 3-methylcytosine (3mC) sites predicted by HAMR in human tRNAs. Researchers can now use our method to identify and characterize RNA modifications using only RNA-seq data, both retrospectively and when asking questions specifically about modified RNA.

Global Analysis of the RNA-Protein Interaction and RNA Secondary Structure Landscapes of the Arabidopsis Nucleus

Posttranscriptional regulation in eukaryotes requires cis- and trans-acting features and factors including RNA secondary structure and RNA-binding proteins (RBPs). However, a comprehensive view of the structural and RBP interaction landscape of nuclear RNAs has yet to be compiled for any organism. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approaches to globally profile these features in Arabidopsis seedling nuclei in vivo. We reveal anticorrelated patterns of secondary structure and RBP binding throughout nuclear mRNAs that demarcate sites of alternative splicing and polyadenylation. We also uncover a collection of protein-bound sequence motifs, and identify their structural contexts, co-occurrences in transcripts encoding functionally related proteins, and interactions with putative RBPs. Finally, using these motifs, we find that the chloroplast RBP CP29A also interacts with nuclear mRNAs. In total, we provide a simultaneous view of the RNA secondary structure and RBP interaction landscapes in a eukaryotic nucleus.

RNase-mediated protein footprint sequencing reveals protein-binding sites throughout the human transcriptome

Although numerous approaches have been developed to map RNA-binding sites of individual RNA-binding proteins (RBPs), few methods exist that allow assessment of global RBP–RNA interactions. Here, we describe PIP-seq, a universal, high-throughput, ribonuclease-mediated protein footprint sequencing approach that reveals RNA-protein interaction sites throughout a transcriptome of interest. We apply PIP-seq to the HeLa transcriptome and compare binding sites found using different cross-linkers and ribonucleases. From this analysis, we identify numerous putative RBP-binding motifs, reveal novel insights into co-binding by RBPs, and uncover a significant enrichment for disease-associated polymorphisms within RBP interaction sites.

Tissue-Specific RNA-Seq in Human Evoked Inflammation Identifies Blood and Adipose LincRNA Signatures of Cardiometabolic Diseases

Objective—Inappropriate transcriptional activation of innate immunity is a pathological feature of several cardiometabolic disorders, but little is known about inflammatory modulation of long intergenic noncoding RNAs (lincRNAs) in disease-relevant human tissues. Approach and Results—We applied deep RNA sequencing (>500 million filtered reads per sample) to blood and adipose during low-dose experimental endotoxemia (lipopolysaccharide) in a healthy human, with targeted replication in separate individuals undergoing endotoxemia (n=6), to identify inflammatory lincRNAs. A subset of these lincRNAs was examined for expression in adipocytes and monocytes, modulation in adipose of obese humans, and overlap with genome-wide association study signals for inflammatory and cardiometabolic traits. Of a stringent set of 4284 lincRNAs, ≈11% to 22% were expressed with 201 and 56 lincRNAs modulated by lipopolysaccharide in blood or adipose, respectively. Tissue-specific expression of a subset of 6 lipopolysaccharide-lincRNAs was replicated with lipopolysaccharide modulation confirmed for all 3 expressed in blood and 2 of 4 expressed in adipose. The broader generalizability of findings in blood of subject A was confirmed by RNA sequencing in 7 additional subjects. We confirmed adipocytes and monocytes as potential cell-sources of selective lipopolysaccharide-regulated lincRNAs, and 2 of these, linc-DMRT2 (P=0.002) and linc-TP53I13 (P=0.01), were suppressed in adipose of obese humans. Finally, we provide examples of lipopolysaccharide-modulated lincRNAs that overlap single nucleotide polymorphisms that are associated with cardiometabolic traits. Conclusions—Our findings provide novel insights into tissue-level, inflammatory transcriptome regulation in cardiometabolic diseases. These are complementary to more usual approaches limited to interrogation of DNA variations.

Evaluating the Impact of Sequencing Depth on Transcriptome Profiling in Human Adipose

Recent advances in RNA sequencing (RNA-Seq) have enabled the discovery of novel transcriptomic variations that are not possible with traditional microarray-based methods. Tissue and cell specific transcriptome changes during pathophysiological stress in disease cases versus controls and in response to therapies are of particular interest to investigators studying cardiometabolic diseases. Thus, knowledge on the relationships between sequencing depth and detection of transcriptomic variation is needed for designing RNA-Seq experiments and for interpreting results of analyses. Using deeply sequenced Illumina HiSeq 2000 101 bp paired-end RNA-Seq data derived from adipose of a healthy individual before and after systemic administration of endotoxin (LPS), we investigated the sequencing depths needed for studies of gene expression and alternative splicing (AS). In order to detect expressed genes and AS events, we found that ∼100 to 150 million (M) filtered reads were needed. However, the requirement on sequencing depth for the detection of LPS modulated differential expression (DE) and differential alternative splicing (DAS) was much higher. To detect 80% of events, ∼300 M filtered reads were needed for DE analysis whereas at least 400 M filtered reads were necessary for detecting DAS. Although the majority of expressed genes and AS events can be detected with modest sequencing depths (∼100 M filtered reads), the estimated gene expression levels and exon/intron inclusion levels were less accurate. We report the first study that evaluates the relationship between RNA-Seq depth and the ability to detect DE and DAS in human adipose. Our results suggest that a much higher sequencing depth is needed to reliably identify DAS events than for DE genes.

Genomic era analyses of RNA secondary structure and RNA-binding proteins reveal their significance to post-transcriptional regulation in plants.

The eukaryotic transcriptome is regulated both transcriptionally and post-transcriptionally. Transcriptional control was the major focus of early research efforts, while more recently post-transcriptional mechanisms have gained recognition for their significant regulatory importance. At the heart of post-transcriptional regulatory pathways are cis- and trans-acting features and factors including RNA secondary structure as well as RNA-binding proteins and their recognition sites on target RNAs. Recent advances in genomic methodologies have significantly improved our understanding of both RNA secondary structure and RNA-binding proteins and their regulatory effects within the eukaryotic transcriptome. In this review, we focus specifically on the collection of these regulatory moieties in plant transcriptomes. We describe the approaches for studying RNA secondary structure and RNA-protein interaction sites, with an emphasis on recent methodological advances that produce transcriptome-wide datasets. We discuss how these methods that include genome-wide RNA secondary structure determination and RNA-protein interaction site mapping are significantly improving our understanding of the functions of these two elements in post-transcriptional regulation. Finally, we delineate the need for additional genome-wide studies of RNA secondary structure and RNA-protein interactions in plants.

The Poly(C) Binding Protein Pcbp2 and Its Retrotransposed Derivative Pcbp1 Are Independently Essential to Mouse Development

ABSTRACT RNA-binding proteins participate in a complex array of posttranscriptional controls essential to cell type specification and somatic development. Despite their detailed biochemical characterizations, the degree to which each RNA-binding protein impacts mammalian embryonic development remains incompletely defined, and the level of functional redundancy among subsets of these proteins remains open to question. The poly(C) binding proteins, PCBPs (αCPs and hnRNP E proteins), are encoded by a highly conserved and broadly expressed gene family. The two major Pcbp isoforms, Pcbp2 and Pcbp1, are robustly expressed in a wide range of tissues and exert both nuclear and cytoplasmic controls over gene expression. Here, we report that Pcbp1-null embryos are rendered nonviable in the peri-implantation stage. In contrast, Pcbp2-null embryos undergo normal development until midgestation (12.5 to 13.5 days postcoitum), at which time they undergo a dramatic loss in viability associated with combined cardiovascular and hematopoietic abnormalities. Mice heterozygous for either Pcbp1 or Pcbp2 null alleles display a mild and nondisruptive defect in initial postpartum weight gain. These data reveal that Pcbp1 and Pcbp2 are individually essential for mouse embryonic development and have distinct impacts on embryonic viability and that Pcpb2 has a nonredundant in vivo role in hematopoiesis. These data further provide direct evidence that Pcbp1, a retrotransposed derivative of Pcpb2, has evolved an essential function(s) in the mammalian genome.

Chemical Modifications Mark Alternatively Spliced and Uncapped Messenger RNAs in Arabidopsis

Global identification of RNA modifications that affect the Watson-Crick base-pairing edge across the Arabidopsis transcriptome provides insights into the many functions of these important additions. Posttranscriptional chemical modification of RNA bases is a widespread and physiologically relevant regulator of RNA maturation, stability, and function. While modifications are best characterized in short, noncoding RNAs such as tRNAs, growing evidence indicates that mRNAs and long noncoding RNAs (lncRNAs) are likewise modified. Here, we apply our high-throughput annotation of modified ribonucleotides (HAMR) pipeline to identify and classify modifications that affect Watson-Crick base pairing at three different levels of the Arabidopsis thaliana transcriptome (polyadenylated, small, and degrading RNAs). We find this type of modifications primarily within uncapped, degrading mRNAs and lncRNAs, suggesting they are the cause or consequence of RNA turnover. Additionally, modifications within stable mRNAs tend to occur in alternatively spliced introns, suggesting they regulate splicing. Furthermore, these modifications target mRNAs with coherent functions, including stress responses. Thus, our comprehensive analysis across multiple RNA classes yields insights into the functions of covalent RNA modifications in plant transcriptomes.

Terpene metabolic engineering via nuclear or chloroplast genomes profoundly and globally impacts off-target pathways through metabolite signalling

Summary The impact of metabolic engineering on nontarget pathways and outcomes of metabolic engineering from different genomes are poorly understood questions. Therefore, squalene biosynthesis genes FARNESYL DIPHOSPHATE SYNTHASE (FPS) and SQUALENE SYNTHASE (SQS) were engineered via the Nicotiana tabacum chloroplast (C), nuclear (N) or both (CN) genomes to promote squalene biosynthesis. SQS levels were ~4300‐fold higher in C and CN lines than in N, but all accumulated ~150‐fold higher squalene due to substrate or storage limitations. Abnormal leaf and flower phenotypes, including lower pollen production and reduced fertility, were observed regardless of the compartment or level of transgene expression. Substantial changes in metabolomes of all lines were observed: levels of 65–120 unrelated metabolites, including the toxic alkaloid nicotine, changed by as much as 32‐fold. Profound effects of transgenesis on nontarget gene expression included changes in the abundance of 19 076 transcripts by up to 2000‐fold in CN; 7784 transcripts by up to 1400‐fold in N; and 5224 transcripts by as much as 2200‐fold in C. Transporter‐related transcripts were induced, and cell cycle‐associated transcripts were disproportionally repressed in all three lines. Transcriptome changes were validated by qRT‐PCR. The mechanism underlying these large changes likely involves metabolite‐mediated anterograde and/or retrograde signalling irrespective of the level of transgene expression or end product, due to imbalance of metabolic pools, offering new insight into both anticipated and unanticipated consequences of metabolic engineering.