Genominlaajuinen systemaattinen karakterisointi bZIP-transkriptiotekijöiden systemaattisesta karakterisoinnista ja niiden ilmentymisprofiileista siementen kehityksen aikana ja vasteena suolastressiin maapähkinässä

kesä 13, 2021
admin

BZIP-geenien tunnistaminen, fylogeneettinen analyysi ja ryhmäluokitus lajikkeessa A. duranensis ja A. ipaensis

Homologiahakujen ja domainien todentamisen perusteella A. duranensis- ja A. ipaensis -genomeista tunnistettiin yhteensä 50 ja 45 ainutlaatuista bZIP-geeniä. Yksityiskohtaiset tiedot näistä geeneistä, mukaan lukien geenin tunniste, genominen sijainti, domeenikoostumus ja ryhmäluokitus, on esitetty lisätiedostossa 1. Olemassa olevan nimikkeistöjärjestelmän mukaisesti annoimme jokaiselle näistä uusista bZIP-geeneistä yksilölliset nimet: AdbZIP1-50 ja AibZIP1-45. Kun bZIP-domeenit oli tarkistettu, 93 geenillä oli tyypillinen bZIP-domeeni, johon kuului muuttumaton N- × 7-R/K-motiivi perusalueella ja Leu-heptad-toisto, joka on sijoitettu täsmälleen yhdeksän aminohappoa R/K:sta ylävirtaan kohti C-päätettä (lisätiedosto 2). Kahdessa muussa bZIP-geenissä, AdbZIP28:ssa ja AibZIP22:ssa, oli epätavallinen substituutio perusalueella: konservoitunut Arg/Lys (R/K) oli korvattu IIe:llä (I). Tämä korvautuminen on raportoitu myös muilla lajeilla .

BZIP-geeniperheen systemaattinen tutkimus suoritettiin ensimmäisen kerran Arabidopsiksessa . Tässä analyysissä bZIP-geenien eri ryhmät erotettiin toisistaan ja nimettiin niiden fylogeneettisten suhteiden ja toiminnallisten eroavaisuuksien perusteella. Tämä luokittelujärjestelmä on sittemmin otettu käyttöön muille lajeille niiden omasta ja Arabidopsiksen genomista saatujen bZIP-geenien ryhmittelyn perusteella . Tässä Arachiksen ja Arabidopsiksen genomeista peräisin olevien bZIP-proteiinien maksimaalisen todennäköisyyden (ML) analyysin perusteella tunnistimme 11 erillistä bZIP-geeniklusteria (ryhmät A-I, S ja U), joilla kaikilla oli korkea bootstrap-tuki (kuva 1). Arachiksen bZIP-geenien alaryhmäluokitus vahvistettiin edelleen fylogeneettisellä puun rekonstruktiolla sen jälkeen, kun siihen oli lisätty soijapavun bZIP-geenit (lisätiedosto 3). Useimpiin bZIP-klaadeihin kuuluvat läheisesti sukua olevat Arachis bZIP:t ja niiden Arabidopsis-ortologit; klaadeissa E ja F ei ole vastaavia jäseniä A. duranensiksessa tai A. ipaensiksessa. Huomionarvoista on, että samaan kladiin kuuluvilla bZIP-geeneillä oli samankaltaisia ryhmäkohtaisia sekvenssiominaisuuksia, mukaan lukien eksoni/intronirakenne, intronivaiheet, MEME-motiivit ja sitoutumiskohdan rakenteen ennuste (analysoidaan tarkemmin jäljempänä). Tämä lajien välisen ryhmäklusteroinnin malli viittaa siihen, että ryhmäkohtaiset piirteet ovat syntyneet ennen Arachisin ja Arabidopsiksen erilaistumista. Eri kasvilajien bZIP-geeneihin on kuitenkin myös kertynyt useita eroja evoluution aikana.

Kuva 1
kuvio1

Pähkinän ja Arabidopsiksen bZIP-geenien fylogeneettinen analyysi. Eri lajien haarojen päissä olevat geenit on merkitty erivärisillä kolmioilla. Maapähkinän bZIP-proteiinit on ryhmitelty yhdeksään erilliseen kladiin (A-D, G-I, S ja U). bZIP-proteiinisekvenssit linjattiin ClustalX-ohjelmalla, ja fylogeneettinen puu rakennettiin PhyML-ohjelmassa maksimikelpoisuusmenetelmää käyttäen. Bootstrap-arvot perustuvat 100 toistoon

Arachiksen bZIP-geenien geenirakenne

Koska intronien ja eksonien järjestäytyminen saattaa viitata bZIP-geenien evolutiiviseen kehityskulkuun , tarkastelimme Arachiksen bZIP-geenien rakennetta, mukaan lukien intronien lukumäärää, pituutta ja splikointivaihetta (Additional file 4). Havaitsimme, että geenien kokonaisrakenteet olivat identtisiä tai samankaltaisia saman fylogeneettisen ryhmän Arachis bZIP-geeneissä. Kun tarkastellaan maapähkinän bZIP-geenien intronien lukumäärää, 24 prosenttia AdbZIP-geeneistä ja 22 prosenttia AibZIP-geeneistä oli intronittomia, ja niitä esiintyi yksinomaan ryhmissä S ja B. Introneja sisältävissä geeneissä intronien lukumäärä vaihteli AdbZIP- ja AibZIP-geeneissä yhdestä 13:een. Ryhmässä G olevissa bZIP-geeneissä oli eniten introneja, mikä on yhdenmukaista muiden palkokasvien genomeissa tehtyjen havaintojen kanssa.

Splikointivaiheet nimettiin kolmeksi splikointivaiheeksi: vaihe 0 (P0), splikointi tapahtui kodonin kolmannen nukleotidin jälkeen; vaihe 1 (P1), splikointi tapahtui kodonin ensimmäisen nukleotidin jälkeen; ja vaihe 2 (P2), splikointi tapahtui toisen nukleotidin jälkeen. Avoimien lukukehysten (ORF) sisällä olevien pilkkoutumiskohtien vaiheet olivat erilaisia, mutta ne olivat hyvin konservoituneita bZIP-domeenin perus- ja sarana-alueilla, koska kaikki muutokset näillä alueilla vaikuttaisivat niiden koodiin ja toimintaan. Intronin sijainnin ja bZIP-domeenissa olevien splikointivaiheiden esiintymisen tai lukumäärän perusteella Arachisin bZIP-geeneissä tunnistettiin neljä intronimallia (a-d) (kuva 2 ja lisätiedosto 2). Kuviossa a oli vain yksi introni, joka oli insertoitunut sarana-alueen -5-asemaan, Gln- ja Ala-aminohappojen väliin; tämä kuvio tunnistettiin kaikissa Arachiksen bZIP-geeneissä ryhmissä A ja G. Kuviossa b oli kaksi introni-insertiota vaiheessa 0, toinen perusalueella ja toinen sarana-alueella; tämä kuvio tunnistettiin kaikissa bZIP-geeneissä ryhmässä D. Tämä kuvio tunnistettiin kaikissa bZIP-geeneissä ryhmässä D. Kuviossa c oli yksi introni -20-asemassa perusalueella vaiheessa 2 (P2), ja se sisältää kaikki ryhmien C ja H bZIP-geenit. Kuviossa d ei ollut introneja perusalueella eikä sarana-alueella, ja se sisältää kaikki ryhmien B ja S bZIP-geenit. Lisäksi useimmat Arachiksen bZIP-geenit, joissa on kuvio d, olivat intronittomia lukuun ottamatta AdbZIP45:tä ja AibZIP40:tä. Kummassakin näistä geeneistä oli yksi introni perus- ja sarana-alueiden ulkopuolella. Tässä havaitut Arachiksen bZIP-domeenin spleikkausvaiheen mallit olivat yhdenmukaisia muissa lajeissa havaittujen mallien kanssa. Geenin rakenteen ja intronien vaiheiden suuri säilyvyys fylogeneettisten klastien sisällä tuki hyväksyttyä ryhmäluokitusta ja viittasi siihen, että näillä erilaisilla eksonien splikointikuvioilla voi olla tärkeä rooli toiminnallisessa evoluutiossa.

Kuva 2
kuvio2

Intronikuviot Arachiksen bZIP- domeenin emäs- ja saranavyöhykkeillä. Kuvan yläosassa näkyy bZIP-domeenin perusrakenne. P0 osoittaa, että intronin pilkkoutumiskohta on kodonien välissä, ja P2 osoittaa, että intronin pilkkoutumiskohta sijaitsee kodonin toisen ja kolmannen nukleotidin välissä. Intronin esiintyvyyden, intronin sijainnin ja pilkkoutumisvaiheen perusteella Arachisin bZIP-geeneissä oli neljä erilaista mallia (a-d). Yksityiskohtaiset tiedot maapähkinän bZIP-proteiinien bZIP-domeenin intronien sijainneista on esitetty lisätiedostossa 2

Arachiksen bZIP-geenien eri ryhmien motiivikoostumukset

M MEME-analyysityökalulla bZIP-geeneissä havaittiin bZIP-domeenin lisäksi monia muitakin konservoituneita motiiveja. Kuten kuvasta 3 käy ilmi, bZIP-domeenin ulkopuolisia konservoituneita motiiveja tunnistettiin yhteensä 18, ja kullekin alaryhmälle muodostettiin konsensusmotiivikoostumukset (lisätiedosto 5). Nämä konsensusmotiivit osoittivat, että motiivien kokonaiskoostumukset olivat samanlaisia saman alaryhmän sisällä mutta erilaisia eri ryhmien välillä. Tämä viittasi siihen, että bZIP-geenien toiminnallinen eroavuus voi määräytyä ryhmäkohtaisten motiivien perusteella. Näiden motiivien yksilöllinen tarkastelu osoitti, että monet niistä olivat ryhmäkohtaisia. Esimerkiksi motiivit 1, 2, 3 ja 10 tunnistettiin vain ryhmässä D, motiivit 5, 14 ja 15 tunnistettiin vain ryhmässä G, motiivi 6 tunnistettiin vain ryhmässä I ja motiivi 9 tunnistettiin vain ryhmässä H. Useat motiivit voivat liittyä tiettyihin biologisiin toimintoihin. Esimerkiksi motiivi 1 on DELAY OF GERMINATION (DOG) 1-domeeni, jota tarvitaan lepotilan indusointiin ja useisiin siementen kypsymisen osatekijöihin, osittain häiritsemällä ABA-signaalin välityskomponentteja. Motiivi 3 sisältää potentiaalisia kaseiinikinaasi II:n (CK II) fosforylaatiokohtia (S/TxxD/E), joilla on keskeinen rooli solujen jakautumisessa ja laajenemisessa ja jotka vaikuttavat erilaisiin kehitykseen ja stressiin reagoiviin reitteihin . Mielenkiintoista on, että nämä ryhmäkohtaiset motiivit on tunnistettu myös saman ryhmän bZIP:eissä muissa palkokasvien genomeissa , mikä viittaa siihen, että motiivikoostumus on konservoitunut kaikissa palkokasveissa.

Kuvio 3
kuvio3

Lisäkonservoitujen motiivien jakautuminen MEME:llä tunnistettuna. Kunkin maapähkinän bZIP-proteiiniryhmän motiivikoostumukset on esitetty bZIP-domeenin ja bZIP-domeenin ulkopuolisten konservoitujen lisämotiivien sijainnin perusteella. bZIP-domeenit on esitetty punaisella, kun taas muut motiivit on korostettu värillisillä laatikoilla, jotka on numeroitu 1-18. Ennustettujen konservoitujen motiivien yksityiskohdat on esitetty lisätiedostossa 6

Arachiksen bZIP:n DNA-sitoutumisalueen rakenne ja dimerisaatio-ominaisuudet

BZIP-domeenin ydinperusalue ja sarana-alue määrittävät toisistaan riippumatta DNA-sitoutumispesifisyyden, kuten useat kokeet ovat osoittaneet . Kahden muuttumattoman paikan, asparagiinin (Asn/N; asento: – 18) ja arginiinin (Arg/R; asento: – 10), epätavallinen korvaaminen muutti DNA:n sitoutumispesifisyyttä . Kohdistimme maapähkinän bZIP-proteiinien perus- ja sarana-alueiden aminohapposekvenssit konservoitujen ja polymorfisten aminohappojäännösten tunnistamiseksi kussakin ryhmässä (lisätiedosto 6). Yhdessäkään maapähkinän bZIP-proteiinissa ei havaittu Asn/N-korvauksia -18-asemassa. Kaikissa ryhmän I jäsenissä oli kuitenkin lysiiniä (Lys/K) arginiinin (R) sijasta – 10 -asemassa, mikä on yhdenmukaista muiden palkokasvilajien ryhmän I bZIP-proteiinien kanssa. Lisäksi AdbZIP28:ssa ja AibZIP22:ssa (ryhmä U) oli hydrofobinen isoleusiinijäännös (Ile/I) arginiinin (Arg/R) sijasta, ja tällaisen korvauksen osoitettiin estävän täysin bZIP:n affiniteetin AP1:n suhteen hiivassa, eikä se tunnista G-bokseja riisissä .

Leu-vetoketjusekvenssi välittää bZIP-proteiinien homo- ja/tai heterodimerisaatiota, ja näiden proteiinien tiedetään sitoutuvan DNA:han dimereina . Leu vetoketju alue koostuu heptadi toistoista, aminohapot viitataan a, b, c, d, e, f ja g kunkin heptadin sisällä . Koska a-, d-, e- ja g-asemissa olevat aminohapot ovat lähellä Leu-vetoketjun rajapintaa, nämä aminohapot määräävät ensisijaisesti Leu-vetoketjun oligomerisaation, dimerisaation vakauden ja dimerin spesifisyyden. Analysoimme maapähkinän bZIP-proteiinien a-, d-, e- ja g-asemissa esiintyvien aminohappojen koostumuksia (kuva 4a).

Kuva 4
kuva4

Arachiksen bZIP-proteiinien dimerisaatio-ominaisuuksien ennustaminen. a Piirakkadiagrammit, jotka osoittavat eri aminohappojen esiintymistiheyden kussakin Arachis bZIP-domeenien Leu-vetoketjun neljässä asennossa (a, d, e ja g). b Histogrammi Asn:n (N) esiintymistiheydestä Leu-vetoketjun a-asennossa kaikissa Arachis bZIP-proteiineissa. c Histogrammi, joka osoittaa vetovoimaisten tai hylkivien g↔e’ -parien esiintymistiheyden heptadia kohti kaikissa Arachis bZIP-proteiineissa. g↔e′-parit on luokiteltu neljään ryhmään g- ja e-asentojen sähköstaattisten varausten mukaan. Oranssilla laatikolla merkitty +/- vetovoima osoittaa, että g-asento on emäksinen ja sitä seuraava e-asento on hapan. Taivaansinisellä laatikolla merkitty -/+ vetovoima osoittaa, että g-asento on hapan ja sitä seuraava e-asento emäksinen. Emäksinen repulsiivinen (vaaleanpunainen laatikko) ja hapan repulsiivinen (vihreä laatikko) osoittavat, että g-asemalla ja sitä seuraavalla e-asemalla on samanlainen varaus, joko molemmilla emäksinen tai molemmilla hapan

A-asemassa, noin 20 % jäännöksistä oli asparagiinia (Asn/N), joka voi muodostaa polaarisen taskun hydrofobiseen rajapintaan, mikä mahdollistaa vakaammat N-N-vuorovaikutukset a↔a′:ssä (vastaava asema vastakkaisessa kierteessä) verrattuna muihin aminohappoihin . Eri heptadeissa toisessa ja viidennessä heptadissa oli eniten Asn/N-jäännöksiä a-asennossa (61,46 % ja 60,22 %; kuva 4b). Leu oli d-asemassa (kuva 4a) 45 prosentissa kaikista maapähkinän bZIP:istä, ja se on yksi dimeeriä eniten stabiloivista alifaattisista aminohapoista. E-asennossa 37 prosentissa kaikista maapähkinän bZIP:istä oli happamia aminohappoja D tai E, kun taas g-asennossa 44 prosentissa kaikista maapähkinän bZIP:istä oli emäksisiä aminohappoja R tai K (kuva 4a). Näiden varautuneiden aminohappojen ajatellaan muodostavan suolasiltoja kierteiden välille sähköstaattisissa vuorovaikutuksissa. Vetovoimaiset tai hylkivät g↔e′ -elektrostaattiset vuorovaikutukset voivat myös muodostaa kierteiden välisiä suolasiltoja, jotka vaikuttavat dimeroitumispesifisyyteen ja vakauteen . Tutkittaessa e- ja g-asemissa olevien varattujen jäännösten osuutta Arachis bZIP -proteiinien dimerisaatio-ominaisuuksien säätelyssä laskettiin vetovoimaisten ja hylkivien g↔e′ -parien frekvenssit kussakin heptadissa (kuva 4c). Kaikissa heptadeissa vetovoimaiset g↔e′-parit keskittyivät toiseen (15,6 %), viidenteen (35 %) ja kuudenteen (30 %) heptadiin, mikä osoittaa, että ne voivat muodostaa täydellisiä vetovoimaisia g↔e′-vuorovaikutussuhteita ja edistää stabiliteettia täydentymisen kautta heterodimeerissä. Kolme ryhmää, jotka käsittivät 28 alaperhettä (BZ1-BZ28), jaettiin edelleen homo- ja heterodimerisaatio-ominaisuuksien, erityisesti dimerisaatiospesifisyyden, perusteella (Lisätiedosto 7).

Koko genomin duplikaation ja tandemduplikaation vaikutus Arachiksen bZIP-geeniperheen laajenemiseen

Identifioimme koko genomin laajuiset kollineaariset duplikoituneet blokit A. duranensiksen (A. duranensis)- ja A. ipaensiksen (A. ipaensis)-genomeista sekä ortologiset kollineaariset blokit kahden genomin välillä. Laskettiin kollineaaristen lohkojen sisällä olevien paralogien ja ortologien pareittaiset synonyymiset etäisyydet (Ks-arvot) ja piirrettiin niiden frekvenssijakaumat (kuva 5a; Ks bin = 0,05). Huippu Ks-taajuus A. duranensiksen ja A. ipaensiksen välillä, joka edustaa keskimääräistä sekvenssivaihtelua, oli 0,035. Tämä edusti näiden kahden läheisesti sukua olevan Arachis-lajin välistä sekvenssieroa, jonka arvioitiin eronneen noin 2,16 miljoonaa vuotta sitten. Lisäksi A. duranensis- ja A. ipaensis -paralogien Ks-huiput olivat vastaavasti 0,90 ja 0,95, mikä vastaa ~ 58 miljoonaa vuotta sitten tapahtuneen varhaisen papilionoidien koko genomin duplikaatiotapahtuman (WGD) aiheuttamaa sekvenssieroa.

Kuvio. 5
kuvio5

Koko genomin duplikaatiosta (WGD) peräisin olevat Arachisin bZIP-geenit. a WGD:stä peräisin olevien paralogien Ks-jakauma WGD:stä peräisin olevista duplikoituneista genomilohkoista A. duranensis ja A. ipaensis. b WGD:stä johdetut duplikoidut bZIP-paralogit on yhdistetty sinisillä (A. duranensis) ja vihreillä (A. ipaensis) viivoilla. A. duranensiksen ja A. ipaensiksen väliset bZIP-ortologit linkitettiin lukuviivoilla

Havaitsimme 35 AdbZIP- ja 32 AibZIP-geeniä, jotka ovat osallisina duplikoituneissa genomilohkoissa, mikä vastaa noin 70 %:a (35/50) ja 71:tä %:a (32/45) bZIP-geeneistä kummassakin lajissa (Kuva 5b ja Additional file 8). Lisäksi monistuneet bZIP-geeniparit esiintyivät joko kromosomin sisällä tai kromosomien välillä, ja jotkin näistä pareista olivat segmentaalisesti monistuneet kerran, kahdesti tai kolmesti. Tämä tulos osoitti, että geenejä säilytetään ensisijaisesti ja että kromosomijärjestelyt ovat yleisiä WGD:n jälkeen. Tandemduplikaatioita havaittiin vain kahdella geeniparilla (AdbZIP33/AdbZIP34 ja AdbZIP41/AdbZIP42) A. duranensiksessa ja vain yhdellä geeniparilla (AibZIP28/AibZIP29) A. ipaensiksessa. Tämä viittaa siihen, että tandemduplikaatioita esiintyi harvoin eikä se ollut segmentaalista duplikaatiota tärkeämpi bZIP-geeniperheen laajentumisessa. Käytimme myös fylogeneettisiä ja synteettisiä analyysejä 35 ortologisen bZIP-geeniparin tunnistamiseksi A. duranensiksen ja A. ipaensiksen välillä. Nämä geenit olivat homeologeja myös tetraploidisen maapähkinän kahden alatyypin välillä.

Ymmärtääksemme Arachisin bZIP-geeneihin vaikuttavia evolutiivisia rajoitteita laskimme Ka/Ks-arvot kullekin duplikoituneelle bZIP-geeniparille kahdessa Arachis-lajeissa (Additional file 9). Useimmissa näistä parittaisista vertailuista Ka/Ks-arvot olivat alle 0,5 (vain yksi parittainen vertailu duplikoituneiden AdbZIP-geenien välillä ja vain kaksi duplikoituneiden AibZIP-geenien välillä oli suurempi kuin 0,5). Tämä viittasi siihen, että voimakas puhdistava valinta vaikutti Arachiksen duplikoituneisiin bZIP-geeneihin poistaakseen haitalliset mutaatiot proteiinitasolla.

Arachiksen bZIP-geenien ekspressioanalyysi maapähkinän siementen kehittymisen aikana

Profiilisoidaksemme bZIP-geenien ilmentymistä käytimme aiemmin julkaisemaamme RNA-seq-dataa , joka dokumentoi geenien ilmentymisen maapähkinän siemenissä eri kehitysvaiheissa: 20, 40 ja 60 päivää kukinnan jälkeen (DAF). Tämän datan avulla tunnistimme kaikkien Arachiksen bZIP-geenien FPKM-arvot ja kaikki eri tavoin ilmentyvät bZIP-geenit kolmessa kehitysvaiheessa. Lukuun ottamatta 24 bZIP-geeniä, jotka eivät ilmentyneet missään kehitysvaiheessa, tunnistettiin neljä ryhmää, jotka sisälsivät vastaavat bZIP-geenit, joilla oli spesifinen ilmentymisprofiili (kuva 6a ja lisätiedosto 10). Ensimmäiseen ryhmään kuului 37 bZIP-geeniä, jotka nousivat ylöspäin kehityksen alkuvaiheessa (20 DAF), mutta laskivat sen jälkeen (40 ja 60 DAF:ssa). Toiseen ryhmään kuului 15 bZIP:tä, jotka olivat ylössäätyneet 40 DAF:n kohdalla, kun taas kolmanteen ryhmään kuului 17 bZIP:tä, jotka olivat alasääntyneet 40 DAF:n kohdalla. Neljänteen ryhmään kuului 22 bZIP:tä, jotka ilmentyivät voimakkaasti kaikissa kolmessa kehitysvaiheessa. Ryhmän neljä voimakkaasti ekspressoituneet bZIP:t jakautuivat pääasiassa kladeihin A, C ja S. Useat näistä bZIP:istä olivat homologisia sellaisten geenien kanssa, joiden on todettu osallistuvan siementen kehitykseen muissa kasveissa, kuten Arabidopsiksessa, riisissä ja maississa. Tässä tutkimuksessa valittiin 12 bZIP-geeniä, jotka ilmentyivät voimakkaasti ja olivat homologisia aiemmin hyvin tutkittujen siementen kehitykseen vaikuttavien geenien kanssa, qRT-PCR-vahvistusta varten, ja havaittiin, että RNA-seq-menetelmällä määritetyt ilmentymismallit olivat yhdenmukaisia qRT-PCR-menetelmällä havaittujen mallien kanssa (Kuva. 6b.

Kuva 6
kuvio6

Arachiksen bZIP-geenien ilmentyminen maapähkinän siementen kehittymisen aikana. a Tunnistettiin neljä ryhmää (ryhmät I-IV), joihin kuuluivat vastaavat bZIP-geenit, joilla oli spesifinen ilmentymiskuva. Kussakin alaryhmässä harmaat viivat osoittivat bZIP-geenien ilmentymisarvot DAF20-, DAF40- ja DAF60-hetkellä. Punainen viiva osoittaa kaikkien bZIP-geenien keskimääräisen FPKM:n. b Siementen kehityksen aikana ilmentyvien 12 bZIP-geenin qRT-PCR-tarkistus. Näytetään qRT-PCR:llä (oranssit histogrammit) ja RNA-seq:llä (siniset viivat) mitatut geenien suhteelliset ilmentymistasot. Tulokset perustuvat kolmeen biologiseen toistoon; virhepalkit edustavat SE:tä. c A. duranensiksen ja A. ipaensiksen bZIP-ortologien ilmentymismalli siementen kehityksen aikana. Ortologien väliset samankaltaiset (merkitty Y:llä) tai poikkeavat (merkitty N:llä) ilmentymismallit on merkitty

Ryhmässä A AdbZIP33 ja AibZIP28 olivat ortologisia Arabidopsis ABA insensitive 5:n (ABI5) kanssa, joka liittyy ABA-signaalin välitykseen sekä siementen kehittymisen ja pitkäikäisyyden säätelyyn Arabidopsiksessa ja palkokasveissa . RNA-seq- ja qRT-PCR-tuloksemme osoittivat, että molemmat ortologiset ABI5-kopiot tetraploidin maapähkinän kahdesta alatyypistä ilmentyivät voimakkaasti kehityksen aikana, mikä viittaa siihen, että näiden geenien toiminta voi olla samanlaista maapähkinässä ja Arabidopsiksessa. qRT-PCR-tuloksemme osoittivat myös, että ryhmän A geenit AdbZIP42, AdbZIP48 ja AibZIP31 ilmentyivät vakaasti kehityksen aikana (kuva 6b ja lisätiedosto 11). Nämä geenit ovat homologisia ABF- ja AREB-geeneille, jotka osallistuvat ABA-välitteiseen siementen kehitykseen, itämiseen ja alkion kypsymiseen . Kolme ryhmän C geeniä (AdbZIP23, AdbZIP37 ja AibZIP30) ekspressoitui myös voimakkaasti, ja ne ovat homologisia maissin bZIP-tekijän Opaque2 kanssa. Opaque2 säätelee proteiinien kertymistä sekä aminohappo- ja sokeriaineenvaihduntaa maissin siemenissä . Lisäksi ryhmän S geenit AibZIP10, AdbZIP12, AdbZIP24, AdbZIP26 ja AdbZIP36 ilmentyivät erittäin voimakkaasti maapähkinän siemenissä (kuva 6b ja lisätiedosto 11). Mielenkiintoista oli, että ryhmän S geeneillä AdbZIP24 ja AdbZIP36 oli samanlainen ilmentymismalli kuin ryhmän C geeneillä AdbZIP37 ja AibZIP30: ilmentymistaso laski siementen kehityksen edetessä.

Tämän jälkeen selvitimme tarkemmin geenien ilmentymisen eroavaisuuksia tetraploidisen maapähkinän AA- ja BB-genomien homeologisten geenien välillä. Heatmap-analyysi osoitti, että yleiset ekspressiomallit siementen kehityksen aikana olivat samankaltaisia 31 AA- ja BB-genomin homeologisten/ortologisten geeniparien osalta. Käytimme differentiaalisen ilmentymisen analyysimenetelmää yhdessä tilastollisten menetelmien kanssa laskeaksemme erot geenien ilmentymisessä näiden geeniparien välillä kunkin näytteen osalta. Havaitsimme, että 3 geeniparia (AdbZIP5 ja AibZIP5, AdbZIP17 ja AibZIP15, AdbZIP46 ja AibZIP41) ilmentyi eri tavoin 20 DAF:ssa, 3 geeniparia (AdbZIP3 ja AibZIP1, AdbZIP4 ja AibZIP4, AdbZIP49 ja AibZIP45) 40 DAF:ssa ja 5 paria (AdbZIP3 ja AibZIP1, AdbZIP33 ja AibZIP28, AdbZIP37 ja AibZIP30, AdbZIP10 ja AibZIP10, AdbZIP1 ja AibZIP3) 60 DAF:ssa. Nämä tulokset osoittivat, että yleinen ilmentyminen säilyi kahden genomin välillä, mutta viittasivat siihen, että 20 %:n geenien ilmentyminen oli erilaistunut kahden genomin rinnakkaisen evoluution ja polyploidisaation aikana (kuva 6c).

Arachiksen bZIP-geenien qRT-PCR-ilmentymisprofiilit suolastressin alaisena

Käyttämällämme qRT-PCR:llä selvitimme bZIP-geenien ilmentymisessä tapahtuvia muutoksia vasteena suolakäsittelyn vaikutuksesta (kuva 7 ja Additional file 12). Emme pystyneet selvästi monistamaan 4 bZIP-geeniä PCR:llä. Kun maapähkinän juuria oli käsitelty suolalla 1 h, 20 geeniä ilmentyi merkitsevästi eri tavoin; 5 h jälkeen 27 geeniä ilmentyi merkitsevästi eri tavoin; ja 10 h jälkeen 41 geeniä ilmentyi merkitsevästi eri tavoin (Kuva 7j; Studentin t-testi: P < 0,05). Kussakin aikapisteessä paljon enemmän geenejä oli ylös- kuin alas-reguloitunut (14 vs. 6 1 h:ssa; 21 vs. 6 5 h:ssa; ja 34 vs. 7 10 h:ssa). Näistä suolakäsittelyn jälkeen eri tavoin ilmentyneistä bZIP-geeneistä monet jakautuivat ryhmiin A ja S (kuva 7k), mikä osoittaa, että näissä ryhmissä olevilla bZIP-geeneillä on tärkeä rooli sokerin signaloinnissa ja abioottisen stressin säätelyssä.

Kuva 7
kuvio7

Arachiksen bZIP-geenien ilmentymistasot maapähkinän juureksissa suolakäsittelyn jälkeen 0, 1, 5 ja 10 tunnin kuluttua. a-i bZIP-geenien ilmentymistasot eri ryhmissä. *: P < 0.05. j Merkitsevästi eri tavoin ilmentyneiden bZIP-geenien määrä kussakin ryhmässä. k Merkitsevästi eri tavoin ilmentyneiden bZIP-geenien määrä 1, 5 ja 10 h suolakäsittelyn jälkeen

Ryhmän A bZIP-geeneillä on stressi- ja/tai ABA-signaalin välittämiseen osallistuvia CKII- ja Ca2 + -riippuvaisten proteiinikinaasien fosforylaatiokohtien motiiveja, ja näillä motiiveilla on tärkeä merkitys kasvien sopeutumiselle erilaisiin abioottisiin ympäristöstressitekijöihin . Monet ryhmän A geenit liittyvätkin suolastressivasteeseen. Arabidopsiksessa ABI5 ja ABFs/AREB ovat keskeisiä ABA-riippuvaisia signaalinsiirtotekijöitä, jotka osallistuvat abioottisen stressin sietoon . GhABF2:n yli-ilmentäminen paransi merkittävästi suolastressin sietokykyä sekä Arabidopsiksessa että puuvillassa . Tomaatissa slAREB1:n ja slbZIP1:n knockout lisäsi suolastressin sietokykyä, kun taas slAREB1:n ja slbZIP1:n yli-ilmentäminen vähensi suolastressin sietokykyä . Tässä geenit AdbZIP42 ja AibZIP35 säätäytyivät merkittävästi suolastressin vaikutuksesta, ja nämä geenit ovat homologisia ABF:n, GhABF2:n, slAREB1:n ja slbZIP1:n kanssa. Lisäksi näiden geenien on raportoitu fosforyloituvan ABA:n aktivoimien SnRK2-proteiinikinaasien toimesta , mikä viittaa siihen, että ABA-vaste-elementtiä sitovien tekijöiden fosforylointi voi olla kriittistä ABA-välitteisen suolastressivasteen kannalta.

Ryhmän B geenit AdbZIP45 ja AibZIP40 nousivat ylöspäin 10 tuntia kestäneen suolastressin jälkeen, ja nämä geenit ovat homologisia AtbZIP17:n kanssa, mikä saattaisi tehostaa useiden suola-stressivasteeseen reagoivien geenien ilmentymää Arabidopsiksessa . Seitsemän ryhmän G bZIP-geeniä (AdbZIP7, AdbZIP15, AdbZIP19, AdbZIP50, AibZIP17, AibZIP21 ja AibZIP38) olivat homologisia Arabidopsiksen AtbZIP41:n ja tomaatin slbZIP38:n kanssa, ja molempien on osoitettu säätelevän suolastressiä negatiivisesti . Näistä seitsemästä geenistä AdbZIP15 säätyi merkittävästi alaspäin 1 tunnin ja 5 tunnin suolastressikäsittelyn jälkeen, kun taas AdbZIP19 ja AibZIP17 säätyivät merkittävästi ylöspäin 10 tunnin suolastressin jälkeen. AdbZIP15, AdbZIP19 ja AibZIP17 saattavat siis lisätä suolastressin vastustuskykyä. AdbZIP15 saattaa olla suolastressin negatiivinen säätelijä, sillä sen ilmentymismalli oli samanlainen kuin slbZIP38:n ilmentymismalli vasteena suolastressiin.

Ryhmän S geenit AdbZIP24 ja AdbZIP36 olivat homologisia AtbZIP1:n, AtbZIP53:n, MtbZIP2:n ja MtbZIP26:n kanssa, ja niiden ilmentymismallit vasteena suolastressiin olivat samankaltaiset (kuva 7). Erityisesti AdbZIP36:n säätely lisääntyi merkittävästi 10 tunnin suolastressin jälkeen. Arabidopsiksen kahden homologisen geenin, AtbZIP1:n ja AtbZIP53:n, osoitettiin ohjelmoivan uudelleen primaarisen hiilihydraatti- ja aminohappoaineenvaihdunnan, jotta juuret sopeutuisivat suolastressiin. Myös homologit MtbZIP2 ja MtbZIP26 indusoituvat transkriptiivisesti suolakäsittelyn vaikutuksesta ja parantavat kasvin sietokykyä suolastressiä vastaan. Erityisesti AdbZIP36:n ilmentymismalli oli samanlainen kuin AtbZIP1:n, MtbZIP2:n ja MtbZIP26:n ilmentymismalli Arabidopsiksessa ja M. truncatulassa , mikä viittaa siihen, että AdbZIP36 saattaa olla positiivinen säätelijä suolastressin sietokyvylle maapähkinässä. Yhteenvetona voidaan todeta, että ekspressioanalyysitutkimuksessamme on tunnistettu useita maapähkinän bZIP-ehdokkaita, jotka saattavat liittyä suolastressivasteeseen ja jotka ovat tulevan tutkimuksen kohteita.

Vastaa

Sähköpostiosoitettasi ei julkaista.