Expression of Linear and Novel Circular Forms of an -Associated Non-Coding RNA Correlates with Atherosclerosis Risk
Human genome-wide association studies have linked single nucleotide polymorphisms (SNPs) on chromosome 9p21.3 near the INK4/ARF (CDKN2a/b) locus with susceptibility to atherosclerotic vascular disease (ASVD). Although this locus encodes three well-characterized tumor suppressors, p16INK4a, p15INK4b, and ARF, the SNPs most strongly associated with ASVD are ∼120 kb from the nearest coding gene within a long non-coding RNA (ncRNA) known as ANRIL (CDKN2BAS). While individuals homozygous for the atherosclerotic risk allele show decreased expression of ANRIL and the coding INK4/ARF transcripts, the mechanism by which such distant genetic variants influence INK4/ARF expression is unknown. Here, using rapid amplification of cDNA ends (RACE) and analysis of next-generation RNA sequencing datasets, we determined the structure and abundance of multiple ANRIL species. Each of these species was present at very low copy numbers in primary and cultured cells; however, only the expression of ANRIL isoforms containing exons proximal to the INK4/ARF locus correlated with the ASVD risk alleles. Surprisingly, RACE also identified transcripts containing non-colinear ANRIL exonic sequences, whose expression also correlated with genotype and INK4/ARF expression. These non-polyadenylated RNAs resisted RNAse R digestion and could be PCR amplified using outward-facing primers, suggesting they represent circular RNA structures that could arise from by-products of mRNA splicing. Next-generation DNA sequencing and splice prediction algorithms identified polymorphisms within the ASVD risk interval that may regulate ANRIL splicing and circular ANRIL (cANRIL) production. These results identify novel circular RNA products emanating from the ANRIL locus and suggest causal variants at 9p21.3 regulate INK4/ARF expression and ASVD risk by modulating ANRIL expression and/or structure.
Vyšlo v časopise:
Expression of Linear and Novel Circular Forms of an -Associated Non-Coding RNA Correlates with Atherosclerosis Risk. PLoS Genet 6(12): e32767. doi:10.1371/journal.pgen.1001233
Kategorie:
Research Article
prolekare.web.journal.doi_sk:
https://doi.org/10.1371/journal.pgen.1001233
Souhrn
Human genome-wide association studies have linked single nucleotide polymorphisms (SNPs) on chromosome 9p21.3 near the INK4/ARF (CDKN2a/b) locus with susceptibility to atherosclerotic vascular disease (ASVD). Although this locus encodes three well-characterized tumor suppressors, p16INK4a, p15INK4b, and ARF, the SNPs most strongly associated with ASVD are ∼120 kb from the nearest coding gene within a long non-coding RNA (ncRNA) known as ANRIL (CDKN2BAS). While individuals homozygous for the atherosclerotic risk allele show decreased expression of ANRIL and the coding INK4/ARF transcripts, the mechanism by which such distant genetic variants influence INK4/ARF expression is unknown. Here, using rapid amplification of cDNA ends (RACE) and analysis of next-generation RNA sequencing datasets, we determined the structure and abundance of multiple ANRIL species. Each of these species was present at very low copy numbers in primary and cultured cells; however, only the expression of ANRIL isoforms containing exons proximal to the INK4/ARF locus correlated with the ASVD risk alleles. Surprisingly, RACE also identified transcripts containing non-colinear ANRIL exonic sequences, whose expression also correlated with genotype and INK4/ARF expression. These non-polyadenylated RNAs resisted RNAse R digestion and could be PCR amplified using outward-facing primers, suggesting they represent circular RNA structures that could arise from by-products of mRNA splicing. Next-generation DNA sequencing and splice prediction algorithms identified polymorphisms within the ASVD risk interval that may regulate ANRIL splicing and circular ANRIL (cANRIL) production. These results identify novel circular RNA products emanating from the ANRIL locus and suggest causal variants at 9p21.3 regulate INK4/ARF expression and ASVD risk by modulating ANRIL expression and/or structure.
Zdroje
1. World Health Organization W 2009 Cardiovascular diseases Fact Sheet
2. BirosE
CooperM
PalmerLJ
WalkerPJ
NormanPE
2010 Association of an allele on chromosome 9 and abdominal aortic aneurysm. Atherosclerosis
3. BownMJ
BraundPS
ThompsonJ
LondonNJ
SamaniNJ
2008 Association between the coronary artery disease risk locus on chromosome 9p21.3 and abdominal aortic aneurysm. Circ Cardiovasc Genet 1 39 42
4. ThompsonAR
GolledgeJ
CooperJA
HafezH
NormanPE
2009 Sequence variant on 9p21 is associated with the presence of abdominal aortic aneurysm disease but does not have an impact on aneurysmal expansion. Eur J Hum Genet 17 391 394
5. HelgadottirA
ThorleifssonG
ManolescuA
GretarsdottirS
BlondalT
2007 A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science 316 1491 1493
6. HelgadottirA
ThorleifssonG
MagnussonKP
GretarsdottirS
SteinthorsdottirV
2008 The same sequence variant on 9p21 associates with myocardial infarction, abdominal aortic aneurysm and intracranial aneurysm. Nat Genet 40 217 224
7. McPhersonR
PertsemlidisA
KavaslarN
StewartA
RobertsR
2007 A common allele on chromosome 9 associated with coronary heart disease. Science 316 1488 1491
8. MatarinM
BrownWM
SingletonA
HardyJA
MeschiaJF
2008 Whole genome analyses suggest ischemic stroke and heart disease share an association with polymorphisms on chromosome 9p21. Stroke 39 1586 1589
9. GschwendtnerA
BevanS
ColeJW
PlourdeA
MatarinM
2009 Sequence variants on chromosome 9p21.3 confer risk for atherosclerotic stroke. Ann Neurol 65 531 539
10. SmithJG
MelanderO
LovkvistH
HedbladB
EngstromG
2009 Common genetic variants on chromosome 9p21 confers risk of ischemic stroke: a large-scale genetic association study. Circ Cardiovasc Genet 2 159 164
11. CluettC
McDermottMM
GuralnikJ
FerrucciL
BandinelliS
2009 The 9p21 myocardial infarction risk allele increases risk of peripheral artery disease in older people. Circ Cardiovasc Genet 2 347 353
12. YeS
WilleitJ
KronenbergF
XuQ
KiechlS
2008 Association of genetic variation on chromosome 9p21 with susceptibility and progression of atherosclerosis: a population-based, prospective study. J Am Coll Cardiol 52 378 384
13. BroadbentHM
PedenJF
LorkowskiS
GoelA
OngenH
2008 Susceptibility to coronary artery disease and diabetes is encoded by distinct, tightly linked SNPs in the ANRIL locus on chromosome 9p. Hum Mol Genet 17 806 814
14. SeidelmannSB
KuoC
PleskacN
MolinaJ
SayersS
2008 Athsq1 is an atherosclerosis modifier locus with dramatic effects on lesion area and prominent accumulation of versican. Arterioscler Thromb Vasc Biol 28 2180 2186
15. KimWY
SharplessNE
2006 The regulation of INK4/ARF in cancer and aging. Cell 127 265 275
16. LiuY
SanoffHK
ChoH
BurdCE
TorriceC
2009 INK4/ARF transcript expression is associated with chromosome 9p21 variants linked to atherosclerosis. PLoS One 4 e5027
17. WesselyR
2010 Atherosclerosis and cell cycle: put the brakes on! Critical role for cyclin-dependent kinase inhibitors. J Am Coll Cardiol 55 2269 2271
18. BoehmM
NabelEG
2003 The cell cycle and cardiovascular diseases. Prog Cell Cycle Res 5 19 30
19. GizardF
AmantC
BarbierO
BellostaS
RobillardR
2005 PPAR alpha inhibits vascular smooth muscle cell proliferation underlying intimal hyperplasia by inducing the tumor suppressor p16INK4a. J Clin Invest 115 3228 3238
20. Gonzalez-NavarroH
Abu NabahYN
VinueA
Andres-ManzanoMJ
ColladoM
2010 p19(ARF) deficiency reduces macrophage and vascular smooth muscle cell apoptosis and aggravates atherosclerosis. J Am Coll Cardiol 55 2258 2268
21. KalininaN
AgrotisA
AntropovaY
IlyinskayaO
SmirnovV
2004 Smad expression in human atherosclerotic lesions: evidence for impaired TGF-beta/Smad signaling in smooth muscle cells of fibrofatty lesions. Arterioscler Thromb Vasc Biol 24 1391 1396
22. ReynisdottirI
PolyakK
IavaroneA
MassagueJ
1995 Kip/Cip and Ink4 Cdk inhibitors cooperate to induce cell cycle arrest in response to TGF-beta. Genes Dev 9 1831 1845
23. GraingerDJ
2004 Transforming growth factor beta and atherosclerosis: so far, so good for the protective cytokine hypothesis. Arterioscler Thromb Vasc Biol 24 399 404
24. JanzenV
ForkertR
FlemingHE
SaitoY
WaringMT
2006 Stem-cell ageing modified by the cyclin-dependent kinase inhibitor p16INK4a. Nature 443 421 426
25. Yvan-CharvetL
PaglerT
GautierEL
AvagyanS
SiryRL
2010 ATP-binding cassette transporters and HDL suppress hematopoietic stem cell proliferation. Science 328 1689 1693
26. ViselA
ZhuY
MayD
AfzalV
GongE
2010 Targeted deletion of the 9p21 non-coding coronary artery disease risk interval in mice. Nature 464 409 412
27. PasmantE
LaurendeauI
HeronD
VidaudM
VidaudD
2007 Characterization of a germ-line deletion, including the entire INK4/ARF locus, in a melanoma-neural system tumor family: identification of ANRIL, an antisense noncoding RNA whose expression coclusters with ARF. Cancer Res 67 3963 3969
28. FolkersenL
KyriakouT
GoelA
PedenJ
MalarstigA
2009 Relationship between CAD risk genotype in the chromosome 9p21 locus and gene expression. Identification of eight new ANRIL splice variants. PLoS One 4 e7677
29. KyriakouT
PalA
PedenJ
GreenF
GloynA
2009 ANRIL, The non-coding RNA present in the chromosome 9 CAD associated locus, has multiple splice variants and a potential regulatory role in CDKN2B expression. Atherosclerosis 207 e3
30. JarinovaO
StewartAF
RobertsR
WellsG
LauP
2009 Functional analysis of the chromosome 9p21.3 coronary artery disease risk locus. Arterioscler Thromb Vasc Biol 29 1671 1677
31. JacobsJJ
KieboomK
MarinoS
DePinhoRA
van LohuizenM
1999 The oncogene and Polycomb-group gene bmi-1 regulates cell proliferation and senescence through the ink4a locus. Nature 397 164 168
32. DhawanS
TschenSI
BhushanA
2009 Bmi-1 regulates the Ink4a/Arf locus to control pancreatic beta-cell proliferation. Genes Dev 23 906 911
33. ChenH
GuX
SuIH
BottinoR
ContrerasJL
2009 Polycomb protein Ezh2 regulates pancreatic beta-cell Ink4a/Arf expression and regeneration in diabetes mellitus. Genes Dev 23 975 985
34. KotakeY
CaoR
ViatourP
SageJ
ZhangY
2007 pRB family proteins are required for H3K27 trimethylation and Polycomb repression complexes binding to and silencing p16INK4alpha tumor suppressor gene. Genes Dev 21 49 54
35. BrackenAP
Kleine-KohlbrecherD
DietrichN
PasiniD
GargiuloG
2007 The Polycomb group proteins bind throughout the INK4A-ARF locus and are disassociated in senescent cells. Genes Dev 21 525 530
36. GilJ
PetersG
2006 Regulation of the INK4b-ARF-INK4a tumour suppressor locus: all for one or one for all. Nat Rev Mol Cell Biol 7 667 677
37. TerranovaR
YokobayashiS
StadlerMB
OtteAP
van LohuizenM
2008 Polycomb group proteins Ezh2 and Rnf2 direct genomic contraction and imprinted repression in early mouse embryos. Dev Cell 15 668 679
38. RinnJL
KerteszM
WangJK
SquazzoSL
XuX
2007 Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell 129 1311 1323
39. ZhaoJ
SunBK
ErwinJA
SongJJ
LeeJT
2008 Polycomb proteins targeted by a short repeat RNA to the mouse X chromosome. Science 322 750 756
40. PandeyRR
MondalT
MohammadF
EnrothS
RedrupL
2008 Kcnq1ot1 antisense noncoding RNA mediates lineage-specific transcriptional silencing through chromatin-level regulation. Mol Cell 32 232 246
41. GuptaRA
ShahN
WangKC
KimJ
HorlingsHM
2010 Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature 464 1071 1076
42. LiuY
SanoffHK
ChoH
BurdCE
TorriceC
2009 Expression of p16(INK4a) in peripheral blood T-cells is a biomarker of human aging. Aging Cell 8 439 448
43. HaraE
SmithR
ParryD
TaharaH
StoneS
1996 Regulation of p16CDKN2 expression and its implications for cell immortalization and senescence. Mol Cell Biol 16 859 867
44. KrishnamurthyJ
TorriceC
RamseyMR
KovalevGI
Al-RegaieyK
2004 Ink4a/Arf expression is a biomarker of aging. J Clin Invest 114 1299 1307
45. AuKF
JiangH
LinL
XingY
WongWH
2010 Detection of splice junctions from paired-end RNA-seq data by SpliceMap. Nucleic Acids Res
46. NielsenGP
Stemmer-RachamimovAO
ShawJ
RoyJE
KohJ
1999 Immunohistochemical survey of p16INK4A expression in normal human adult and infant tissues. Lab Invest 79 1137 1143
47. SchmidM
SenM
RosenbachMD
CarreraCJ
FriedmanH
2000 A methylthioadenosine phosphorylase (MTAP) fusion transcript identifies a new gene on chromosome 9p21 that is frequently deleted in cancer. Oncogene 19 5747 5754
48. ChenK
WallisJW
McLellanMD
LarsonDE
KalickiJM
2009 BreakDancer: an algorithm for high-resolution mapping of genomic structural variation. Nat Methods 6 677 681
49. ZaphiropoulosPG
1996 Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: correlation with exon skipping. Proc Natl Acad Sci U S A 93 6536 6541
50. SuzukiH
ZuoY
WangJ
ZhangMQ
MalhotraA
2006 Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing. Nucleic Acids Res 34 e63
51. DixonRJ
EperonIC
HallL
SamaniNJ
2005 A genome-wide survey demonstrates widespread non-linear mRNA in expressed sequences from multiple species. Nucleic Acids Res 33 5904 5913
52. LiH
DurbinR
2009 Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25 1754 1760
53. LiR
LiY
KristiansenK
WangJ
2008 SOAP: short oligonucleotide alignment program. Bioinformatics 24 713 714
54. WangZ
RolishME
YeoG
TungV
MawsonM
2004 Systematic identification and analysis of exonic splicing silencers. Cell 119 831 845
55. FairbrotherWG
YehRF
SharpPA
BurgeCB
2002 Predictive identification of exonic splicing enhancers in human genes. Science 297 1007 1013
56. WangZ
BurgeCB
2008 Splicing regulation: from a parts list of regulatory elements to an integrated splicing code. Rna 14 802 813
57. FrazerKA
BallingerDG
CoxDR
HindsDA
StuveLL
2007 A second generation human haplotype map of over 3.1 million SNPs. Nature 449 851 861
58. NigroJM
ChoKR
FearonER
KernSE
RuppertJM
1991 Scrambled exons. Cell 64 607 613
59. HoldtLM
BeutnerF
ScholzM
GielenS
GabelG
2010 ANRIL expression is associated with atherosclerosis risk at chromosome 9p21. Arterioscler Thromb Vasc Biol 30 620 627
60. CunningtonMS
Santibanez KorefM
MayosiBM
BurnJ
KeavneyB
2010 Chromosome 9p21 SNPs Associated with Multiple Disease Phenotypes Correlate with ANRIL Expression. PLoS Genet 6 e1000899
61. Messaoudi-AubertSE
NichollsJ
MaertensGN
BrookesS
BernsteinE
2010 Role for the MOV10 RNA helicase in polycomb-mediated repression of the INK4a tumor suppressor. Nat Struct Mol Biol 17 862 868
62. YapKL
LiS
Munoz-CabelloAM
RaguzS
ZengL
2010 Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a. Mol Cell 38 662 674
63. Mancini-DinardoD
SteeleSJ
LevorseJM
IngramRS
TilghmanSM
2006 Elongation of the Kcnq1ot1 transcript is required for genomic imprinting of neighboring genes. Genes Dev 20 1268 1282
64. KanduriC
ThakurN
PandeyRR
2006 The length of the transcript encoded from the Kcnq1ot1 antisense promoter determines the degree of silencing. Embo J 25 2096 2106
65. DixonRJ
EperonIC
SamaniNJ
2007 Complementary intron sequence motifs associated with human exon repetition: a role for intragenic, inter-transcript interactions in gene expression. Bioinformatics 23 150 155
66. CocquerelleC
DaubersiesP
MajerusMA
KerckaertJP
BailleulB
1992 Splicing with inverted order of exons occurs proximal to large introns. Embo J 11 1095 1098
67. BonenL
VogelJ
2001 The ins and outs of group II introns. Trends Genet 17 322 331
68. KosA
DijkemaR
ArnbergAC
van der MeidePH
SchellekensH
1986 The hepatitis delta (delta) virus possesses a circular RNA. Nature 323 558 560
69. TsagrisEM
Martinez de AlbaAE
GozmanovaM
KalantidisK
2008 Viroids. Cell Microbiol 10 2168 2179
70. ZegginiE
WeedonMN
LindgrenCM
FraylingTM
ElliottKS
2007 Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science 316 1336 1341
71. ScottLJ
MohlkeKL
BonnycastleLL
WillerCJ
LiY
2007 A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science 316 1341 1345
72. SheteS
HoskingFJ
RobertsonLB
DobbinsSE
SansonM
2009 Genome-wide association study identifies five susceptibility loci for glioma. Nat Genet 41 899 904
73. WrenschM
JenkinsRB
ChangJS
YehRF
XiaoY
2009 Variants in the CDKN2B and RTEL1 regions are associated with high-grade glioma susceptibility. Nat Genet 41 905 908
74. UnoS
ZembutsuH
HirasawaA
TakahashiA
KuboM
2010 A genome-wide association study identifies genetic variants in the CDKN2BAS locus associated with endometriosis in Japanese. Nat Genet 42 707 710
75. BeiJX
LiY
JiaWH
FengBJ
ZhouG
2010 A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci. Nat Genet 42 599 603
76. EmanueleE
FontanaJM
MinorettiP
GeroldiD
2010 Preliminary evidence of a genetic association between chromosome 9p21.3 and human longevity. Rejuvenation Res 13 23 26
77. SebastianiP
SolovieffN
PucaA
HartleySW
MelistaE
2010 Genetic Signatures of Exceptional Longevity in Humans. Science
78. HankerAB
MoritaS
RepaskyGA
RossDT
SeitzRS
2008 Tools to study the function of the Ras-related, estrogen-regulated growth inhibitor in breast cancer. Methods Enzymol 439 53 72
79. BrookesS
RoweJ
RuasM
LlanosS
ClarkPA
2002 INK4a-deficient human diploid fibroblasts are resistant to RAS-induced senescence. Embo J 21 2936 2945
80. ShieldsJM
ThomasNE
CreggerM
BergerAJ
LeslieM
2007 Lack of extracellular signal-regulated kinase mitogen-activated protein kinase signaling shows a new type of melanoma. Cancer Res 67 1502 1512
81. LiH
HandsakerB
WysokerA
FennellT
RuanJ
2009 The Sequence Alignment/Map format and SAMtools. Bioinformatics 25 2078 2079
Štítky
Genetika Reprodukčná medicínaČlánok vyšiel v časopise
PLOS Genetics
2010 Číslo 12
- Je „freeze-all“ pro všechny? Odborníci na fertilitu diskutovali na virtuálním summitu
- Gynekologové a odborníci na reprodukční medicínu se sejdou na prvním virtuálním summitu
Najčítanejšie v tomto čísle
- Functional Comparison of Innate Immune Signaling Pathways in Primates
- Expression of Linear and Novel Circular Forms of an -Associated Non-Coding RNA Correlates with Atherosclerosis Risk
- Genome-Wide Interrogation of Mammalian Stem Cell Fate Determinants by Nested Chromosome Deletions
- Histone H2A C-Terminus Regulates Chromatin Dynamics, Remodeling, and Histone H1 Binding