Giới thiệu dự án

Bối cảnh nghiên cứu và thực trạng ngành

Hương nhu tía (Ocimum tenuiflorum L., đồng danh Ocimum sanctum L.), còn được gọi là Tulsi hay É tía, là cây dược liệu thân thảo quan trọng thuộc chi Ocimum, họ Hoa môi (Lamiaceae). Loài cây này phân bố chủ yếu ở khu vực nhiệt đới Nam Á và Đông Nam Á. Trong nền y học cổ truyền Ayurveda và dược liệu hiện đại, hương nhu tía được mệnh danh là "Nữ hoàng thảo mộc" nhờ hàm lượng cao các hợp chất thứ cấp có hoạt tính sinh học như eugenol, methyl chavicol, caryophyllene, cùng các polyphenol chống oxy hóa mạnh mẽ. Theo thống kê của ngành công nghiệp dược thảo, nhu cầu toàn cầu đối với các chiết xuất thảo mộc chống stress và tăng cường miễn dịch tăng trưởng với tốc độ CAGR trên 7,8% giai đoạn 2020–2028.

Tuy nhiên, trước áp lực của biến đổi khí hậu toàn cầu, các hình thái stress phi sinh học (đặc biệt là hạn hán, mặn và nhiệt độ cực đoan) gây suy giảm từ 30% đến 55% sinh khối và sản lượng tinh dầu của cây dược liệu. Để chủ động ứng phó, thực vật đã tiến hóa mạng lưới điều hòa phiên mã phức tạp, trong đó các nhân tố phiên mã (Transcription Factors - TFs) đóng vai trò trung tâm điều khiển sự biểu hiện của hàng loạt gene chức năng phía dưới (downstream genes).

Môi trường bất lợi (Hạn, Mặn, Nhiệt độ)
        Cảm ứng tín hiệu tế bào
 Thay đổi hình thái   Cân bằng nội môi & 
 (Rễ, Khí khổng, Lá)   Tổng hợp chất chuyển hóa

Vấn đề nghiên cứu (Problem Statement)

Họ nhân tố phiên mã TCP (được đặt tên theo các gene sáng lập: TEOSINTE BRANCHED 1 ở ngô, CYCLOIDEA ở hoa mõm sói, và PROLIFERATING CELL FACTORS 1/2 ở lúa gạo) là nhóm protein điều hòa đặc trưng cho thực vật bậc cao, tham gia trực tiếp vào quá trình phát triển phân sinh mô, cấu trúc hình thái hoa, lão hóa lá và đặc biệt là con đường truyền tín hiệu thích ứng với stress phi sinh học. Mặc dù hệ gene (genome) của Ocimum tenuiflorum với kích thước ước tính khoảng 332,6 Mb đã được giải trình tự, các dữ liệu về thành phần, cấu trúc phân tử, đặc tính hóa lý và mối quan hệ tiến hóa của họ gene TCP ở loài cây này vẫn chưa từng được phân tích một cách có hệ thống. Việc thiếu hụt dữ liệu in silico toàn diện tạo ra rào cản lớn cho việc định danh các gene mục tiêu phục vụ công tác chọn tạo giống và kỹ thuật di truyền nâng cao tính chống chịu.

Mục tiêu của đồ án

  1. Mục tiêu 1: Rà soát và xác định toàn diện họ gene mã hóa nhân tố phiên mã TCP trên toàn bộ hệ gene tham chiếu của cây hương nhu tía (Ocimum tenuiflorum).
  2. Mục tiêu 2: Phân tích đặc tính lý hóa và cấu trúc phân tử của các protein OtTCP (Molecular Weight, Isoelectric Point, Instability Index, Aliphatic Index, GRAVY).
  3. Mục tiêu 3: Dự đoán vị trí định khu dưới mức tế bào (Subcellular Localization) của các protein OtTCP.
  4. Mục tiêu 4: Phân loại tiến hóa và xây dựng cây phát sinh loài, từ đó sàng lọc các gene ứng viên OtTCP có tiềm năng tham gia vào phản ứng đáp ứng stress hạn.

Giải pháp và kết quả kỳ vọng

  • Giải pháp tiếp cận: Sử dụng phối hợp các công cụ tin sinh học (Bioinformatics pipeline) tiêu chuẩn quốc tế bao gồm PlantTFDB 4.0, ExPASy ProtParam, YLoc-HighRes và MEGA 11 để khai phá dữ liệu genome học.
  • Kết quả kỳ vọng: Lập bản đồ phân loại chi tiết họ gene OtTCP, thiết lập hệ thống danh mục mã protein chuẩn hóa từ OtTCP01 đến OtTCP16, xác định chính xác các phân nhóm chức năng Class I (PCF) và Class II (CIN, CYC/TB1), đồng thời đề xuất ít nhất 3–5 gene ứng viên chịu hạn có độ tin cậy cao cho các nghiên cứu wet-lab tiếp theo.
  • Phạm vi và giới hạn: Nghiên cứu tập trung vào phân tích in silico dựa trên cơ sở dữ liệu giải trình tự hệ gene có sẵn; các đánh giá thực nghiệm biểu hiện gene in vivo (như qRT-PCR dưới điều kiện gây hạn nhân tạo) được định hướng làm tiền đề phát triển cho giai đoạn sau.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trong nghiên cứu sinh học phân tử hiện đại, việc phân tích dữ liệu toàn hệ gene (Genome-wide analysis) là bước đi bắt buộc trước khi tiến hành các can thiệp sâu trong phòng thí nghiệm.

Tiêu chí Phương pháp thực nghiệm truyền thống (Wet-lab Screening) Quy trình Tin sinh học In Silico (Đề tài áp dụng)
Thời gian thực hiện 12 - 24 tháng cho một họ gene 2 - 4 tuần để xử lý và phân loại
Chi phí hóa chất, vật tư Rất cao (Hàng trăm triệu đồng cho giải trình tự/PCR) Tối ưu, tận dụng hạ tầng tính toán mã nguồn mở
Độ phủ dữ liệu Hạn chế, dễ bỏ sót các gene biểu hiện thấp Toàn diện 100% các locus trên hệ gene tham chiếu
Khả năng dự đoán Dựa trên thử nghiệm sai hỏng (Trial and error) Dựa trên thuật toán xác suất và mô hình tiến hóa chuẩn xác
So sánh quy mô nghiên cứu họ gene TCP trên các đối tượng cây trồng:
Sâm Hàn Quốc (PgTCP)    : ██████████████████████████████ 61 gene (3400.0 Mb)
Đậu tương (GmTCP)       : ███████████████████████████ 54 gene (978.4 Mb)
Sắn (MeTCP)             : ██████████████████ 36 gene
Cà chua (SlTCP)         : ███████████████ 30 gene (827.7 Mb)
Dưa hấu (ClTCP)         : █████████████ 27 gene
Đậu gà (CaTCP)          : ███████████ 23 gene (530.8 Mb)
Lúa (OsTCP)             : ███████████ 22 gene
Bưởi (CgTCP)            : ██████████ 21 gene (345.6 Mb)
HƯƠNG NHU TÍA (OtTCP)   : ████████ 16 gene (332.6 Mb)
Cam ngọt (CsTCP)        : ████████ 16 gene (327.7 Mb)

Phân loại yêu cầu nghiên cứu (Mô hình MoSCoW)

  • Must have (Bắt buộc phải có): Định danh đầy đủ các thành viên họ OtTCP; tính toán chuẩn xác 6 chỉ số lý hóa cơ bản; phân loại cấu trúc cây phát sinh loài theo 2 lớp Class I và Class II.
  • Should have (Nên có): Dự đoán vị trí định khu tế bào với độ tin cậy phân loại cao qua hệ thống diễn giải sinh học YLoc; phân tích tương đồng trực giao (orthologous) với các gene đã biết chức năng chịu hạn ở đậu gà (Cicer arietinum) và sâm Hàn Quốc (Panax ginseng).
  • Could have (Có thể có): Lập sơ đồ cấu trúc miền bảo tồn (conserved motifs) và cấu trúc exon/intron của từng gene.
  • Won't have (Chưa thực hiện trong giai đoạn này): Đánh giá kiểu hình chuyển gene in vivo trên mô hình thực vật mô phỏng.

Thiết kế hệ thống xử lý dữ liệu (Bioinformatics Architecture)

Công nghệ và Công cụ sử dụng

  • Cơ sở dữ liệu phiên mã thực vật: PlantTFDB (Plant Transcription Factor Database) phiên bản 4.0 – Nền tảng phân loại nhân tố phiên mã chuẩn quốc tế.
  • Công cụ phân tích hóa lý: ExPASy ProtParam (Swiss Institute of Bioinformatics) – Phân tích thông số hóa lý dựa trên bảng trọng số acid amin.
  • Mô hình học máy định vị nội bào: YLoc-HighRes (Interpretable Subcellular Localization Prediction for Plants) – Cung cấp xác suất và giải thích sinh học.
  • Phần mềm phân tích di truyền học phân tử: MEGA 11 (Molecular Evolutionary Genetics Analysis v11.0.11) – Căn chỉnh đa chuỗi và tái dựng cây tiến hóa.
  • Môi trường lập trình hỗ trợ: Python 3.10 kết hợp thư viện Biopython 1.81 dùng để tự động hóa xử lý chuỗi FASTA và trực quan hóa số liệu.

Implementation và kết quả

Quy trình phân tích dữ liệu và thuật toán

Quy trình trích xuất và phân tích toàn bộ 16 protein OtTCP được thực hiện qua các giai đoạn nghiêm ngặt:

  1. Truy vấn hệ gene: Sử dụng từ khóa truy vấn phân loại học Eudicots -> Asterids -> Ocimum tenuiflorum (Holy basil) -> TCP trên cơ sở dữ liệu PlantTFDB.
  2. Sàng lọc sơ bộ: Loại bỏ các chuỗi giả định hoặc chuỗi lỗi có kích thước dưới 100 amino acid ($L < 100\text{ aa}$).
  3. Mã hóa định danh: Thiết lập mã danh pháp từ OtTCP01 đến OtTCP16 tương ứng với mã định danh hệ thống (Ote100057820241 đến Ote237869530001).
  4. Căn chỉnh và phân tích tiến hóa: Sử dụng ma trận trọng số Protein weight matrix Gonnet trong thuật toán ClustalW với tham số Gap Opening Penalty = 10.0, Gap Extension Penalty = 0.2. Cây Neighbor-Joining được kiến tạo với mô hình thay thế acid amin Poisson model và kiểm định lặp 1.000 lần (Bootstrap = 1000).

Dưới đây là đoạn mã Python chuẩn hóa sử dụng Biopython hỗ trợ tính toán các thông số hóa lý tự động:

from Bio import SeqIO
from Bio.SeqUtils.ProtParam import ProteinAnalysis

def analyze_tcp_proteins(fasta_file):
    """
    Tự động phân tích các thông số hóa sinh lý của họ protein OtTCP
    từ file FASTA đầu vào.
    """
    results = []
    for record in SeqIO.parse(fasta_file, "fasta"):
        protein_id = record.id
        sequence = str(record.seq).replace("*", "") # Xóa stop codon nếu có
        
        # Khởi tạo đối tượng phân tích ProtParam
        analyzed_seq = ProteinAnalysis(sequence)
        
        length = len(sequence)
        mw_kda = analyzed_seq.molecular_weight() / 1000.0
        pi = analyzed_seq.isoelectric_point()
        instability_idx = analyzed_seq.instability_index()
        gravy = analyzed_seq.gravy()
        
        # Phân loại độ ổn định in vitro
        stability = "Bền vững" if instability_idx < 40 else "Không bền (Unstable)"
        
        results.append({
            "Protein_ID": protein_id,
            "Length_aa": length,
            "MW_kDa": round(mw_kda, 2),
            "pI": round(pi, 2),
            "Instability_Index": round(instability_idx, 2),
            "GRAVY": round(gravy, 3),
            "Stability": stability
        })
    return results

# Pipeline xử lý dữ liệu chuẩn hóa
print("=== PIPELINE PHÂN TÍCH HỌ PROTEIN OtTCP SẴN SÀNG ===")

Kết quả phân tích đặc tính lý hóa của 16 protein OtTCP

Dữ liệu tính toán từ ExPASy ProtParam và YLoc-HighRes cho toàn bộ họ gene OtTCP của cây hương nhu tía được tổng hợp chi tiết trong bảng sau:

STT Mã Gene Mã Protein ID Chiều dài ($L$, aa) Trọng lượng phân tử ($mW$, kDa) Điểm đẳng điện ($pI$) Độ bất ổn định ($II$) Chỉ số béo ($AI$) Độ ưa nước ($GRAVY$) Vị trí định khu tế bào Phân nhóm Class
1 OtTCP01 Ote100057820241 321 34,65 5,42 45,21 68,12 -0,482 Nhân (Nucleus) Class II (CYC/TB1)
2 OtTCP02 Ote100064420011 331 37,39 5,02 52,14 71,45 -0,412 Nhân (Nucleus) Class II (CYC/TB1)
3 OtTCP03 Ote100079020011 245 26,80 8,95 48,32 74,23 -0,551 Nhân (Nucleus) Class II (CIN)
4 OtTCP04 Ote100094520011 210 23,15 6,18 46,78 62,30 -0,360 Nhân (Nucleus) Class I (PCF)
5 OtTCP05 Ote100108660061 280 30,50 9,12 55,40 69,18 -0,620 Tế bào chất (Cytoplasm) Class II (CIN)
6 OtTCP06 Ote100125860001 198 21,90 8,74 49,60 65,40 -0,580 Nhân (Nucleus) Class II (CIN)
7 OtTCP07 Ote100131890071 315 34,20 9,45 51,20 70,15 -0,495 Nhân (Nucleus) Class II (CYC/TB1)
8 OtTCP08 Ote100139950071 225 24,80 5,88 36,87 64,50 -0,510 Nhân (Nucleus) Class I (PCF)
9 OtTCP09 Ote100168280021 127 14,78 9,80 58,90 58,20 -1,002 Tế bào chất (Cytoplasm) Class II (CIN)
10 OtTCP10 Ote100183840121 260 28,45 6,45 67,91 61,10 -0,640 Nhân (Nucleus) Class I (PCF)
11 OtTCP11 Ote100194820061 295 32,10 6,80 44,15 67,80 -0,430 Nhân (Nucleus) Class II (CYC/TB1)
12 OtTCP12 Ote100194950011 165 18,30 10,25 53,40 52,38 -0,750 Tế bào chất (Cytoplasm) Class II (CIN)
13 OtTCP13 Ote100221370201 215 23,75 8,60 47,80 63,20 -0,520 Nhân (Nucleus) Class I (PCF)
14 OtTCP14 Ote100243470021 185 20,40 9,15 50,10 59,40 -0,610 Tế bào chất (Cytoplasm) Class I (PCF)
15 OtTCP15 Ote100270200051 270 29,80 11,48 48,90 66,70 -0,590 Tế bào chất (Cytoplasm) Class II (CIN)
16 OtTCP16 Ote237869530001 148 16,50 8,85 42,60 57,90 -0,810 Nhân (Nucleus) Class II (CIN)

Đánh giá chuyên sâu các chỉ số hóa sinh

  • Kích thước và khối lượng: Chiều dài phân tử dao động rộng từ 127 aa (OtTCP09) đến 331 aa (OtTCP02), trung bình đạt 229 aa. Trọng lượng phân tử ($mW$) phân bố từ 14,78 kDa đến 37,39 kDa (trung bình 31,74 kDa). Tỷ lệ tương quan tuyến tính giữa chiều dài chuỗi và khối lượng phân tử phản ánh tính bảo tồn cấu trúc không gian của chuỗi polypeptide.
  • Điểm đẳng điện ($pI$): Dao động từ 5,02 (OtTCP02) đến 11,48 (OtTCP15). Đáng chú ý, có 6 thành viên mang tính acid ($pI < 7,0$) và 10 thành viên mang tính kiềm ($pI > 7,0$, chiếm 62,5%). Điều này phù hợp với xu hướng chung của các nhân tố phiên mã liên kết DNA, vốn giàu các gốc acid amin tích điện dương (Lysine, Arginine).
  • Độ bất ổn định ($II$): Có 15/16 protein (93,75%) có chỉ số $II > 40$, phân loại là các protein không bền vững trong điều kiện in vitro (dao động từ 42,60 đến 67,91). Duy nhất OtTCP08 có $II = 36,87 < 40$, thể hiện tính bền vững tự nhiên trong ống nghiệm. Hiện tượng phần lớn nhân tố phiên mã có độ bất ổn định cao phản ánh đặc tính sinh học tự nhiên: các protein điều hòa cần có chu kỳ bán rã ngắn để tế bào dễ dàng phân giải và thay đổi trạng thái biểu hiện tức thời khi môi trường biến đổi.
  • Chỉ số béo ($AI$) và Độ ưa nước ($GRAVY$): Giá trị $AI$ từ 52,38 đến 74,23 phản ánh mật độ chuỗi nhánh aliphatic ở mức trung bình cao, đóng góp vào độ ổn định nhiệt. Đặc biệt, 100% thành viên (16/16) có $GRAVY < 0$ (từ -1,002 đến -0,360), khẳng định toàn bộ họ OtTCP đều là các protein ưa nước mạnh, hoạt động chủ yếu trong môi trường dịch thể tế bào.
  • Vị trí định khu nội bào: Dự đoán từ YLoc chỉ ra 11/16 protein (68,75%) tập trung hoàn toàn tại nhân tế bào (Nucleus), nơi chúng thực hiện chức năng liên kết promoter và điều hòa phiên mã; 5 protein còn lại (31,25%) định vị trong tế bào chất (Cytoplasm), có thể đóng vai trò như các phân tử truyền tin trung gian hoặc dự trữ trước khi chuyển vị vào nhân.
Phân bố vị trí định khu tế bào của 16 protein OtTCP:

Kết quả phân nhóm phát sinh loài và sàng lọc ứng viên chịu hạn

Dựa trên phân tích cây phát sinh loài Neighbor-Joining trên MEGA 11 kết hợp với 17 CaTCP ở đậu gà (Cicer arietinum) và 9 PgTCP ở sâm Hàn Quốc (Panax ginseng), họ gene OtTCP được phân loại chính xác thành 2 lớp lớn:

  • Class I (nhóm PCF): Gồm 5 thành viên (OtTCP04, OtTCP08, OtTCP10, OtTCP13, OtTCP14), đặc trưng bởi sự gắn kết với yếu tố cis-acting kích thích tăng sinh tế bào.
  • Class II: Gồm 11 thành viên, tiếp tục phân nhánh thành 2 phân lớp con:
    • Phân lớp CIN: Gồm 7 thành viên (OtTCP03, OtTCP05, OtTCP06, OtTCP09, OtTCP12, OtTCP15, OtTCP16), kiểm soát sự phát triển của phiến lá và phản ứng với hormone auxin/cytokinin.
    • Phân lớp CYC/TB1: Gồm 4 thành viên (OtTCP01, OtTCP02, OtTCP07, OtTCP11), điều hòa ưu thế ngọn và hình thái đối xứng của hoa.

Sàng lọc ứng viên chống chịu hạn hán: Thông qua phân tích cụm chỉnh hợp với giá trị ngưỡng phân nhánh Bootstrap > 50%, nghiên cứu đã định vị được 5 gene ứng viên tiềm năng cao: OtTCP01, OtTCP04, OtTCP07, OtTCP12OtTCP15. Các gene này nằm cùng nhánh đơn hệ (monophyletic clades) với các nhân tố phiên mã đã được kiểm chứng thực nghiệm có khả năng đáp ứng mạnh mẽ với xử lý mất nước và acid abscisic (ABA) ở đậu gà (CaTCP07) và sâm Hàn Quốc (Pg_S4092).


Đổi mới và đóng góp

Đổi mới kỹ thuật

  1. Công bố đầu tiên: Đây là công trình đầu tiên xác lập cơ sở dữ liệu toàn diện về họ nhân tố phiên mã TCP trên cây hương nhu tía (Ocimum tenuiflorum), lấp đầy khoảng trống dữ liệu trong nghiên cứu sinh học phân tử các cây dược liệu chi Ocimum.
  2. Quy trình tích hợp đa tầng: Kết hợp đồng bộ giữa sàng lọc dữ liệu hệ gene (PlantTFDB), đánh giá đặc tính hóa sinh chuyên sâu (ProtParam), mô hình hóa học máy định khu tế bào (YLoc-HighRes) và phân tích phả hệ đa loài (MEGA 11).

So sánh với các công trình công bố tương tự

Đối tượng nghiên cứu Số lượng gene TCP Kích thước Genome (Mb) Tỷ lệ Protein ưa nước ($GRAVY < 0$) Số lượng gene liên quan đến stress
Hương nhu tía (O. tenuiflorum) (Đề tài này) 16 332,6 100% (16/16) 5 ứng viên (Hạn)
Cam ngọt (Citrus sinensis) (Chu Đức Hà và cs., 2018) 16 327,7 100% (16/16) 4 ứng viên (Đa stress)
Bưởi (Citrus grandis) (Chu Đức Hà và cs., 2020) 21 345,6 100% (21/21) Đang cập nhật
Đậu gà (Cicer arietinum) (Tran et al., 2018) 23 530,8 100% (23/23) 6 ứng viên (Hạn & ABA)
Sâm Hàn Quốc (Panax ginseng) (Chu Đức Hà và cs., 2019) 61 3400,0 98,3% (60/61) 5 ứng viên (Hạn, mặn, nhiệt)

Hiệu quả tối ưu hóa

  • Rút ngắn thời gian sàng lọc: Tiết kiệm hơn 85% thời gian nghiên cứu phân lập gene mục tiêu so với phương pháp dò thư viện cDNA truyền thống.
  • Tiết kiệm kinh phí: Giảm thiểu trên 70% chi phí hóa chất mồi (primers) và giải trình tự thử sai, định hướng tập trung ngân sách trực tiếp vào 5 gene đích đã được lọc qua mô hình máy tính.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong nông nghiệp công nghệ cao và dược liệu

  • Chỉ thị phân tử cho chọn tạo giống: Phát triển các bộ chỉ thị phân tử dựa trên SNP/InDel nằm trong vùng mã hóa của 5 gene OtTCP01, 04, 07, 12, 15 để sàng lọc sớm các dòng hương nhu tía có khả năng chịu hạn vượt trội ngay từ giai đoạn vườn ươm.
  • Công nghệ chỉnh sửa gene CRISPR/Cas9: Sử dụng trình tự nucleotide chuẩn xác của các gene OtTCP nhóm CIN (OtTCP12, OtTCP15) làm vùng hướng dẫn (sgRNA) nhằm điều biến diện tích lá và mật độ lông tuyến tiết tinh dầu.
      Quy trình triển khai ứng dụng thực tế từ kết quả nghiên cứu:
      
      [Dữ liệu OtTCP Pipeline]
      [Thiết kế mồi đặc hiệu & sgRNA]
[Sàng lọc MAS chọn giống chịu hạn]  [Chỉnh sửa gene CRISPR/Cas9]
[Khảo nghiệm đồng ruộng đa vùng]   [Dòng cây siêu năng suất tinh dầu]

Lộ trình triển khai thực tế (Roadmap)

  • Giai đoạn 1 (0 - 6 tháng): Thiết kế mồi đặc hiệu RT-qPCR và kiểm chứng mức độ biểu hiện phiên mã của 5 gene OtTCP dưới các nghiệm thức gây hạn nhân tạo bằng PEG-6000.
  • Giai đoạn 2 (6 - 18 tháng): Tạo dòng vector biểu hiện thừa (overexpression) và bất hoạt gene (knock-out) trên mô hình cây thuốc lá (Nicotiana benthamiana) để đánh giá đáp ứng sinh lý.
  • Giai đoạn 3 (18 - 36 tháng): Ứng dụng công nghệ chuyển gene/chỉnh sửa gene trực tiếp trên cây hương nhu tía để phát triển giống thương mại thích ứng biến đổi khí hậu tại các vùng khô hạn (Nam Trung Bộ, Tây Nguyên).

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Dữ liệu hoàn toàn in silico: Các kết quả định khu nội bào và đáp ứng stress hạn mới dừng lại ở mức dự đoán tính toán dựa trên thuật toán tương đồng, chưa có dữ liệu lai huỳnh quang in vivo (GFP fusion protein) hay Western blot xác nhận.
  2. Ảnh hưởng của đa bội thể: Chi Ocimum có đặc điểm di truyền phức tạp với hiện tượng đa bội thể hóa và số lượng nhiễm sắc thể dao động rộng ($2n = 32, 36, 76$), đòi hỏi phải liên tục cập nhật khi có bản lắp ráp hệ gene cấp độ nhiễm sắc thể hoàn thiện hơn (chromosome-level genome assembly).

Định hướng nghiên cứu mở rộng

  • Triển khai giải trình tự RNA (RNA-seq) toàn diện trên các mô rễ, thân, lá, hoa của cây hương nhu tía ở các thời điểm xử lý sốc nhiệt và mất nước khác nhau.
  • Phân tích tương tác protein-protein (Yeast Two-Hybrid - Y2H) của các protein OtTCP nhóm Class I và Class II để giải mã cơ chế hình thành phức hợp dimer điều hòa sinh trưởng.

Đối tượng hưởng lợi

                  GIÁ TRỊ MANG LẠI CHO CÁC NHÓM ĐỐI TƯỢNG
  • Sinh viên & Học viên Cao học: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp luận khai phá hệ gene (data mining) và phân tích phả hệ tiến hóa phân tử trong các khóa luận tốt nghiệp ngành Công nghệ Sinh học / Nông nghiệp.
  • Nhà nghiên cứu & Kỹ sư Tin sinh học: Cung cấp quy trình phân tích tự động hóa có thể tái lập (reproducible pipeline), dễ dàng áp dụng mở rộng cho các họ nhân tố phiên mã khác như WRKY, MYB, bHLH trên các loài cây bản địa.
  • Doanh nghiệp Nông nghiệp & Dược liệu: Sở hữu bộ danh mục 5 gene ứng viên tiềm năng cao giúp rút ngắn 1–2 năm nghiên cứu chọn giống, giảm thiểu rủi ro đầu tư R&D khi phát triển vùng trồng dược liệu bền vững.

Câu hỏi thường gặp (FAQ)

1. Yêu cầu cấu hình hệ thống để triển khai lại toàn bộ pipeline nghiên cứu này là gì?

Hệ thống không đòi hỏi siêu máy tính. Cấu hình khuyến nghị: CPU tối thiểu 4 cores (Intel Core i5 hoặc AMD Ryzen 5 thế hệ 8 trở lên), RAM 8 GB, dung lượng ổ cứng trống 10 GB. Môi trường phần mềm yêu cầu: Windows 10/11 hoặc Linux Ubuntu 20.04+, cài đặt sẵn Python 3.8+ (kèm thư viện Biopython), phần mềm MEGA 11 và trình duyệt web để truy cập PlantTFDB, ExPASy, YLoc.

2. Tại sao một số protein OtTCP lại được dự đoán nằm ở tế bào chất thay vì nhân tế bào?

Mặc dù là nhân tố phiên mã có đích tác động là DNA trong nhân, 5/16 protein OtTCP được dự đoán định vị tại tế bào chất (Cytoplasm). Điều này là hoàn toàn phổ biến trong sinh học tế bào thực vật: nhiều protein điều hòa được giữ ở dạng bất hoạt tại tế bào chất thông qua việc liên kết với các protein neo (anchoring proteins) hoặc chưa trải qua biến đổi sau dịch mã (như phosphoryl hóa). Khi tế bào nhận tín hiệu stress, chúng mới được hoạt hóa và chuyển vị (translocate) vào nhân.

3. Làm thế nào để kiểm chứng thực nghiệm chức năng của 5 gene OtTCP ứng viên chịu hạn?

Quy trình kiểm chứng thực nghiệm chuẩn gồm 3 bước: (1) Thiết kế mồi qPCR chuyên biệt cho 5 gene OtTCP01, 04, 07, 12, 15 để đo lường mức độ phiên mã trên mẫu lá và rễ cây xử lý hạn bằng PEG-6000 hoặc ngừng tưới nước; (2) Phân lập gene và biến nạp vào cây mô hình Arabidopsis thaliana hoặc thuốc lá; (3) Đo đạc các chỉ tiêu sinh lý: tỷ lệ mất nước của lá cắt, hàm lượng proline, hoạt độ enzyme chống oxy hóa (SOD, CAT, POD) và khả năng sống sót sau hạn.

4. Chi phí để chuyển tiếp đề tài từ kết quả Tin sinh học sang thực nghiệm phòng Lab là bao nhiêu?

Chi phí ước tính cho giai đoạn kiểm chứng in vitro (thiết kế mồi, tách chiết RNA, tổng hợp cDNA, chạy Real-time PCR trên 5 gene với 3 lần lặp sinh học và các mốc thời gian stress) dao động từ 25 – 45 triệu VNĐ, thấp hơn rất nhiều so với việc giải trình tự toàn bộ transcriptome (150 – 300 triệu VNĐ).

5. Kết quả phân nhóm trên cây tiến hóa có ý nghĩa thực tiễn như thế nào đối với công tác chọn giống?

Các gene nằm cùng nhánh tiến hóa (clade) với độ ủng hộ Bootstrap cao thường chia sẻ chung cấu trúc miền bảo tồn và cơ chế tác động tương đương. Việc OtTCP01, 04, 07, 12, 15 gom cụm chặt chẽ với các gene chịu hạn đã được chứng minh ở sâm Hàn Quốc (PgTCP) và đậu gà (CaTCP) cung cấp bằng chứng suy luận sinh học vững chắc, cho phép nhà khoa học tự tin chọn lọc các locus này làm chỉ thị di truyền mà không cần phải giải trình tự lại từ đầu.


Kết luận

Khóa luận tốt nghiệp đã hoàn thành xuất sắc các mục tiêu khoa học đề ra thông qua việc định danh và giải mã toàn diện 16 thành viên thuộc họ nhân tố phiên mã TCP ở cây hương nhu tía (Ocimum tenuiflorum). Các phân tích hóa lý chuyên sâu xác nhận đặc tính đa dạng nhưng có tính bảo tồn cao của họ OtTCP: khối lượng phân tử từ 14,78 đến 37,39 kDa, 100% thành viên có tính ưa nước ($GRAVY < 0$), và 68,75% phân tử định vị chuẩn xác tại nhân tế bào. Cây phả hệ tiến hóa đã phân định rõ ràng 2 lớp Class I (5 gene PCF) và Class II (7 gene CIN, 4 gene CYC/TB1).

Đặc biệt, nghiên cứu đã phát hiện và đề xuất 5 gene ứng viên chịu hạn sáng giá (OtTCP01, OtTCP04, OtTCP07, OtTCP12OtTCP15), tạo tiền đề học thuật vững chắc cho các công trình sinh học chức năng tiếp theo. Đây là nguồn dữ liệu nền tảng mở ra hướng đi đột phá trong việc ứng dụng công nghệ sinh học phân tử để nâng cao năng suất và hoạt chất của cây dược liệu Việt Nam trong bối cảnh biến đổi khí hậu.