بررسی بیان چهار همسانه ژنی عضو خانواده MADS- box در گیاه دوپایه ترشک (Rumex acetosa L.)

بررسی بیان چهار همسانه ژنی عضو خانواده MADS- box در گیاه دوپایه ترشک (Rumex acetosa L.)

علی محمد شکیب¹*، چارلز اینزورث² 

¹ پژوهشکده بیوتکنولوژی کشاورزی ایران- کرج

² آزمایشگاه بیولوژی ملکولی گیاهی، ایمپریال کالج، انگلستان

چکیده

بررسی ژن­هایی که در تشکیل گل دخالت دارند می­توانند ما را در فهم بیشتر مکانیسم­های تعیین  جنسیت کمک نماید. در این تحقیق، الگوی بیان چهار ژن از اعضای خانواده ژنی MADS- box شامل RaA2 ، RaB17 ، RaG24 و RaK16 جداسازی شده از گیاه دو پایه ترشک ( Rumex acetosa L.) مورد مطالعه قرار گرفتند. برای استخراج RNA از بافت­های مختلف (گل آذین، ساقه، برگ و ریشه) گیاهان نر و ماده برای آزمون نورترن بلات استفاده گردید. ژن­های RaA2 و RaB17 در اعضای زایشی گیاه بیان می­شوند اما ژن­های RaG24 و RaK16 هم در اعضای رویشی و هم در اعضای زایشی بیان می­شوند. هم ردیفی توالی پروتئینی این همسانه­های ژنی با پروتئین­های MADS- box شناخته شده سایر گیاهان نشان داد که آنها متعلق به زیر خانواده­های مجزایی هستند. توالی پروتئینی ژن­های RaA2 و RaB17 بیشترین تشابه را به ترتیب با پروتئین­های MADS- box سیب (Malus domestica) و CMB1  گیاه Diantus caryophyllus نشان دادند. توالی پروتئینی ژن­های RaG24و RaK16 بیشترین تشابه را به ترتیب با پروتئین TDR3 گوجه فرنگی (Lycopersicon esculentum) و پروتئین MADS5  گیاه Betula pendula و پروتئین POTM1 سیب زمینی (Solanum tuberosum) نشان داد که هر دو در مراحل رویشی و زایشی بیان می­شوند. در آزمون سادرن بلات الگوی دورگه سازی متفاوتی بین نر و ماده وجود داشت. الگوی دورگه سازی سادرن بلات نشان می دهد که این ژن­ها به حالت تک کپی و به صورت خانواده­های ژنی کوچک در ژنوم قرار دارند.

کلمات کلیدی: گیاه دو پایه، بیان ژن،  MADS- box، Rumex acetosa L

مقدمه

ژن­های MADS- box خانواده ژنی بزرگی هستند که تقریبا در تمام موجودات پر سلولی یافت می­شوند (Krogan et al., 2000; Martinez et al., 2003; Singer et al., 2007; Svensson and Engstron, 2002; Tanabe et al., 2005; Tanabe et al., 2003). در بازدانگان، خیلی از ژن­های خانواده MADS- box در مراحل مختلف نمو گل، به ویژه در تعیین هویت مریستم گل و اعضاء گل نقش دارند (Davies et al., 1996; Parenicova et al., 2003; Reiechman and Meyerowitz, 1997 ). بیشتر اعضای این خانواده دارای بخش مشترک بنام MADS- box هستند (Schwarz-Sommer et al., 1990) که به DNA متصل می­شوند (Huang et al., 1995; Shiraishi et al., 1993; Shore and Sharrocks, 1995). جداسازی و تعیین نقش ژن­های این خانواده در موجودات مختلف در حال بررسی است (Henschel et al., 2002; Singer et al., 2007; Sevensson and Engestron, 2002). بخش مشترک بین ژن­های MADS- box به عنوان ابزاری برای جداسازی ژن­های این خانواده در گونه­های مختلف گیاهی مورد استفاده قرار می­گیرد (Kater et al., 1998; Puneli et al., 1991). این روش به ویژه در گیاهانی که جهش یافته­های فنوتیپی در آنها مشاهده نشده است، ارزشمند است. یکی از موضوعات مورد علاقه محققین، مطالعه مکانیزم­های تظاهر جنسیت و نحوه تکامل آن در گونه­های گیاهی است (Charlesworth, 2002; Singer et al., 2007; Tanurdzic and Banks, 2004). دو گیاه دو پایهSilene latifolia  و Rumex acetosaبه عنوان گیاهان مدل در مطالعات تظاهر جنسیت مورد استفاده قرار می­گیرند (Ainsworth et al., 2005). در گیاه دو پایه ترشک (Rumex acetosa L. )، که در آن گل­های نر و ماده به صورت جدا گانه روی گیاهان منفرد تشکیل می­شود، مطالعه ژن­های دخیل در تمایز اعضای گل می­تواند در فهم مکانیزم­های تعیین جنسیت گیاهان کمک نماید. قبلا با هدف جداسازی ژن­های MADS- box، یک خزانه cDNA ساخته شده از گل آذین­های نر و ماده ترشک با استفاده از دو ژن DEFICIENS و PLENA به عنوان کاوشگر غربال گردیدند (Ainsworth et al., 1995). هفت همسانه با طول کامل انتخاب شدند که سه مورد از این همسانه­ها به نام­های RAD1 (Rumex acetosa DEFICIENS-like)، RAD2 و RAP1 (Rumex acetosa PLENA-like) مورد بررسی تفصیلی قرار گرفته (Ainsworth et al., 1995) و نشان داده­اند که در پریموردیاهای اعضای جنسی گل در پایه­های نر و ماده ترشک بیان می­شوند. با اینکه تا به حال چندین ژن از ترشک جداسازی و مشخصات مولکولی آنها تعیین شده اما نقش آنها بدلیل نبود روش انتقال ژن به ترشک تعیین نشده است. تعدادی از ژن­های این خانواده از گونه­های گیاهی دیگر جداسازی شده و نقش آنها تعیین گردیده است. هم اکنون تعداد بسیار زیادی توالی ژنی در بانک اطلاعات DNA وجود دارد که می­توان در حد امکان، از آنها برای کسب اطلاعات بیشتر در باره تشابه توالی و پیش بینی کارکرد ژن­های جدیداً جداسازی شده یا ژن­هایی که کارکرد ناشناخته دارند، استفاده کرد. این اطلاعات می­تواند در بحث­های تکاملی و جنبه­های مهندسی ژنتیک کاربرد داشته باشند. در اینجا خصوصیات چهار همسانه ژنی جدا شده از Rumex acetosa شامل: (RaA2)، (RaB17)، (RaG24) و  (RaK16)از نظر تشابه توالی و بررسی بیان تشریح شده است. 

مواد و روش­ها

گیاهان نر و ماده ترشک Rumex acetosa L.)) به عنوان مواد گیاهی مورد استفاده قرار گرفتند. بررسی­ها در آزمایشگاه زیست شناسی مولکولی، امپریال کالج مستقر در وای، انگلستان صورت گرفت.   

استخراج DNA، RNA: DNA برگی گیاهان نر و ماده و RNA کل ریشه، ساقه، برگ و گل آذین گیاهان نر و ماده به طور جداگانه مطابق روش­های توصیه شده توسط اینزورث و همکاران (Ainsworth et al., 1995; Ainsworth, 1994 ) جداسازی گردیدند.

تهیه کاوشگر: توالی کاوشگرهای مورد استفاده از بخش´3 همسانه­ها، بدون ناحیه MADS، تهیه شدند. برای همسانه RaA2 یک قطعه 480 جفت بازی با استفاده از برش آنزیمی Sst1 بکار رفت. برای همسانه Ra17 یک قطعه 780 جفت بازی با برش آنزیم­های Kpn1/ EcoR1 و برای همسانه RaG24 و همسانه RaK16 به ترتیب یک قطعه 740 جفت بازی با برش آنزیم  HindIIIو یک قطعه 530 جفت بازی با برش آنزیم­های Kpn1/ Xba1 مورد استفاده قرار گرفت. قطعات DNA بالا سپس با ³²P نشاندار شدند.

دورگه سازی سادرن و نوردرن بلات: RNA کل تهیه شده از بافت­های مختلف روی ژل آگارز فرمالدهید 2/1 درصد جداسازی گردید.  سپس RNA به غشاء نایلونی منتقل و با کاوشگر نشاندار شده دورگه شد. غشاء­ها با بافر شستشو با شدت بالا (0.2 x SSC) و دمای  ^oC65 شستشو و مجاور فیلم X-ray  قرار داده شدند.

نتایج و بحث

توالی­های اسید آمینه­ای پیش­ بینی شده برای چهار همسانه MADS- box ترشک با استفاده از جستجویBLAST  برای تعیین اینکه آیا شباهتی با توالی­های موجود در بانک اطلاعاتی دارند، مورد مقایسه قرار گرفتند. پیدا کردن شبیه­ترین توالی­ها برای همردیفی توالی، با استفاده از برنامه MegAlign انجام و درخت فیلوژنتیک با استفاده از روش کلاستال (Clustal) تهیه گردید. بر اساس این بررسی، پروتئین­های MADS- box ترشک در گروه های اصلی SQUA/AP1، MdMADS1/CMB1،  TDR3/SaMADSA و POTM1/SLM5 قرار گرفتند (شکل1). همسانه RaA2 شامل یک cDNA به طول 1113 جفت باز است که یک پروتئین با 243 اسید امینه به وزن 2/28 کیلو دالتون رمز می­کند. مقایسه انجام گرفته با استفاده از بانک اطلاعات پروتئین نشان داد که ژن رمز کننده این پروتئین، یعنی RaA2 به گروه SQUA/AP1 نزدیک بوده (شکل1) و تشابه توالی 9/61 درصدی با پروتئین MADS- box از سیب (Malus domestica) (شماره بازیابی (AJ000759 و تشابه توالی 3/61 درصدی با پروتئین (SQUA) SQUAMOSA از گل میمون (Antirrhinum majus) (Huijser et al., 1992) و پروتئین (AP1) APETALA1 ازArabidopsis thaliana  (Martinez et al., 2003) دارد. همچنین تشابه توالی 1/60 درصدی با SaMADS C ، هومولوگ AP1، از Sinapis alba (Menzel et al., 1995) و تشابه توالی 3/59 درصدی با پروتئینBoiAP1 ازBrassica oleracea  (شماره بازیابی U67451) (Car and Irish, 1997) نشان داد. در این گروه، بالاترین تشابه (3/95 درصد و 5/94 درصد) بین پروتئینی ازB. oleracea  با شماره بازیابی U67451 و پروتئین­های AP1 از  S. albaو A. thaliana  مشاهده گردید. SQUA و AP1با همدیگر 5/64 درصد تشابه نشان می­دهند. اعضای این گروه عموما در تعیین هویت مریستم گل دخالت دارند. شکل­های 2 و 3 به ترتیب همردیفی و درصد تشابه توالی پروتئین RaA2 با دیگر پروتئین­های این گروه را نشان می­دهد. همسانهRaB17  دارای 1140 جفت باز است که پیش بینی می­شود یک پروتئین با 251 اسید آمینه به وزن 8/28 کیلو دالتون با کدون آغاز ترجمه (ATG) در موقعیت 159، از آن رمز شود. مقایسه پروتئین RaB17 با دیگر پروتئین­های MADS نشان داد کهRaB17  می­تواند در گروه MdMADS1/CMB1  قرار گیرد (شکل1). در این گروه، RaB17 تشابه توالی 9/63 درصدی با پروتئین CMB1 از Dianthus caryophyllus  با شماره بازیابی Q39685 و تشابه توالی 7/54 درصدی با PrMADS1 از  Pinus radiata با شماره بازیابی U42399، تشابه توالی 2/52 درصدی با پروتئین MADS- box گل میمون با شماره بازیابی X95467 و تشابه توالی 6/51 درصدی با پروتئین MdMADS1 سیب (Malus domestica) (Sung and An, 1997) دارد. همردیفی توالی پروتئین پیش بینی شده RaB17 با سایر پروتئین­های MADS در شکل4 نشان داده شده است. مانند دیگر پروتئین­های  MADSگیاهی تشابه بالایی در بخش MADS بین پروتئین  B17با سایر اعضای خانواده MADS و تشابه کمی در توالی­های پائین دست بخش MADS وجود دارد. در این گروه دو پروتئین از سیب (MdMADS8 و MdMADS1) بالاترین تشابه (6/99 درصد) را دارند. پروتئین MADS از گل میمون (با شماره بازیابی X9546) به ترتیب دارای تشابه 7/66 درصدی و 3/66 درصدی با پروتئین­های سیب هستد. پروتئین CMB1 تشابه 9/63 درصدی با PrMADS1 دارد. ژن MdMADS1 در اعضای گل و میوه های جوان بیان می­شود (Sung and An, 1997). با اینحال، اطلاعات زیادی در باره الگوی بیان یا کارکرد ژن­های این گروه وجود ندارد.

 همسانهRaG24  دارای یک cDNA به طول 1122 جفت باز است که پیش بینی می­شود یک پروتئین 214 اسید آمینه­ای با وزن 4/32 کیلو دالتون رمز کند که از نوکلئوتید 5 آغاز و در نوکلئوتید 646 خاتمه یابد. مقایسه توالی پروتئینRaG24  با دیگر پروتئین­های MADS نشان داد که این پروتئین متعلق به گروه TDR3/SaMADS A می­باشد (شکل1). اعضای این گروه عموما در هر دو اعضای رویشی و زایشی بیان می­شوند و کارکرد­های مختلفی مانند دخالت در تشکیل مریستم گل آذین و گل دارند. پروتئین RaG24 دارای تشابه 2/64 درصدی با TDR3 از گوجه فرنگی (Lycopersicun esculentum) (Pnueli et al., 1991)، تشابه 6/62 درصدی با  AGL20از Arabidopsis thaliana (با شماره بازیابی AC003680)، تشابه 9/61 درصدی با TobMADS1 از توتون (Mandel et al., 1994) و تشابه 5/61 درصدی با SaMADS A از Sinapis alba (Menzel et al., 1996) می­باشد. در این گروه، AGL20 تشابه 4/94 درصدی با  SaMADS Aو TRD3 تشابه 5/95 درصدی با  TobMADS1دارند. پائین­ترین درصد تشابه (3/43 درصد) بین پروتئین PrMADS9 از Pinus radiate (با شماره بازیابی U90344) و پروتئین AGL14 از Arabidopsis thaliana (با شماره بازیابی U20184) می­باشد. بررسی همردیفی پروتئین RaG24 با سایر پروتئین­های MADS گیاهی در شکل5 نشان داده شده است.

شکل 1- درخت فیلوژنی بر اساس همردیفی پروتئین­های MADS ترشک (به صورت برجسته) با شبیه­ترین پروتئین­های MADS با استفاده از روش کلاستال (Clustal). BoiAP1 از Brassica oleracea (شماره بازیابی U67451) (Car and Irish, 1997)، SaMADS C از Sinapis alba  (شماره بازیابی X81480) (Menzel et al., 1995)،AP1  ازArabidopsis thaliana  (شماره بازیابی S27109) (Mendel et al., 1992)، SQUA از Antirrhinum majus  (شماره بازیابی S20886) (Huijser et al., 1992)، پروتئینMADS  از Malus domestica (شماره بازیابی AJ000759)، RaA2; ازRumex acetosa ، همسانAGL8  از Solanum commersonii (شماره بازیابی AF002666)،POTM1  ازSolanum tuberosum  (شماره بازیابی Q42429) (Kang and Hannapel, 1995)، NAP1 از Nicotiana tabacum  (شماره بازیابی AF009126)،MADS5  از Betula pendula  (شماره بازیابی X99655)،NAP1  از Nicotiana tabacum (شماره بازیابی AF009126)، SLM5  از Silene latifolia  (شماره بازیابی X80492)، (Hardenack et al., 1994)،RaK16  از Rumex acetosa ،MdMADS8  ازMalus domestica  (شماره بازیابی AJ001681)،MdMADS1  از Malus domestica   (شماره بازیابی U78947) (Sung and An, 1997)، پروتئینMADS  از Antirrhinum majus  (شماره بازیابی X95467) ،MdMADS1  از Malus domestica (شماره بازیابی U78947)،PrMADS1  از Pinus radiata  (شماره بازیابی U42399)، CMB1از Dianthus caryophyllus (شماره بازیابی Q39685)، RaB17; از Rumex acetosa ، TDR3 از Lycopersicon esculentum  (شماره بازیابی X60756) )(Pnueli et al., 1991)، TobMADS1 از Nicotiana tabacum  (شماره بازیابی X76188) )(Mandel et al., 1994)، پروتئین  MADS ازPimpinella brachycarpa  (شماره بازیابی AF082531)،AGL20  از Arabidopsis thaliana  (شماره بازیابی AC003680)، SaMADS A  از Sinapis alba  (شماره بازیابی U25696) (Menzel et al., 1996)، RaG24  از Rumex acetosa ،AGL14  از Arabidopsis thaliana  (شماره بازیابی U20184) (Rounsley et al., 1995)، PrMADS9 از  Pinus radiata (شماره بازیابی U90344)، RAP1 از Rumex acetosa  (شماره بازیابی X89107)( Ainsworth et al., 1995)، RAD1  از Rumex acetosa  (شماره بازیابی X89113) (Ainsworth et al., 1995).                       

Figure 1- Phylogenetic tree produced from an alignment of sorrel MADS-box proteins (in bold) with the most similar MADS-box proteins using Clustal method.

شکل 2 – همردیفی توالی­های اسید امینه همسانه ژنی RaA2 ترشک با سایر پروتئن­های MADS. SQUA  از  Antirrhinum majus(شماره بازیابی S20886) (Huijser et al., 1992)، AP1 از Arabidopsis thaliana (شماره بازیابی S27109)،SaMADS  از Sinapis albas (شماره بازیابی X81480) (Menzel et al., 1995)، پروتئین MADS از Malus domestica (شماره بازیابی AJ000759)، BoiAP1  از Brassica oleracea (شماره بازیابی U67451) (Car and Irish, 1997).

Figure 2- Alignment of amino acid sequence of the sorrel RaA2 clone and other MADS-box proteins.

شکل 3- درصد تشابه پروتئین همسانه ژنی RaA2 ترشک با پروتئن­هایSQUA ، AP1،   SaMADS ، MADSC، پروتئینMADS  از M. domestica و BoiAP1 .

Figure 3- Sequence similarity percentage of the sorrel RaA2 clone with other MADS proteins.

شکل 4 – همردیفی توالی­های اسید امینه همسانه ژنی RaB17 ترشک با سایر پروتئن­های MADS.  CMB1  از Dianthus caryophyllus (شماره بازیابی Q39685)،PrMADS1  از Pinus radiata  (شماره بازیابی U42399)، MdMADS8 از Malus domestica (شماره بازیابی AJ001681) MdMADS1 از Malus domestica (شماره بازیابی U78947) (Sung and An, 1997)، پروتئینMADS  از Antirrhinum majus (شماره بازیابی X95467).  

Figure 4- Alignment of amino acid sequence of the sorrel RaB17 clone and other MADS-box proteins.

شکل 5 – همردیفی توالی­های اسید امینه همسانه ژنی RaG24 ترشک با سایر پروتئن­های MADS.  TDR3 ازLycopersicon esculentum  (شماره بازیابی X60756) (Penueli et al., 1991)، پروتئینAGL20  از Arabidopsis thaliana (شماره بازیابی AC003680)،TobMADS1  از Nicotiana tabacum  (شماره بازیابی  X76188) (Mandel et al., 1994)، SaMADS A از Sinapis alba  (شماره بازیابی U25696) (Menzel et al., 1996)، پروتئینMADS  از Pimpinella brachycarpa   (شماره بازیابی AF082531)، PrMADS9 ازPinus radiata  (شماره بازیابی U90344)، AGL14 از Arabidopsis thaliana (شماره بازیابی U20184) (Rounsley et al., 1995).

Figure 5- Alignment of amino acid sequence of the sorrel RaG24 clone and other MADS-box proteins.

همسانه RaK16 ترشک دارای 969 جفت باز است که می­تواند پروتئینی با 225 اسید آمینه با وزن 8/25 کیلودالتون را رمز نماید. مقایسه پروتئین  RaK16 با سایر پروتئین­های MADS  گیاهی نشان داد که این پروتئین به پروتئین­های گروهPOTM1/SLM5   خیلی نزدیک است (شکل1). RaK16 تشابه توالی1/63 درصدی با MADS5 از Betula pendula (با شماره بازیابی X99655) و POTM1 از Solanum tuberosum (Kang and Hannapel, 1995) دارد. پروتئین RaK16 تشابه 2/62 درصدی با NAP1 از توتون Nicotiana tabacum (با شماره بازیابی AF009126) و هومولوگ AGL8 از Solanum commersonii (با شماره بازیابی AF002666) نیز نشان داد. این پروتئین همچنین تشابه 2/60 درصدی با  SLM5 از Silene latifolia (Hardenack et al., 1994) دارد. در بین اعضای این گروه، POTM1 بالاترین تشابه (96 درصد) را با پروتئینی از Solanum commersonii و نیز تشابه 3/76  درصدی با NAP1 از N. tabacum دارد. MADS5 تشابه 3/69 درصدی با POTM1 و نیز تشابه 9/66 درصدی با SLM5 دارد.

برای بررسی الگوی بیان همسانه­های ژنی ترشک، RNA از بافت­های مختلف شامل ریشه، ساقه، برگ و گل آذین گیاهان نر و ماده استخراج و پس از جداسازی روی ژل به غشاء نایلونی منتقل و غشاء­ها با کاوشگر مربوطه دورگه سازی شدند. همسانه­های RaA2 و RaB17 در نمونه­های RNA­ی بافت گل آذین هر دو گیاه نر و ماده بیان بالایی نشان دادند (شکل­های A6 و B6) ولی بیان آنها در سایر بافت­ها دیده نشد. با اینحال، در همسانه RaB17 علایم بیان خیلی ضعیفی در بافت­های رویشی دیده شد، که می­تواند ناشی از تشابه توالی با دیگر mRNAها باشد. همسانه­هایRaG24  وRaK16  در همه بافت­ها بیان شدند (شکل­های C6 و D6). اگرچه بسته به نوع بافت، سطح بیان همسانه­ها تفاوت­هایی نشان داد. بیان RaG24 در بافت گل آذین ماده نسبت به گیاهان نر شدیدتر بود. همه بافت­های گیاهان نر در مقایسه با بافت­های مشابه در گیاهان ماده سطح بیان پائین نشان دادند. سطح بیانRaG24  در ساقه بالاتر از میزان آن در برگ­ها بود و میزان بیان آن در ریشه­ها کمترین مقدار را داشت. بیان  RaK16در بافت ساقه و گل آذین گیاهان ماده بالاتر از سایر بافت­ها بود.

به منظور تعیین تعداد کپی این ژن­ها، آنالیز سادرن بلات پس از هضم با آنزیم HindIII یا EcoRI انجام گرفت (شکل7). الگوهای دورگه سازی بین گیاهان نر و ماده تفاوت دارند. این تفاوت­ها ممکن است مربوط به چند شکلی­هایی باشد که معمولا در یک گونه دگرگشن انتظار می­رود. الگوی دورگه سازی نشان از وجود یک کپی منفرد و خانواده ژنی کوچک از هر یک از ژن­ها در ژنوم می­دهد. 

شکل 6 – آنالیز نوردرن بلات. مقدار 20 میکروگرم RNA از بافت­های گل آذین، برگ، ساقه و ریشه گیاهان نر (M) و ماده (F) الکتروفورز شده وسپس به غشاء نایلونی منتقل شدند و با کاوشگر­های RaA2 (A) ، RaB17 (B) ، RaG24 (C) و RaK16 (D) دورگه سازی گردید.

Figure 6- Northern blot analysis. 20 µg RNA from inflorescence, leaf, stem and root tissues of male (M) and female (F) plants after electrophoresis were transferred to membrane. Hybridizations were carried out using RaA2 (A), RaB17 (B), RaG24 (C) and RaK16 (D) as probes.

شکل 7 – آزمون سادرن بلات. مقدار 10 میکروگرم DNA از گیاهان نر (M) و ماده (F) پس از هضم با HindIII (A، B و C) یا EcoRI (D) برای دورگه سازی مورد استفاده قرار گرفت. کاوشگر­های مورد استفاده شامل RaA2 (A)، RaB17 (B)، RaG24 (C) و RaK16 (D) بودند.

Figure 7- Southern blot analysis. 10 µg DNA from male (M) and female (F) plants after digestion with HindIII (A, B and C) or EcoRI (D) were used for hybridization. RaA2 (A), RaB17 (B), RaG24 (C) and RaK16 (D) were used as probes.

بحث

الگوی بیان ژن­های ترشک و مقایسه توالی پروتئینی آنها با ژن­های شناخته شده از سایر گیاهان بررسی گردید. بر اساس مقایسه توالی پروتئین pRaA2، ترشک تشابه توالی بالایی با پروتئین  MADSسیب (Malus domestica)، پروتئین SQUA گل میمون (Antirrhinum majus) و پروتئین AP1 علف تال دارد. در گل میمون، پروتئین SQUA در مریستم­های گل به میزان بالایی بیان می شود و بیان آن تا مراحل بعدی ریخت زایی اعضای گل بجز در تمایز پرچم تداوم می­یابد (Huijser et al., 1992). جهش در ژن SQUA باعث ایجاد گل­هایی با تعداد زیاد براکت و گل­های تغییر شکل یافته و ناقص می­شود (Huijser et al., 1992). در علف تال، ژنAP1  در پریموردیاهای گل و اعضای گل مانند کاسبرگ­ها و گلبرگ­ها بیان می­شود. بیان این ژن برای تبدیل مریستم گل آذین به مریستم گل و نیز تعیین هویت کاسبرگ و گلبرگ لازم است (Mandel et al., 1992). محصول ژنRaA2  نیز با پروتئین هومیوتیکBoiAP1  کلم Brassica oleracea   (Shiraishi et al., 1993) و پروتئین SaMADS C  گیاهSinapis alba  (Menzel et al., 1995) مشابهت نشان می­دهد. در کلم، ژنBoiAP1  در کاسبرگ­ها، گلبرگ­ها و در پریموردیاهای پرچم بیان می­شود. ژن SaMADS C نیز همانند ژنAP1  گیاه  S. albaابتدا در پریموردیاهای گل بیان می­شود (Menzel et al., 1995). آزمون­های نوردرن بلات نشان دادند که ژنA2  ترشک به میزان بالایی در گل آذین­های نر و ماده بیان می­شود چرا که mRNA مربوط به ژن A2 در ریشه­ها، ساقه­ها و برگ­ها قابل ردیابی نبود. بر اساس الگوی بیان و تشابه توالی، به نظر می­رسد که ژن RaA2 در مرحله زایشی بیان و احتمالاً کارکردی همانند ژن SQUA داشته باشد. پروتئین RaB17 تشابه توالی بالایی با پروتئین­هایCMB1  از Dianthus caryophyllus، PrMADS1 از Pinus radiate، یک پروتئین MADS گل میمون و پروتئین­های MdMADS1 وMdMADS8  از سیب (M. domestica) نشان داده است. ژنMdMADS1  در همه اعضای گل و میوه­های جوان بجز در برگ­ها بیان می­شود (Sung and An, 1997). بیان بالاتر این ژن در مراحل اولیه نمو گل و میوه دلالت بر آن دارد که این ژن در شروع نمو اعضای زایشی نقش عمده­ای ایفا می­کند. آنالیز نوردرن بلات نشان داد که ژنRaB17  در زمان نمو گل به شدت بیان می­شود و بیان آن در هر دو گیاه نر و ماده بیانگر آن است که این ژن در نمو گل گیاهان نر و ماده نقش دارد. پروتئینRaB17  با بعضی از ژن­های MADS  که کارکرد آنها روشن نشده است، تشابه توالی دارد. پروتئین RaG24 با پروتئین­هایTDR3  گوجه فرنگی (Lycopersicon esculentum ) (Pnueli et al., 1991)،AGL20  علف تال،TobMADS1  توتون (Mandel et al., 1994) و  SaMADS A گیاهSinapis alba  (Menzel et al., 1996) بیشترین تشابه را دارد. ژنTDR3  در برگ­ها و مریستم­های گل اما نه در اعضای گل گوجه فرنگی بیان می­شود (Pnueli et al., 1991). ژنTobMADS1  در هردو اعضای رویشی و زایشی توتون بیان می­شود ولی بیان آن در برگ­های بالغ و گلبرگ­ها نسبتا پائین است (Mandel et al., 1994). بیان ژن  SaMADS Aپس از القا گلدهی در مریستم­های انتهایی و در مریستم­های گل پس از نمو پریموردیاهای کاسبرگ، مشاهده شده است (Menzel et al., 1996). ژن RaG24 با ژنTobMADS1  فوق الذکر نیز تشابه بالایی دارد. پروتئین این ژن همچنین مشابه پروتئین  AGLعلف تال می­باشد که تشابه بالایی با پروتئین AGL14 داشته و به طور اختصاصی در ریشه بیان می­شود (Rounsley et al., 1995). آزمون نوردرن بلات نشان داد که ژن RaG24 در بافت­های مختلف اعضای رویشی و زایشی بیان می­شود. معمولا، بیان ژن  RaG24در بافت­های گیاه ماده شدیدتر از بافت­های گیاه نر بوده و بالاترین میزان نسخه­برداری را در ساقه و گل آذین نشان می­دهد. در گیاه سیلن (White Campion) بیان هر دو ژن  SLM4و SLM5 در گل­های گیاه ماده نسبت به گیاه نر بیشتر بود اما میزان بیان ژن­هایSLM2  وSLM3  در گل­های نر نسبت به ماده بالاتر بوده است (Hardenack et al., 1994). الگوی بیانی ژنRaG24  در اعضای رویشی و زایشی و تشابه توالی با بعضی از ژن­های خانواده MADS بیانگر آن است که این ژن در هر دو مرحله رویشی و مرحله نمو گل نقش دارد. پروتئین RaK16، همانند پروتئین AGL8، بالاترین تشابه را با پروتئین POTM1 از سیب زمینی (Solanum tuberosum) دارد که بیان آن در اعضای رویشی بالا است (Kang and Hannapel, 1995). پروتئین RaK16 نیز همانندSQUA/AP1 با پروتئین MADS5 گیاه Betula pendula، پروتئین  NAP1توتون و SLM5 سیلن (Silene latifolia) تشابه نشان داد. ژنSLM5  مشابه ژن­های  SQUAو AP1 بوده و در تمام مریستم گل، و در ادامه در کاسبرگ­ها و گلبرگ­ها بیان می­شود (Kempin et al., 1993). بیان  RaK16در اعضای گل دیده نشد و در بافت­های گل آذین، برگ، ساقه و ریشه مشاهده گردید. این نتایج نشان می­دهد که این ژن احتمالا کارکرد­های مختلفی طی نمو رویشی و زایشی دارد. آنالیز الگو­های بیان و تشابه توالی 4 ژن عضو خانواده MADS- box جدا شده از ترشک نشان داد که ژن­های RaA2 و RaB17 در نمو گل نقش داشته و ژن­های RaG24 و RaK16 در هر دو مرحله رویشی و زایشی کارکردهایی دارند، با اینحال، نقش واقعی آنها را فقط می­توان با آنالیز این ژن­ها در گیاهان جهش یافته یا گیاهان تراریخته تعیین نمود.