PLOS ONE: Karakterisering av differentiellt uttryckta gener som är involverade i Pathways associerad med magcancer

Abstrakt

För att undersöka mönster av genuttryck i magcancer, totalt 26 parade magcancer och noncancerous vävnader från patienter rekryterades för genuttryck microarray analyser. Limma metoder användes för att analysera data och gener ansågs vara signifikant differentiellt uttryckt om False Discovery Rate (FDR) värde var < 0,01, P
-värdet var < 0,01 och fold change (FC) var > 2. Därefter tillsattes Gene ontologi (GO) kategorier som används för att analysera de viktigaste funktionerna hos de differentiellt uttryckta gener. Enligt Kyoto Encyclopedia of gener och genom (Kegg) databank har vi hittat vägar signifikant samband med differential generna. Gene-Act nätverk och samexpression nätverk byggdes respektive baserade på relationerna mellan de gener, proteiner och föreningar i databasen. 2371 mRNA och 350 lncRNAs anses väsentligt differentiellt uttryckta gener valdes för ytterligare analys. GO kategorier, väg analyserar och Gene-Act nätverk visade en konsekvent resultat som uppregleras gener var ansvariga för tumörbildning, migration, angiogenes och mikrobildning, medan nedregleras gener involverade i metabolismen. Dessa resultat av denna studie ge några nya rön om kodande RNA, lncRNAs, vägar och samuttryck nätverk i magcancer som kommer att vara användbart för att styra ytterligare utredning och rikta terapi för denna sjukdom

Citation. Li H Yu B, Li J, Su L, Yan M, Zhang J, et al. (2015) Karakterisering av differentiellt uttryckta gener som är involverade i Pathways associerad med magcancer. PLoS ONE 10 (4): e0125013. doi: 10.1371 /journal.pone.0125013

Academic Redaktör: Francisco J. Esteban, University of Jaen, Spanien

Mottagna: 9 november 2014. Accepteras: 6 mars 2015, Publicerad: 30 April, 2015

Copyright: © 2015 Li et al. Detta är en öppen tillgång artikel distribueras enligt villkoren i Creative Commons Attribution License, som tillåter obegränsad användning, distribution och reproduktion i alla medier, förutsatt den ursprungliga författaren och källan kredit

datatillgänglighet: Alla relevanta uppgifter är inom pappers- och dess stödinformationsfiler. Alla microarray-filer är tillgängliga från NCBI Gene Expression Omnibus (GEO) databas (nummer "GSE65801", http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE65801).

Finansiering: Detta arbete har finansierats med bidrag för analys från National Natural Science Foundation i Kina [No. 81172324, No. 91229106, nr 81272749 och nr 81372231], vetenskap och teknik kommissionen i Shanghai kommun [No. 13ZR1425600], och projekt i National Science & Teknik pelare Program Kina (nr 2014BAI09B03). Finansiärerna hade ingen roll i studiedesign, datainsamling och analys, beslut att publicera, eller beredning av manuskriptet

Konkurrerande intressen:.. Författarna har förklarat att inga konkurrerande intressen finns

Introduktion

Gastric cancer (GC) är en av de vanligaste cancer i världen, och dess förekomst är särskilt hög i östra Asien, särskilt i Kina. Cirka 952 tusen nya fall av magcancer diagnostiserades över hela världen under 2012, och hälften av dem inträffade i östra Asien (främst i Kina) [1]. I Kina, är majoriteten av patienterna med GC diagnosen i ett sent skede med dålig prognos. Därför belysa de molekylära mekanismerna bakom GC progression är viktigt att identifiera viktiga biomarkörer och utveckling av effektiva riktade terapier.

Under det senaste årtiondet har genuttryck microarrays blivit ett vanligt verktyg för att undersöka gen-transkriptnivåer inom cancerforskningen. Microarray data används för en mängd olika analyser, såsom oövervakad klustring, klassificering, differentiellt uttryck analys och uttryck kartläggning av kvantitativa drag loci [2]. Det är inte bara hjälper till att identifiera viktiga dysfunktionella gener i cancer men ger genomet hela information om genuttryck samtidigt samt [3,4]. I denna studie, genomförde vi en genomet hela undersökning av uttrycket av lncRNAs och mRNA från parade prover av primära magcancer vävnader och noncancerous vävnader, att profilera differentiellt uttryckta lncRNAs och kodande RNA. Studie av dessa uppgifter kommer att ge värdefull information om mekanismen för cancer och tillåta upptäckten av viktiga gener som kan fungera som framtida mål för cancerbehandling.

Metoder och material

Etiska uttalande

Skrivet informerat samtycke erhölls från alla deltagare. Studien godkändes av den mänskliga forsknings etikkommitté Ruijin Hospital, Shanghai Jiao Tong University, School of Medicine.

Vävnadsprov

Vävnader togs från primärgastriska karcinom från obehandlade patienter som genomgick D2 radikal gastrektomi i Shanghai Ruijin Hospital. För varje cancervävnad, var ett parat noncancerous vävnadsprov från det intilliggande område på samma gång. Storleken på varje prov var omkring 0,1 cm ^{3. Alla prover placerades i RNAlater inom 15 minuter efter excision och lagrades i flytande kväve tills RNA-extraktion. I denna studie var 32 parade vävnader samlas in för microarray och 26 parade prover var inskrivna för analysen nästa steg i GO, väg och nätverket efter kvalitetskontroll med hjälp av 3D Principal komponentanalys (3D-PCA) och klusteranalys.}

Microarray experiment

Agilent SurePrint G3 Human GE 8x60K Microarray (Design ID: 028004) användes i denna studie. Totalt RNA isolerades och förstärks med hjälp av en låg ingång Quick Amp Labeling Kit, en färg (Catȇ0-2305, Agilent Technologies, USA). Därefter tillsattes de märkta CRNAs renas genom en RNeasy mini kit (Cat˥06, QIAGEN, Tyskland).

Baserat på tillverkarens anvisningar, ades varje objektglas hybridiseras med 600ng Cy3-märkt cRNA med användning av en Gene Expression Hybridization Kit (CatȆ8-5242, Agilent Technologies, USA) och tvättas med Gene Expression Wash Buffer Kit (CatȆ8-5327, Agilent Technologies, USA).

En Agilent microarray scanner (Cat # G2565CA, Agilent teknik, USA) och Feature Extraction programvara 10,7 (Agilent Technologies, USA) anbringades på scanna varje bild med samma inställningar som visas som följer, Dye kanal: Grön, Skanningsupplösning upplösning~~POS=HEADCOMP = 3 pm, 20bit. Rådata normaliserades av kvantiluppskattaren algoritmen, Gene Spring Software 11,0 (Agilent Technologies, USA) (beskrivs i S5 tabell).

limma

Linjära modeller och empiriska Bayes metoder tillämpas för att analysera data i denna studie. De resulterande P
-värden justerades med hjälp av BH FDR algoritmen. Det fanns tre standarder för oss att anse att en gen var signifikant differentiellt uttryckta, var FDR värde < 0,01, P
-värdet var < 0,01 och faldig förändring var > 2. (Beskrivs i S5 tabell) katalog

GO kategori

Vi utförde Gene ontologi (GO) analyser för att analysera funktionerna hos de differentiellt uttryckta generna i vår microarray enligt den fördelningsnyckel funktionella klassificeringen av The National Center for Biotechnology Information (NCBI). Generellt Fishers exakta test och χ
^{2 prov applicerades på klassificera GO kategori och falska upptäckten hastighet (FDR,) beräknades för att korrigera P
-värde ( N
_{k
avser antalet Fishers test P
-värden mindre än χ
2 prov P
-värden). Anrikningen Re gavs av: Re = ( n
f
/ N
) /( N
f
/ N
) i de betydande kategorier ( N
f
är antalet differential gener inom viss kategori, n
är det totala antalet gener inom samma kategori, n
f
är antalet differential gener i hela microarray, och N
är det totala antalet gener i mikromatrisen) (beskrivs i detalj i S5 tabell) Review}}

Pathway analyser

Pathway annoteringar hos differential exressed generna erhölls från Kegg (http: //www .genome.jp /Kegg /). Pathway kategorier med en FDR < 0,01 markerades. Anrikningen av betydande vägar gavs av: anrikning = /, som hjälpt oss att hitta mer betydelsefulla vägar i vår studie ( n
_{g
är antalet differential gener inom särskilt väg, n
en
är det totala antalet gener i samma väg, N
g
är antal differential gener som har åtminstone en väg anteckning, och N
en
är antalet gener som har åtminstone en väg anteckning i hela microarray.) (detaljerade i S5 tabell).}

Gene-Act nätverk

Enligt Kegg databasen, kan en gen involverad i flera vägar eller interagera med flera andra gener. Alla gen-gen interaktion slogs samman för att bygga Gene-lagen nät baserat på differential vägar, som hjälpte oss att avslöja de signalvägar och viktiga reglerande gener i GC.

Samuttryck nätverks

Gene samuttryck nätverk byggdes enligt den normaliserade signalintensiteten av specifika uttryck gener. Grad centra definieras som antalet länkar en nod har till en annan, som bestämmer den relativa betydelsen av gener. Vad mer, var k-kärnor används som en metod för att förenkla diagrammet topologi analyser. Kärn reglerande faktorer (gener) som har de högsta grader ansluta mest närliggande gener och bygga upp strukturen i nätverket (i detalj i S5 tabell).

realtid kvantitativ PCR

Totalt RNA extraherades från vävnader med hjälp av Trizol-reagens (Invitrogen) i enlighet med tillverkarens instruktioner. Den kvantitativa realtids-polymeraskedjereaktion (PCR) utfördes genom användning av SYBR-grön PCR Master Mix i en Snabb Realtids-PCR 7500 System (Applied Biosystems). Primrarna av de 10 generna visade i S4 Tabell. PCR-reaktioner utfördes vid 50 ° C under 2 min, följt av 40 cykler av 95 ° C under 15 s och 60 ° C under 1 min. ΔCt beräknades genom att subtrahera Ct av β-aktin-RNA (kontroll) från Ct av RNA av provet, respektive. ΔΔCt beräknades sedan genom att subtrahera ΔCt av styrningen från ΔCt av provet. Faldig förändring beräknades genom ekvationen 2-ΔΔCt.

Statistisk analys

SPSS 19 och Microsoft Excel 2010 användes för att analysera data. Uttrycksnivåer mellan cancervävnader och angränsande noncancerous vävnader analyserades med parade prov t-test. P
-värden under 0,05 betraktades som statistiskt signifikant.

Resultat

microarray analyser

Totalt 42,405 mänskliga gener profileras i vår studie av med användning av en Agilent G3 Human GE 8x60K microarray. Vi har lämnat vårt dataset i förvaret av "Gene Expression Omnibus" och numret var "GSE65801" (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE65801). Vi använde linjära modeller och empiriska Bayes metoder för att analysera data (se Metoder). Det fanns 2371 mRNA och 350 lncRNAs betraktas som de differentiellt uttryckta gener genom limma för analysen nästa steg (Fig 1A).

Bland alla 2371 differential mRNA, finns det 1142 mRNA nedregleras och 1229 mRNA upp- regleras i vår observation på förändringar av genuttryck mellan magcancer och kontrollvävnader (fig 1c). De flesta av de differentiella mRNA har visat sig vara korrelerade med carcinogenes och metastas i de flesta typer av cancer (tabell 1). Generna som GKN2, PGC, MUC6, Chia, PSCA och FBP2 var bland de 20 nedregleras gener, medan KLK8, SFRP4, INHBA, CLDN1, CST1, FAP, SPP1, OLFM4 och KRT17 var bland de 20 upp -regulated gener (tabell 1). Dock har vissa gener som HOXC9, FNDC1, STRA6, KCNE2, PGA3 och KCNJ16 inte rapporterats i magcancer och deras roller förblir okänd (tabell 1).

Dessutom fann vi 193 nedregleras lncRNAs och 156 upp-reglerade lncRNAs bland totalt 350 differential lncRNAs baserade på profilering (bild 1B). De flesta av de lncRNAs har inte fått en officiell namn och deras funktioner är fortfarande okända. Men några har rapporterats spela kritiska roller i cancer, såsom H19, GUCY1B2, MEG3 och AKR7L (tabell 2).

I vår tidigare rapport [36], den fold change (FC) i H19 i 74 magcancer kontra parade noncancerous vävnader var 6,015, med en P
-värde av 0,017. Detta resultat överensstämde med data från H19 (Absolute FC = 6.06) i detta microarray analyser. Dessutom överuttryck av H19 bidrar till spridning, migration, invasion och metastas av magcancer.

Gene Ontology kategorier

Alla differentiellt uttryckta generna klassificeras i olika funktionella kategorier beroende på Gene ontologi (GO) projekt för biologiska processer. Baserat på våra microarray uppgifter, analyser GO indikerade att 208 GO termer anrikades ( P Hotel < 0,01, FDR < 0,01) (S1 tabell). De primära GO kategorier för 170 upp-reglerade GO termer var fokuserade på celladhesion, angiogenes, flercellig organism utveckling, Axon vägledning, skelett systemutveckling, kollagen fibriller organisation, positiv reglering av angiogenes, såra och negativ reglering av celltillväxt (fig 2A) . De viktigaste GO kategorier för nedregleras gener var matsmältning, xenobiotiska metabolisk process, trans transport, jontransport, liten molekyl metabolisk process, negativ reglering av tillväxt, glutation metabolisk process, cellulärt svar för kadmium jon och metabolisk process (Fig 2B).

Enligt differential gener och funktioner, byggde vi en GO Tree att utforska samspelet mellan alla olika GO kategorier de. Mångfalden i dessa kategorier när man jämför cancer och kontrollvävnader föreslog att magcancer kan associeras med signifikant uppreglerat cellmigration, celltillväxt, angiogenes, cell-cell adhesion och cellytereceptor signalvägar, medan cellernas ämnesomsättning processer och jon trans transport är nedreglerade (fig 3).

Pathway analyser

Pathway analyser användes för att identifiera de betydande vägar som är förknippade med de differentiellt uttryckta generna enligt Kegg. Det fanns 32 upp-reglerade vägar och 31 nedregleras vägar baserade på våra data (Fig 4). Dessutom vägen profilering överensstämde med resultaten för GO kategorierna i cancerrelaterade biologiska funktioner. Våra data visade några differential gener högt upp-regleras som föreslog sina inblandade vägar var activiated. Till exempel var SFRP4, WNT11, FZD2, MYC starkt uttryckt i cancervävnader som representerar Wnt banan var activiated och BCL2A1, ICM1, TNFSF14 i NF-kB-vägen var starkt uttryckt liksom. De flesta av de cancerrelaterade signalvägar som JAK /STAT, Wnt, NF-kB, PI3K, mTOR, igelkott och Notch vägar aktiverades i magsäckscancer jämfört med noncancerous vävnader baserat på våra data (S2 tabell). De uppreglerat vägar som var inriktade på celladhesion, transkriptions dysreglering, cancer och differentieringen var korrelerade med tumorogenesis och metastaser (fig 4A). Men nedregleras vägar var i allmänhet ansvarig för metabolism (Fig 4B).

Gene-Act nätverk

Baserat på GO kategorier och väg analys kan en gen involverad i flera vägar eller interagera med flera andra gener. Vi poolade differential generna och byggt upp ett nätverk av interaktioner mellan differentiellt uttryckta gener. En hög grad proteinet reglerar eller regleras av många andra proteiner, vilket innebär en viktig roll i Gene-Act nätverket (S3 Table). Glutation S-transferas (GST) familj, cytokrom P450 (CYP) familj, UDP glukuronosyltransferas 2 (UGT2) familj, epidermal tillväxtfaktorreceptor (EGFR) familj och cAMP-beroende proteinkinas katalytisk beta (PRKACB) var i centrum för den gen-gen interaktion nätverk. De kan spela en nyckelroll i nätverket eftersom de hade den starkaste graden (grad > 25) centralities (gen-gen interaktioner) (Fig 5). Det har rapporterats att GST, EGFR och PRKACB är ansvariga för signaltransduktionsvägar involverade i tumörtillväxt och differentiering i olika typer av cancrar [42,43].

Gene samuttryck nätverks

vi producerade en gen samuttryck nät baserat på det differentiellt uttryckta gener, proteiner och proteinkomplex i cancervävnader och noncancerous (kontroll) vävnader, respektive. Jämfört med kontrollen anslutningarna mellan gener i cancervävnader var mindre, vilket tyder på att de flesta av de fysiologiska gen-gen interaktion och kopplingar i normala vävnader hade brutits eller förlorade i cancervävnader (Fig 6A och 6B). De gener med hög grad och k-kärna vilket betyder att de hade de flesta av interaktioner med andra geneswere kallas nyckelgener i samspelet nätverk (Fig 6B) inklusive TRO, GPR124, TIMP2, EMCN, SLIT3, HTRA1, SPARC, LAMA4 och MEOX2 (tabell 3). De var ansvariga för cellsignalering, adhesion, angiogenes, migration, tillväxt och metastaser.

Att microarray resultat genom qPCR

Vi utförde kvantitativ realtids-PCR (qPCR) den 6 uppregleras gener (COL1A, BGN, SPP1, Melk, IGFBP4, SPARC) och 4 nedregleras gener (PGC, SST, MT1X, S100P) för att verifiera våra data i magcancer vävnader (tumör) och noncancerous vävnader (Normal). Uttrycksförhållandena för dessa 10 gener (Tumör /Normal) från qPCR överensstämmer med dem från microarray (S4 tabell). Den föreslog data från differential gener uttryck från microarray var tillförlitliga. Dessutom har vårt team har arbetat med några av differential gener som PHF10 [55], CEACAM6 [56], SFRP1 [57], SOX11 [58], CLDN1 [59] för att undersöka deras uttryck och funktioner i magcancer och resultaten perfekta visat vår microarray data.

Diskussion

microarray genuttryck analyser av magcancer har tidigare använts för att förutsäga diagnostiska markörer [60] och för att identifiera gener uttrycksmönster i samband med prognosen [ ,,,0],61,62], men det har inte använts för att avslöja molekylära interaktioner mellan lncRNAs och mRNA i GC. I denna studie analyserade vi 26 gastriska cancervävnader med parade noncancerous vävnader och profilerad generna differentiellt uttryckt enligt deras GO kategorier, vägar, Gene-Act-nätverk och Co-Expression nätverk.

genuttryck resultat erhölls genom att använda en Agilent G3 Human GE 8x60K microarray, som inte bara omfattar transkriptom databaser för mRNA-mål men även prober för lncRNAs (långa icke-kodande RNA). Med kombinationen av mRNA och lncRNAs, kan den utföra två experiment på en enda mikromatris och förutsäga lncRNA funktion och interaktion med mRNA. Analyserna visade en uppsättning av gener som differentiellt uttryckta mellan magcancer och normal vävnad. Några av dem har tidigare redovisats i magsår eller andra cancerformer. Till exempel, uttryck av gastrokine-2 (GKN2) var betydligt nedregleras eller frånvarande i magcancer cellinjer, gastric intestinal metaplasi och tumörvävnad. Överuttryck av GKN2 bidragit till celltillväxt, migration och invasion av magcancer och grep cellcykeln vid G1-S övergångsfasen [6]. Däremot var nivåerna av uttryck av inhibin beta A (INHBA) signifikant högre i cancervävnad än i angränsande normal slemhinna, och det anses vara en oberoende prognostisk faktor i magcancer [22]. Dessutom upptäckte vi några nya gener, såsom TMEM184A, PSAPL1, KIAA1199, CLRN3 och FNDC1, som inte har rapporterats i magcancer tidigare, och deras roll i cancer är fortfarande okända.

En av fördelarna med vår genuttryck microarray analys är att det representerade ett uttryck för lncRNAs och mRNA så att båda kunde undersökas tillsammans. Vår tidigare rapport om den roll som lncRNA H19 och dess nätverk i GC [36] baserades på denna microarray data. Men de flesta av de lncRNAs såsom DRD5, FMO6P, SNAR-A3 och TPRXL visade i vår microarray inte har identifierats och behöver ytterligare utredning för att klargöra sina roller i magcancer.

Baserat på vår genuttryck profilering data, generna och deras funktioner aktiveras i magcancer var ansvariga för spridning, adhesion, migration och metastasering, vilket var i linje med resultaten från vägen analyser. Intressant nog upptäckte vi att de flesta av de cancerrelaterade signalvägar som rapporterats tidigare, såsom Notch, mTOR och Hedgehog aktiverades i GC baserat på våra data. Dessa resultat stöder synpunkten att heterogenitet är kännetecknande för GC. Jämförelse av samexpression nätverk mellan normala vävnader och cancer tyder på att uttrycket, funktioner och interaktioner hos majoriteten av fysiologiska genen förloras eller skadas i magcancer, medan proliferation, migration och metastaser var onormalt förbättras. Dessa intressanta fynd matcha egenskaperna hos cancer, såsom anaplasi och dedifferentiering. Dessa differentiellt uttryckta gener som är involverade i signalvägar fungerade som nyckelgener i samexpression nätverk kan vara potentiella mål för anti-cancerterapi eller diagnostiska markörer i framtiden.

Bakgrundsinformation
S1 tabell. Pathway analyser av differential gener
doi: 10.1371 /journal.pone.0125013.s001
(XLSX) Review S2 tabell. GO-analyser av differential gener
doi: 10,1371 /journal.pone.0125013.s002
(XLSX) Review S3 Table. Gene-Act nätverk av differential gener
doi: 10.1371 /journal.pone.0125013.s003
(XLSX) Review S4 Tabell. Primers och verifiering
doi: 10.1371 /journal.pone.0125013.s004
(DOCX) Review S5 tabell. Metoder och material
doi: 10.1371 /journal.pone.0125013.s005
(DOCX) Review

Tack till

Vi vill tacka Dr Fred Bogott på Austin Medical Center , University of Minnesota, och Dr. Joshua Liao vid Hormel institutet, Austin Minnesota, för deras engelska redigering av detta manuskript. Vi skulle vilja tacka finansiärer att stödja denna studie. Detta arbete har finansierats med bidrag från National Natural Science Foundation i Kina [No. 81172324, No. 91229106, nr 81272749 och nr 81372231], vetenskap och teknik kommissionen i Shanghai kommun [No. 13ZR1425600], och projekt i National Science & Teknik pelare Program Kina (nr 2014BAI09B03).

• PLOS ONE: Survival analys av patienter med intervall cancerpatienter som genomgår Gastric Cancer Screening av Endoscopy
• PLOS ONE: Prognostic Betydelsen av Matrix metalloproteinas-7 i Gastric Cancer Survival: A Meta-Analysis