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
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa