شناسایی مسیرها و ژنهای دارای بیان افتراقی در تنش سرما در نیشکر، با استفاده از شاخصهای مورفولوژیکی، بیوشیمیایی و آنالیز rna-seq

با توجه به حساسیت نیشکر به سرما، به‌منظور شناسایی مسیرها و ژن‌های دارای بیان افتراقی در ارقام نیشکر در هنگام تنش سرما و استفاده از آن‌ها در برنامه اصلاح نبات نیشکر، پس از وقوع سرمای 1.2 درجه سانتی گراد در دی‌ماه 1394، تعداد 454 رقم نیشکر با استفاده از شاخص‌های مورفولوژیکی، بیوشیمیایی و آنالیز rnaseq، در مزرعه تحقیقاتی موسسه تحقیقات نیشکر خوزستان موردبررسی قرار گرفت. در مرحله اول با استفاده از شاخص‌های مورفولوژیکی، ارقام متحمل و حساس در برابر سرما، انتخاب شدند. در مرحله دوم در سال 1395، پس از تنش سرما، شاخص‌های بیوشیمیایی از قبیل پرولین و مالون دی آلدئید، در ارقام متحمل افزایش بیشتری پیدا کرد و بر این اساس دو رقم br0001 و tuc66107 به ترتیب به‌عنوان متحمل‌ترین و حساس‌ترین رقم نسبت به سرما شناسایی شدند. سپس این ارقام در سردخانه تحت تنش دمای 2 درجه سانتی گراد به مدت 12 ساعت قرار گرفتند. در مرحله سوم، استخراج rna، توالی یابی به روش rnaseq و ساخت کتابخانه cdna انجام گردید. سپس آنالیز بیان و تفسیر ژن‌های دارای بیان افتراقی انجام گردید. نتایج بررسی واکنش مورفولوژیکی و بیوشیمیایی ارقام نیشکر به تنش سرما نشان داد ارقامی که در مرحله مورفولوژیکی میزان تحمل بیشتری به تنش سرما دارند، همین ارقام در مرحله بررسی بیوشیمیایی نیز ازنظر میزان پرولین و mda در سطح بالاتری نسبت به ارقام حساس قرار می‌گیرند. در بررسی مولکولی نیز در هنگام تنش سرما، مسیرهای بیوسنتز فعال شامل بیوسنتز آرژنین، اسیدهای آمینه پرولین، فنیل پروپانوئید، ریبوفلاوین و مسیر پیام‌رسان mapk شناسایی شدند که در هنگام تنش سرما، بیان ژن‌های تنظیم‌کننده آن‌ها، در رقم مقاوم افزایش بیان بیشتری نسبت به رقم حساس داشت. همچنین ژن pr1 در مسیر پیام‌رسان mapk شناسایی گردید که با کاهش دما، در رقم مقاوم افزایش بیان بیشتری نسبت به رقم حساس از خود نشان داد.

Identification of pathways and genes with differential expression in sugarcane after cold stress, through RNA-seq analysis

IntroductionIn order to identify cold tolerant cultivars in sugarcane, in December 2015 and one week after the occurrence of 1.2 °C, by morphological study of 454 sugarcane cultivars, the reaction of these cultivars and their tolerance to cold were investigated.Materials and methodsIn the first stage of the experiment (morphology), 54 cultivars were introduced as resistant and tolerant cultivars and 400 cultivars were introduced as cold and very sensitive cultivars. In 2016, to evaluate cultivars based on biochemical indices, out of 54 resistant cultivars, 5 cultivars, and out of 400 susceptible cultivars, 5 cultivars were selected and the amount of free amino acids, Proline and Malondialdehyde, in them before and after the onset of cold, it was measured. Among the selected cultivars of the previous stage, BR0001 cultivar was selected as the most tolerant cultivar, and TUC66107 cultivar as the most sensitive cultivar to cold, to apply cold stress in the cold storage and continue the project in the molecular stage. After planting them in plastic pots, leaf sampling was performed in the refrigeration under stress and to extract TotalRNA. After RNA extraction, the quantity and quality of RNA samples were measured. Using the Illumina Hiseq250 platform in the RNAseq technique, after confirming the quantity and quality of the extracted RNAs, the TruSeq cDNA library was constructed by the Chinese company Novogene and sequencing for RNAseq was done as a pairend with a length of 150 nucleotides. After sequencing, the data were downloaded in compressed files in Fastq format from the Novogene website. The quality of the raw reads was checked using FastQC software. Trimmomatic software was used in the trimming step. After aligning the reads to the sugarcane reference genome, quantified reads were normalized, differential expression analysis was performed, and genes with differential expression were analyzed and pathways were mapped. After analyzing the differential expression of genes and obtaining the values, GO enrichment analysis was performed to classify the genes based on the placement of their products. Next, the path analysis was performed using the kobas database. Then, in order to enrich the pathways, significant pathways were identified with Corrected PValue <0.05 for genes with increased and decreased expression.Results and discussionhe results of differential expression of genes showed that out of a total of 62285 expressed genes, 12674 genes have differential expression in different comparisons between treatments, of which 6939 genes decreased expression and 5735 genes increased expression. After enriching the pathways, genes that could be 95% probable in 18 treatment comparisons to attribute their change in expression to cold stress were identified, and marker genes were identified in these comparisons. According to the results of the metabolic study, increasing the concentration of free amino acids, Proline and Malondialdehyde increased the plant’s tolerance to cold. In this study, due to cold stress, various biosynthetic and metabolic pathways were activated under the influence of the activity of genes controlling these pathways to produce their products. These pathways include the biosynthesis of the amino acids arginine and Proline, the biosynthesis of Phenylpropanoid, the riboflavin metabolism, and the MAPK messenger pathway, in which the genes involved increased expression during cold stress. Meanwhile, the PR1 gene was detected in the MAPK messenger pathway, which showed an increase in expression with decreasing temperature Increased expression of genes in MAPK, Phenylpropanoid, sulfur, Glycerophospholipid, arachidonic acid, lipid ether, riboflavin and Proline signaling pathways indicates that the products of some genes controlling these pathways play a role in increasing cold tolerance of sugarcane.ConclusionsThe results of morphological, metabolic and molecular studies showed that there is a correlation between these indices and their values in resistant cultivar. When screening sugarcane clones for cold, the best clones can be selected by measuring the amount of secondary metabolites. In this study, genes affecting cell physiology pathways were identified and their role in increasing sugarcane tolerance to cold was elucidated. Using these genes, you can use them in a sugarcane breeding program to screen for coldresistant clones. In the sugarcane breeding program, these genes can be used as markers in the screening of clones.