DOI: https://doi.org/10.1038/s41467-025-61527-1
PMID: https://pubmed.ncbi.nlm.nih.gov/40624070
تاريخ النشر: 2025-07-07
المؤلف: Han Cheng وآخرون
الموضوع الرئيسي: التvariations الكروموسومية والجينية
طرق
في هذه الدراسة، تم استخدام شجرة المطاط المستنسخة CATAS 7-33-97 لإجراء تحليلات شاملة للترانسكريبتوم والميتابولوم. تم جمع أوراق برونزية شابة وعينات من اللاتكس الطازج من الأشجار التي تم تقطيعها باستخدام طريقة التقطيع الحلزوني نصف، مع تجميد العينات بسرعة في النيتروجين السائل للتحليل اللاحق. تم إنتاج بيانات تسلسل الجينوم الكامل (WGS) باستخدام منصات متعددة، بما في ذلك PacBio Sequel II لقراءات عالية الدقة، وOxford Nanopore لتسلسلات ONT، وIllumina لتسلسلات Hi-C. تم إجراء تسلسل مناعة الكروماتين (ChIP-Seq) على نوى مرتبطة بالفورمالديهايد من الأوراق الشابة، مع توفير بروتوكولات مفصلة في المواد التكميلية.
أشارت النتائج إلى تغييرات ملحوظة في التعبير على مر الزمن، تم تحليلها من خلال تجميع K-means، مع تحديد فترات التقطيع الرئيسية على أنها T1 وT7 وTN. كشفت تحليل إثراء الوظائف عن مجموعات جينية متميزة، وأبرزت الشبكة التنظيمية عوامل النسخ الرئيسية المشاركة في تخليق الجزيئات الكبيرة. من الجدير بالذكر أن تطبيق حمض الجاسمونيك (JA) أظهر أنه يعزز نشاط تخليق المطاط وعائد المطاط الجاف خلال مراحل التقطيع المبكرة. تم استخدام تحليلات إحصائية، بما في ذلك اختبارات هايبرجيومترية أحادية الجانب واختبارات t لطلاب ذات جانبين، للتحقق من النتائج، مع تقديم البيانات كقيم متوسطة ± SD. كما شملت الدراسة اختبارات ثنائية اللمعان واختبارات هجين الخميرة لتوضيح التفاعل بين عوامل النسخ ومناطق المحفز المتعلقة بتخليق المطاط.
نتائج
تقدم البحث تجميعًا شاملاً للجينوم شجرة المطاط (CATAS 7-33-97) مع تحديد الهياكل الجينية، باستخدام 193.19 جيجابايت من قراءات HiFi، و654.04 جيجابايت من قراءات Hi-C، و199.24 جيجابايت من قراءات ONT فائقة الطول. أسفر التجميع الأولي عن 256 قطعة إجمالية totaling 3.16 جيجابايت، والتي تم تنقيحها لإنتاج نوعين من الهياكل الجينية بدون فجوات، hapA وhapB، كل منهما بحجم حوالي 1.56 جيجابايت. احتوى الهيكل الجيني A على 38,635 جينًا مشفرًا للبروتين، بينما كان لدى hapB 38,618، مع ظهور كلا الهيكلين الجينيين بدرجة عالية من التسلسلات المتكررة (84.63% لـ hapA و84.53% لـ hapB). تم تأكيد جودة التجميع مع اكتمال بنسبة 97.5% ودرجات BUSCO تبلغ 0.988 لـ hapA و0.99 لـ hapB، إلى جانب معدلات خطأ قاعدة منخفضة.
بالإضافة إلى ذلك، حددت الدراسة المناطق المركزية من خلال تحليل ChIP-seq، كاشفة عن قمم ربط متميزة والعديد من تسلسلات التكرار المتتالي، بينما تم تمييز التيلوميرات بواسطة تسلسل تكرار مكون من 7 قواعد. من الجدير بالذكر أنه تم توثيق تنوع جيني واسع، بما في ذلك 2,476,725 SNPs و6330 InDels، مع تحديد متغير هيكلي كبير (sv33M) على الكروموسوم 8، مما يشير إلى تنوع جيني كبير بين سكان شجرة المطاط المزروعة والبرية. قد تلعب هذه الاختلافات الهيكلية الكبيرة دورًا حاسمًا في الديناميات التطورية لـ Hevea brasiliensis.
مناقشة
تسلط قسم المناقشة في ورقة البحث الضوء على النتائج المهمة من الجينوميات المقارنة للجينات المنتجة للمطاط ضمن عائلة Euphorbiaceae، مع التركيز بشكل خاص على جينوم T2T لشجرة المطاط (Hevea brasiliensis). حددت الدراسة 27,311 مجموعة أورثوغرافية عبر ثمانية أنواع، مع 8,364 مجموعة أورثوغرافية مشتركة بين الجميع، مما يشير إلى توسعات وانكماشات واسعة في عائلات الجينات تتوافق مع إنتاج المطاط. من الجدير بالذكر أن المستنسخ المزروع من المطاط CATAS 7-33-97 أظهر 1,297 مجموعة أورثوغرافية موسعة و662 مجموعة متقلصة مقارنة بجينوم المطاط البري ذو العائد المنخفض، مع كشف تحليلات إثراء علم الجينات (GO) أن العائلات الموسعة تشارك في عمليات حاسمة لإنتاج المطاط، مثل تخليق التيربين واستجابة الجروح.
توضح الدراسة أيضًا دور عائلات الجينات المحددة، وخاصة عامل إطالة المطاط (REF) وجينات بروتين الجسيمات المطاطية الصغيرة (SRPP)، في تعزيز قدرات إنتاج المطاط. وجدت الدراسة زيادة كبيرة في عدد جينات REF/SRPP في أشجار المطاط المزروعة مقارنة بالأقارب البرية، مما يشير إلى أساس جيني لتحسين عائد اللاتكس. بالإضافة إلى ذلك، كشفت تحليل التقطيع المتتالي أن مسار الميفالونات يلعب دورًا حاسمًا في تخليق المطاط، مع تراكم الميفالونات كخزان رئيسي للكربون. تقترح الدراسة نموذجًا حيث ينظم إشارات حمض الجاسمونيك (JA)، الذي يتم تفعيله بواسطة الجروح الميكانيكية أثناء التقطيع، تخليق المطاط من خلال زيادة التعبير عن الجينات الرئيسية، بما في ذلك MVK1، مما يعزز عائد اللاتكس. بشكل عام، توفر هذه الدراسة فهمًا شاملاً للشبكات الجينية والهرمونية التنظيمية التي تدعم إنتاج المطاط، مما يمهد الطريق لاستراتيجيات التربية المستقبلية التي تهدف إلى تطوير مستنسخات مطاطية عالية العائد.
DOI: https://doi.org/10.1038/s41467-025-61527-1
PMID: https://pubmed.ncbi.nlm.nih.gov/40624070
Publication Date: 2025-07-07
Author(s): Han Cheng et al.
Primary Topic: Chromosomal and Genetic Variations
Methods
In this study, the rubber tree clone CATAS 7-33-97 was utilized for comprehensive transcriptome and metabolome analyses. Young bronze leaves and fresh latex samples were collected from trees tapped using a half spiral tapping method, with samples flash-frozen in liquid nitrogen for subsequent analysis. The whole genome sequencing (WGS) data were generated using multiple platforms, including PacBio Sequel II for high-fidelity reads, Oxford Nanopore for ONT sequences, and Illumina for Hi-C sequences. Chromatin immunoprecipitation sequencing (ChIP-Seq) was conducted on formaldehyde-crosslinked nuclei from young leaves, with detailed protocols provided in supplementary materials.
The results indicated significant expression changes over time, analyzed through K-means clustering, with key tapping intervals denoted as T1, T7, and TN. Functional enrichment analysis revealed distinct gene clusters, and the regulatory network highlighted key transcription factors involved in macromolecular synthesis. Notably, the application of jasmonic acid (JA) was shown to enhance rubber biosynthesis activity and dry rubber yield during early tapping stages. Statistical analyses, including one-sided hypergeometric tests and two-sided Student’s t-tests, were employed to validate findings, with data presented as mean values ± SD. The study also included dual-luciferase assays and yeast one-hybrid assays to elucidate the interaction between transcription factors and promoter regions related to rubber biosynthesis.
Results
The research presents a comprehensive haplotype-resolved assembly of the rubber tree genome (CATAS 7-33-97), utilizing 193.19 Gb of HiFi reads, 654.04 Gb of Hi-C reads, and 199.24 Gb of ultra-long ONT reads. The initial assembly yielded 256 contigs totaling 3.16 Gb, which were refined to produce two gap-free haplotypes, hapA and hapB, each approximately 1.56 Gb in size. Haplotype A contained 38,635 protein-coding genes, while hapB had 38,618, with both haplotypes exhibiting a high degree of repetitive sequences (84.63% for hapA and 84.53% for hapB). The assembly quality was confirmed with a completeness of 97.5% and BUSCO scores of 0.988 for hapA and 0.99 for hapB, alongside low base error rates.
Additionally, the study identified centromeric regions through ChIP-seq analysis, revealing distinct binding peaks and multiple tandem repeat sequences, while telomeres were characterized by a 7-base repeat sequence. Notably, extensive genomic variations were documented, including 2,476,725 SNPs and 6330 InDels, with a significant structural variant (sv33M) identified on chromosome 8, suggesting substantial genetic diversity between the cultivated and wild rubber tree populations. This large structural variation may play a critical role in the evolutionary dynamics of Hevea brasiliensis.
Discussion
The discussion section of the research paper highlights significant findings from the comparative genomics of rubber-producing genes within the Euphorbiaceae family, particularly focusing on the T2T genome of the rubber tree (Hevea brasiliensis). The study identified 27,311 orthogroups across eight species, with 8,364 orthogroups common to all, indicating extensive gene family expansions and contractions that correlate with rubber yield. Notably, the cultivated rubber clone CATAS 7-33-97 exhibited 1,297 expanded and 662 contracted orthogroups compared to a low-yield wild rubber genome, with gene ontology (GO) enrichment analyses revealing that expanded families are involved in processes critical for rubber production, such as terpene biosynthesis and wounding response.
The research further elucidates the role of specific gene families, particularly the rubber elongation factor (REF) and small rubber particle protein (SRPP) genes, in enhancing rubber production capabilities. The study found a significant increase in the number of REF/SRPP genes in cultivated rubber trees compared to wild relatives, suggesting a genetic basis for improved latex yield. Additionally, the analysis of consecutive tapping revealed that the mevalonate pathway plays a crucial role in rubber biosynthesis, with mevalonate accumulation serving as a major carbon reservoir. The study proposes a model where jasmonic acid (JA) signaling, activated by mechanical wounding during tapping, regulates rubber biosynthesis by upregulating key genes, including MVK1, thereby enhancing latex yield. Overall, this research provides a comprehensive understanding of the genetic and hormonal regulatory networks that underpin rubber production, paving the way for future breeding strategies aimed at developing high-yield rubber clones.