Transcript shortening via alternative polyadenylation promotes gene expression during fracture healing

Abstract

Maturation of the 3′ end of almost all eukaryotic messenger RNAs (mRNAs) requires cleavage and polyadenylation. Most mammalian mRNAs are polyadenylated at different sites within the last exon, generating alternative polyadenylation (APA) isoforms that have the same coding region but distinct 3′ untranslated regions (UTRs). The 3′UTR contains motifs that regulate mRNA metabolism; thus, changing the 3′UTR length via APA can significantly affect gene expression. Endochondral ossification is a central process in bone healing, but the impact of APA on gene expression during this process is unknown. Here, we report the widespread occurrence of APA, which impacts multiple pathways that are known to participate in bone healing. Importantly, the progression of endochondral ossification involves global 3′UTR shortening, which is coupled with an increased abundance of shortened transcripts relative to other transcripts; these results highlight the role of APA in promoting gene expression during endochondral bone formation. Our mechanistic studies of transcripts that undergo APA in the fracture callus revealed an intricate regulatory network in which APA enhances the expression of the collagen, type I, alpha 1 (Col1a1) and Col1a2 genes, which encode the 2 subunits of the abundantly expressed protein collagen 1. APA exerts this effect by shortening the 3′UTRs of the Col1a1 and Col1a2 mRNAs, thus removing the binding sites of miR-29a-3p, which would otherwise strongly promote the degradation of both transcripts. Taken together, our study is the first to characterize the crucial roles of APA in regulating the 3′UTR landscape and modulating gene expression during fracture healing.

References

[1] Claes L, Recknagel S, Ignatius A. Fracture healing under healthy and inflammatory conditions. Nat. Rev. Rheumatol., 2012, 8: 133-143

[2] Bahney CS et al. Cellular biology of fracture healing. J. Orthop. Res., 2019, 37: 35-50

[3] Einhorn TA, Gerstenfeld LC. Fracture healing: mechanisms and interventions. Nat. Rev. Rheumatol., 2015, 11: 45-54

[4] Richard P, Manley JL. Transcription termination by nuclear RNA polymerases. Genes Dev., 2009, 23: 1247-1269

[5] Shi Y, Manley JL. The end of the message: multiple protein-RNA interactions define the mRNA polyadenylation site. Genes Dev., 2015, 29: 889-897

[6] Proudfoot NJ. Transcriptional termination in mammals: stopping the RNA polymerase II juggernaut. Science, 2016, 352: aad9926

[7] Tian B, Manley JL. Alternative polyadenylation of mRNA precursors. Nat. Rev. Mol. Cell Biol., 2017, 18: 18-30

[8] Gruber AJ, Zavolan M. Alternative cleavage and polyadenylation in health and disease. Nat. Rev. Genet., 2019, 20: 599-614

[9] Derti A et al. A quantitative atlas of polyadenylation in five mammals. Genome Res., 2012, 22: 1173-1183

[10] Mayr C. Evolution and biological roles of alternative 3’UTRs. Trends Cell Biol., 2016, 26: 227-237

[11] Bartel DP. Metazoan microRNAs. Cell, 2018, 173: 20-51

[12] Ha M, Kim VN. Regulation of microRNA biogenesis. Nat. Rev. Mol. Cell Biol., 2014, 15: 509-524

[13] Tian B, Pan Z, Lee JY. Widespread mRNA polyadenylation events in introns indicate dynamic interplay between polyadenylation and splicing. Genome Res., 2007, 17: 156-165

[14] Singh I et al. Widespread intronic polyadenylation diversifies immune cell transcriptomes. Nat. Commun., 2018, 9

[15] Xia Z et al. Dynamic analyses of alternative polyadenylation from RNA-seq reveal a 3’-UTR landscape across seven tumour types. Nat. Commun., 2014, 5

[16] The Cancer Genome Atlas Research Network, et al The Cancer Genome Atlas Pan-Cancer analysis project. Nat. Genet., 2013, 45: 1113-1120

[17] Sheppard S, Lawson ND, Zhu LJ. Accurate identification of polyadenylation sites from 3’ end deep sequencing using a naive Bayes classifier. Bioinformatics, 2013, 29: 2564-2571

[18] Zhang Y et al. Alternative polyadenylation: methods, mechanism, function, and role in cancer. J. Exp. Clin. Cancer Res., 2021, 40: 51

[19] Wang R, Tian B. APAlyzer: a bioinformatics package for analysis of alternative polyadenylation isoforms. Bioinformatics, 2020, 36: 3907-3909

[20] Wang R, Nambiar R, Zheng D, Tian B. PolyA_DB 3 catalogs cleavage and polyadenylation sites identified by deep sequencing in multiple genomes. Nucleic Acids Res., 2018, 46: D315-D319

[21] Wang R, Zheng D, Yehia G, Tian B. A compendium of conserved cleavage and polyadenylation events in mammalian genes. Genome Res., 2018, 28: 1427-1441

[22] Nourse, J., Spada, S. & Danckwardt, S. Emerging roles of RNA 3’-end cleavage and polyadenylation in pathogenesis, diagnosis and therapy of human disorders. Biomolecules 10 915 (2020).

[23] Khajuria DK et al. Aberrant structure of fibrillar collagen and elevated levels of advanced glycation end products typify delayed fracture healing in the diet-induced obesity mouse model. Bone, 2020, 137: 115436

[24] Brown ML et al. Delayed fracture healing and increased callus adiposity in a C57BL/6J murine model of obesity-associated type 2 diabetes mellitus. PLoS One, 2014, 9: e99656

[25] Kaback LA et al. Osterix/Sp7 regulates mesenchymal stem cell mediated endochondral ossification. J. Cell Physiol., 2008, 214: 173-182

[26] Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol., 2014, 15

[27] Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinforma., 2009, 10

[28] Ge SX, Son EW, Yao R. iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data. BMC Bioinforma., 2018, 19

[29] Jiao X et al. DAVID-WS: a stateful web service to facilitate gene/protein list analysis. Bioinformatics, 2012, 28: 1805-1806

[30] San Antonio, J. D., Jacenko, O., Fertala, A. & Orgel, J. Collagen structure-function mapping informs applications for regenerative medicine. Bioengineering 8, 3 (2020).

[31] McGeary, S. E. et al. The biochemical basis of microRNA targeting efficacy. Science 366, eaav1741 (2019).

[32] Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res., 2020, 48: D127-D131

[33] van Rooij E et al. Dysregulation of microRNAs after myocardial infarction reveals a role of miR-29 in cardiac fibrosis. Proc. Natl. Acad. Sci. USA, 2008, 105: 13027-13032

[34] Dey, S. et al. Loss of miR-29a/b1 promotes inflammation and fibrosis in acute pancreatitis. JCI Insight 6, e149539 (2021).

[35] Li Z et al. Biological functions of miR-29b contribute to positive regulation of osteoblast differentiation. J. Biol. Chem., 2009, 284: 15676-15684

[36] Ulrich V et al. Chronic miR-29 antagonism promotes favorable plaque remodeling in atherosclerotic mice. EMBO Mol. Med., 2016, 8: 643-653

[37] Kramer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics, 2014, 30: 523-530

[38] Aghajanian P, Mohan S. The art of building bone: emerging role of chondrocyte-to-osteoblast transdifferentiation in endochondral ossification. Bone Res., 2018, 6: 19

[39] Lianoglou S, Garg V, Yang JL, Leslie CS, Mayr C. Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression. Genes Dev., 2013, 27: 2380-2396

[40] Gruber AJ et al. Discovery of physiological and cancer-related regulators of 3’ UTR processing with KAPAC. Genome Biol., 2018, 19

[41] Gruber AJ, Gypas F, Riba A, Schmidt R, Zavolan M. Terminal exon characterization with TECtool reveals an abundance of cell-specific isoforms. Nat. Methods, 2018, 15: 832-836

[42] Zhou, H. et al. IRAK2 directs stimulus-dependent nuclear export of inflammatory mRNAs. Elife 6, e29630 (2017).

[43] Cushing L et al. miR-29 is a major regulator of genes associated with pulmonary fibrosis. Am. J. Respir. Cell Mol. Biol., 2011, 45: 287-294

[44] Hanna J, Hossain GS, Kocerha J. The potential for microRNA therapeutics and clinical research. Front. Genet., 2019, 10: 478

[45] Le Bleu HK et al. Extraction of high-quality RNA from human articular cartilage. Anal. Biochem., 2017, 518: 134-138

[46] Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods, 2015, 12: 357-360

[47] Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics, 2015, 31: 166-169

[48] Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, 2010, 26: 139-140

[49] Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res., 2002, 30: 207-210

[50] Leclerc GJ, Leclerc GM, Barredo JC. Real-time RT-PCR analysis of mRNA decay: half-life of Beta-actin mRNA in human leukemia CCRF-CEM and Nalm-6 cell lines. Cancer Cell Int., 2002, 2

Funding

U.S. Department of Health & Human Services | NIH | National Institute of Diabetes and Digestive and Kidney Diseases (National Institute of Diabetes & Digestive & Kidney Diseases)(DK121327)