Kết quả giải trình trình tự chuỗi gen EBNA-1 của EBV thực hiện trên máy ABI-3100 Avant Genetic Analyzer, được biểu thị bằng chương trình chuyên biệt, thể hiện ở hình 3.11.
Hình 3.11. Hình ảnh giản đồ (chromatogram) một phần trình tự gen EBNA-1 (Ghi chú: (A): EB17, (B): EB18, (C): EB19, (D): EB20)
Hình 3.11 cho thấy, mỗi màu đều xác nhận một loại nucleotide tương ứng,
chứng tỏ chuỗi gen EBNA-1 của EBVthu nhận có chất lượng tốt, được giải trình
trình tự hoàn chỉnh và không bị nhiễu trong quá trình giải trình trình tự. A A C B C C B D D
Bằng các chương trình xử lí chuỗi trình tự nucleotide, chúng tôi thu được
chuỗi nucleotide của gen EBNA-1, sử dụng làm các chuỗi gen đích để truy cập
Ngân hàng gen.
Hình 3.12. Kết quả chuỗi gen EBNA-1 của EB17 được thu nhận sau giải trình tự.
3.2.4. Kết quả truy cập Ngân hàng gen xác định chuỗi gen EBNA-1 của EBV Kết quả xác nhận trình tự nucleotide của các chuỗi gen EBNA-1 của EBV thu nhận sử dụng chương trình Blast thông qua truy cập Ngân hàng gen được thể hiện ở bảng 3.3
Bảng 3.3.Kết quả truy cập Ngân hàng gen sử dụng chuỗi nucleotide gen EBNA-1 của EB17
Số đăng kí Chuỗi gen đích truy cập
Tỉ lệ tương
đồng
AY961628.3 Human herpesvirus 4 strain GD1, complete genome 99% V01555.2 Epstein-Barr virus (EBV) genome, strain B95-8 97% AJ507799.2 Human herpesvirus 4 complete wild type genome 97%
* 20 * 40 * 60 * --//GGGGCATTCAGATAGGAAGCCATTTTTCCACCCTGTAGGGGATGCCGATTATTTTGAATACCTCCAAGAAG 80 * 100 * 120 * 140 GTGGCCCAGATGGTGAGCCTGACGTGCCCCCGGGAGCGATAGAGCAGGGCCCCACAGATGACCCAGGAGA * 160 * 180 * 200 * GGCCCAAGCACTGGACCCCGGGGTCAGGGTGATGGAGGCAGGCGCAAAAAAGGAGGGTGGTTTGGAAAGC 220 * 240 * 260 * 280 ATCGTGGTCAAGGAGGTTCCAACCCGAAATTTGAGAACATTGCAGAAGGTTTAAGAGTTCTCCTGGCTAG * 300 * 320 * 340 * GAGTCACGTAGAAAGGACTACCGAGGAAGGAAATTGGGTCGCCGGTGTGTTCGTATATGGAGGTAGTAAG 360 * 380 * 400 * 420 ACCTCCCTTTACAACCTCAGGCGAGGAATTGCCCTTGCTGTTCCACAATGTCGTATTACACCATTGAGTC * 440 * 460 * 480 * GTCTCCCCTTTGGAATGGCCCCTGGACCCGGCCCACAACCTGGCCCACTAAGGGAGTCCATTGTCTGTTA 500 * 520 * 540 * 560 TTTCATGGTCTTTTTACAAACTCATATATTTGCTGAGGTTTTGAAGGATGCGATTAAGGACCTTGTTAT * 580 * 600 * 620 * GACAAAGCCCGCTCCTACCTGCAATATCAAGGTGACTGTGTGCAGCTTTGACGATGGAGTAGATTTGCCT 640 * 660 * 680 * 700 CCCTGGTTTCCACCTATGGTGGAAGGGGCTGCCGCGGAGGGTGATGACGGAGATGACGGAGATGAAGGAG * 720 GTGATGGAGATGAAGGGTTGAGGAA //--
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
SEQUENCE OF EB17 EBNA1 >
Số đăng kí Chuỗi gen đích truy cập
Tỉ lệ tương
đồng
M80517.1
Epstein-Barr virus, artifactual joining of B95-8 complete (adsbygoogle = window.adsbygoogle || []).push({});
genome and the sequences from Raji of the large deletion found in B95-8
97%
AY825078.1 Human herpesvirus 4 EBNA-1 mRNA, partial cds 96% AJ507799.2 Human herpesvirus 4 complete wild type genome 97%
Bảng 3.3 cho thấy trình tự các chuỗi gen EBNA-1 của EB17 chúng tôi thu nhận đều có tỉ lệ tương đồng cao (từ 97% đến 99%) so với các gen EBNA-1 đã được đăng kí trong Ngân hàng gen của chủng EBV và các chủng EBV khác trên thế giới. Các trình tự gen của chúng tôi thu nhận là gen EBNA-1 của chủng EBV đang nghiên cứu. Điều này cũng khẳng định toàn bộ quá trình kĩ thuật, nhằm thu
nhận gen EBNA-1 của EBV, được chúng tôi thực hiện là hoàn toàn tin cậy và là
cơ sở cho việc xây dựng phương pháp chuẩn đoán nhanh bằng kỹ thuật PCR.
3.2.5. Kết quả xác định phân týp EBV qua phân tích vị trí acid amin số 487 của EBNA-1 thu nhận với chuỗi gen tương ứng của các chủng Việt Nam và Thế giới EBNA-1 thu nhận với chuỗi gen tương ứng của các chủng Việt Nam và Thế giới
Vị trí acid amin số 487 được sử dụng để so sánh đối chiếu EBNA-1 của 4 mẫu EBV thu nhận từ bệnh nhân NPC nhập viện tại Hà Nội (EB17, EB18, EB19, EB20) và một số chủng Việt Nam mới thu nhận gần đây là chủng HP6 (Hải Phòng); chủng EB9 (Hà Nội) đại diện bệnh nhân nhập viện tại miền Bắc và chủng MT7 (Miền Trung) đại diện miền Trung. Một số chủng quốc tế được so sánh, đặc biệt là chủng cổ điển B95-8 (1) và B95-8 (2) được sử dụng để phân biệt prototype; chủng C18, GD1 và HHV (Trung Quốc) đại diện Nam Trung Quốc của châu Á và chủng AG876 (Ghana) đại diện EBV týp 1 và châu Phi. (Hình 3.13).
Kết quả so sánh trình bày ở Hình 3.13 cho thấy, tại vị trí số 487, 4 chủng EBV của Hà Nội chứa 1 loại acid amin, đó là V (valine), từ đó xác định thuộc phân týp: V-val (Variant valine).
Kết quả xác định phân týp của 13 chủng, bao gồm cả 4 chủng nghiên cứu và 9 chủng tham chiếu của Việt Nam và thế giới được liệt kê ở Bảng 3.4. Bảng 3.4 cũng cho thấy chi tiết, cả 4 mẫu nghiên cứu thuộc V-val (EB17, EB18, EB19, EB20); mẫu EB9 (Hà Nội) thuộc V-val; mẫu Miền trung thuộc P-thr (MT7); mẫu Hải Phòng thuộc mẫu P-ala (HP6); 3 mẫu đại diện Nam Trung Quốc: mẫu C18 thuộc P-thr, GD1 thuộc V-val và mẫu HHV thuộc P-ala; 2 chủng cổ điển B95-8 (1) và B95-8 (2) thuộc P-ala; chủng AG876 (Ghana) đại diện EBV týp 1 thuộc V-leu
Hình 3.13. So sánh một phần trình tự acid amin EBNA-1 của 13 chủng nghiên cứu, gồm 4 chủng đang nghiên cứu là EB17; EB18; EB 19; EB20; chủng EB9 của Hà Nội, chủng MT7 của miền trung và chủng HP6 của Hải Phòng. 2 biến thể của chủng cổ điển B95-8; chủng AG876 và 3 chủng của Trung Quốc (HHV, GD1 và C18). Vị trí aa 487 được chỉ dẫn bằng mũi tên; aa ở vị trí số 487 của các chủng EBV so sánh được đóng khung dọc.`
Bảng 3.4. Kết quả xác định phân týp của các chủng EBV qua phân tích acid amin tại vị trí số 487 của EBNA-1
TT Chủng Năm phân lập aa 487 Kết quả xác
định phân týp
1 EB17- EBNA-1 2000 V V-val
aa 487
* 20 * 40 * 60 * 80 * 100 EB17-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPTDDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRVLLARSHVERTTEEGNWVAGVF : 105 EB18-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPTDDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRVLLARSHVERTTEEGNWVAGVF : 105 EB19-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPTDDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRVLLARSHVERTTEEGNWVAGVF : 105 EB20-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRVLLARSHVERTTEEGNWVAGVF : 105 EB9-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPTDDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRVLLARSHVERTTEEGNWVAGVF : 105 MT7-EBNA1 : PFFHPVGEADYFEYHQEGGPDGEPDMPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNQKFENIAEGLRTLLARCHVERTTDEGTWVAGVF : 105 HP6-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPTDDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGHGGSNQKFENIAEGIKALLARSHVERTTEEGNWVAGVF : 105 AG 876-EBN : PFFHPVAEADYFEYHQEGGPDGEPDMPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWYGKHRGEGGSSQKFENIAEGLRLLLARCHVERTTEDGNWVAGVF : 105 B95-8(1)-E : PFFHPVGEADYFEYHQEGGPDGEPDVPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRALLARSHVERTTDEGTWVAGVF : 105 B95-8(2)-E : PFFHPVGEADYFEYHQEGGPDGEPDVPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRALLARSHVERTTDEGTWVAGVF : 105 HHV-EBNA1 : PFFHPVGEADYFEYHQEGGPDGEPDVPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRALLARSHVERTTDEGTWVAGVF : 105 GD1-EBNA1 : PFFHPVGDADYFEYLQEGGPDGEPDVPPGAIEQGPTDDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNPKFENIAEGLRVLLARSHVERTTEEGNWVAGVF : 105 C18-EBNA1 : ----PVGQADYFEYHQEGGPDGEPDMPPGAIEQGPADDPGEGPSTGPRGQGDGGRRKKGGWFGKHRGQGGSNQKFENIADGLRTLLARCHVERTTDEGTWVAGVF : 101
2 EB18- EBNA-1 2000 V V-val
3 EB19- EBNA-1 2000 V V-val
4 EB20- EBNA-1 2000 V V-val
5 EB9- EBNA-1 2000 V V-val
6 MT7- EBNA-1 2009 T P-thr
7 HP6- EBNA-1 2009 A P-ala
8 AG876- EBNA-1 2007 L V-leu 9 B95-8 (1)-EBNA-1 2006 A P-ala 10 B95-8 (2)-EBNA-1 1996 A P-ala
11 HHV- EBNA-1 2006 A P-ala
12 GD1- EBNA-1 2005 V V-val
13 C18 – EBNA-1 2004 T P-thr
Kết quả xác định acid amin vị trí aa số 487 của 4 chủng EBV Hà Nội cho thấy đều thuộc phân týp V-val; các mẫu này được phân lập năm 2000. Theo nghiên cứu của tác giả Lê Thanh Hòa và cs (2010), có khi phân tích 22 mẫu nghiên cứu thu thập các năm 2000 và gần đây (2008-2009), 14 mẫu có phân týp V-val; trong đó tất cả các mẫu năm 2000 đều thuộc phân týp V-val [10]. Như vậy cho đến nay, chưa phát hiện phân týp khác ngoài V-val ở các mẫu NPC thu nhận trước năm 2000. Một điều lưu ý là EBV có đặc tính đa hình biến đổi giữa các chủng phân lập ở các vùng địa lý thuộc Nam và Đông Nam châu Á [68] và nghiên cứu của Mai et al [43], [44] chứng minh các chủng EBV thuộc V-val ở vùng Nam Trung Quốc, có khả năng tăng cường sao chép và có vai trò quan trọng trong tiến triển NPC. Xác định được phân týp (V-val) trong số mẫu EBV chúng tôi nghiên cứu (EB17, EB18, EB19, EB20) và so sánh (EB9) của Hà Nội, giúp chúng tôi khẳng định các mẫu nghiên cứu năm 2000 đều thuộc V-val, còn các mẫu MT7, HP6 năm 2009 xuất hiện thuộc P-ala và P-thr và các chủng so sánh khác của thế giới thuộc V-leu, P-ala và P-thr . Tuy số mẫu chưa nhiều, nhưng nghiên cứu của chúng tôi cho thấy EBV Việt Nam, về mặt dịch tễ học, có bị ảnh hưởng của vệt dịch tễ EBV từ Quảng Đông - Nam Trung Quốc tới Đông Nam châu Á.
Thảo luận chung
- Với kết quả thu được, trên cơ sở trình tự nucleotide và acid amin của HPV-16 và HPV-18 là nguyên nhân gây bệnh ung thư CTC ở các mẫu bệnh phẩm làm đối tượng nghiên cứu thu nhận tại Hà Nội. Mối liên quan giữa loại/týp HPV-16 và HPV-18 đã được xác định và bệnh ung thư CTC ở một số bệnh nhân tại Hà Nội, bước đầu tìm hiểu chính xác nguyên nhân gây bệnh này ở bệnh nhân phụ nữ. Trên cơ sở bước đầu phát hiện HPV-16/HPV-18 bằng phương pháp đa mồi multiplex-PCR chúng ta xây dựng phương pháp chẩn đoán phân tử dựa trên phản ứng PCR nhằm hỗ trợ các phương pháp chẩn đoán tế bào và mô bệnh học. (adsbygoogle = window.adsbygoogle || []).push({});
- Kết quả nghiên cứu đặc điểm sinh học phân tử các chủng EBV của Việt Nam thu nhận dữ liệu từ gen EBNA-1 (gen gây ung thư giai đoạn sớm) trong báo cáo phân tích là một sản phẩm có ý nghĩa quốc gia và quốc tế, đặc biệt về virus EBV của vùng châu Á có xu hướng gây sơ nhiễm, nhiễm tiềm tàng, tái hoạt gây tiến triển NPC. Các trình tự gen của chúng tôi thu nhận là gen EBNA-1 của của chủng EBV đang nghiên cứu. Điều này cũng khẳng định toàn bộ quá trình kĩ thuật, nhằm thu nhận gen EBNA-1 của EBV, được chúng tôi thực hiện là hoàn toàn tin cậy. Xác định được phân týp (V-val) trong số mẫu EBV chúng tôi nghiên cứu (EB17, EB18, EB19, EB20) và so sánh (EB9) của Hà Nội, giúp chúng tôi khẳng định các mẫu nghiên cứu năm 2000 đều thuộc V-val, còn các mẫu MT7, HP6 năm 2009 xuất hiện thuộc P-ala và P-thr và các chủng so sánh khác của thế giới thuộc V-leu, P-ala và P-thr. Tuy số mẫu chưa nhiều, nhưng nghiên cứu của chúng tôi cho thấy EBV Việt Nam, về mặt dịch tễ học, có bị ảnh hưởng của vệt dịch tễ EBV từ Quảng Đông - Nam Trung Quốc tới Đông Nam châu Á là cơ sở cho việc xây dựng phương pháp chuẩn đoán nhanh bằng kỹ thuật PCR.
KẾT LUẬN VÀ KIẾN NGHỊ 1. KẾT LUẬN
Từ những kết quả phân tích và bàn luận ở trên, chúng tôi rút ra một số kết luận chính như sau:
1. Trình tự gen E6 có kích thước 321 bp, gen L1 có kích thước 556 bp, trong
số các mẫu nghiên cứu PK10, PK11, PK12, PS16 đã được thu nhận và phân tích. Kết quả thực hiện phản ứng Multiplex-PCR trên 4 mẫu HPV nghiên cứu là các mẫu song nhiễm cả hai týp HPV-16 và HPV-18.
2. Thực hiện thành công phản ứng Multiplex-PCR phát hiện HPV16/HPV18
trên 11 mẫu nghiên cứu.
3. Trình tự gen EBNA-1 của các mẫu EB17, EB18, EB19, EB20 có kích
thước khoảng 800 bp phân lập tại Hà Nội đã được thu nhận và phân tích.
4. Xác định được phân týp (V-val) trong số mẫu EBV (EB17, EB18, EB19,
EB20) và so sánh (EB9) của Hà Nội, giúp chúng tôi khẳng định các mẫu nghiên cứu năm 2000 đều thuộc V-val, còn các mẫu MT7, HP6 năm 2009 xuất hiện thuộc P-ala và P-thr và các chủng so sánh khác của thế giới thuộc V-leu, P-ala và P-thr.
2. KIẾN NGHỊ
Tiếp tục thu nhận, giải trình tự các mẫu bệnh phẩm bệnh nhân EBV và HPV nhằm thu được dữ liệu sinh học phân tử gen EBNA-1 của chủng EBV và hai gen E6 của chủng HPV-16 và L1 của chủng HPV-18, làm cơ sở cho chẩn đoán, giám định phân loại và xây dựng phương pháp chuẩn đoán nhanh các mẫu bệnh phẩm chủng EBV và HPV.
TÀI LIỆU THAM KHẢO Tiếng Việt
1. Phạm Thị Hoàng Anh, Nguyễn Mạnh Quốc, Nguyễn Bá Đức, Nguyễn Chấn
Hùng (2001), “Tình hình bệnh ung thư ở Việt Nam năm 2000”, Tạp chí thông tin Y dược 2/2001, Bộ Y tế - Viện thông tin Y học Trung ương, tr. 3- 11.
2. Bộ môn Vi Sinh (2001), Vi Sinh Y học, Nhà xuất bản Y học, Hà Nội, tr. 386.
3. Trịnh Quang Diện (1995), Phát hiện các dị sản, loạn sản và ung thư cổ tử
cung bằng phương phát tế bào học, Luận án Tiến sỹ Y học, Trường Đại học
Y Hà Nội.
4. Trịnh Quang Diện (2002), Theo dõi bằng tế bào học và mô bệnh học các tế
bào bào vẩy không điển hình ý nghĩa chưa xác định (ASCUS) gặp trong
phát hiện tế học các tổn thương nội biểu mô và ung thư cổ tử cung, Chuyên
đề ung thư học, Y học thực hành, Nhà xuất bản Y học, số 431, tr. 266- 269.
5. Đào Trung Dũng (2003), Chẩn đoán tế bào học ASCUS trong phát hiện
sớm ung thư cổ tử cung, Luận văn Bác sỹ Y khoa, Trường Đại học Y Hà Nội.
6. Phan Trường Duyệt, Đinh Thế Mỹ (2003), Lâm sàng sản phụ khoa, Nhà
xuất bản Y học, Hà Nội, tr. 444- 448.
7. Nguyễn Văn Đô, Phan Thị Phi Phi (2003), Xác định đột biến mất đoạn 30bp và đột biến điểm của gen EBV-LMP ở bệnh nhân ung thư vòm mũi họng Việt Nam, Hội nghị sinh học phân tử và hoá sinh toàn quốc, tr. 129-134.
8. Phạm Thị Hồng Hà (2000), Giá trị của phiến đồ âm đạo - cổ tử cung, soi cổ (adsbygoogle = window.adsbygoogle || []).push({});
tử cung và mô bệnh học trong việc phát hiện sớm ung thư cổ tử cung, Luận
văn Thạc sỹ Y học, Trường Đại học Y Hà Nội.
9. Lê Thanh Hoà (2006), Y - Sinh học phân tử, Quyển I, Nhà xuất bản y học,
Hà Nội, tr. 64- 85.
10. Lê Thanh Hòa, Vũ Thị Tiến, Hoàng Thị Minh Châu, Lê Thanh Hà, Nguyễn
Đình Phúc (2010), Xác định 3 phân typP-ala, P-thr và V-val của EBNA-1 ở EBV phân lập tại Việt Nam, Tạp chí Nghiên cứu Y học 70(5), tr. 8-12.
11.Vương Tiến Hoà (2004), Một số vấn đề bệnh lý cổ tử cung, Nhà xuất bản Y học, Hà Nội.
12.Trần Phương Mai (1999), Nhận xét 89 trường hợp ung thư cổ tử cung tại
Viện BVBM TSS trong 6 năm (1992- 1997), Tạp chí Y học thực hành số
1/1999, tr. 37- 39.
13.Phan Thị Phi Phi (2005), Báo cáo tổng quan các đặc điểm sinh học chủ yếu
gặp ở bệnh nhân ung thư vòm mũi họng Việt Nam và một số ứng dụng lâm sàng. Hội thảo chuyên đề ung thư vòm mũi họng, Hà Nội (07-09.11.2005).
14.Nguyễn Đình Phúc (2006), “Nghiên cứu chẩn đoán lâm sàng và gen virus
Epstein – Barr trong ung thư vòm mũi họng’’, Luận án TS khoa học Y học
– ĐHYHN.
15.Nguyễn Quốc Trực, Lê Văn Xuân (2000), Chẩn đoán và điều trị các thương
tổn tiền ung thư cổ tử cung, Tạp chí thông tin Y - Dược 8/2000, Viện thông
tin Y học Trung ương, tr. 220- 224.
Tiếng Anh
16.Apt D, Watts RM, Suske G, U Bernard U (1996), “High Sp1/Sp3 ratios in
epithelial cells during epithelial differentiation and cellular transcription correlate with the activation of the HPV-16 promoter”, Virology, 224,
pp. 281– 291.
17.Bosch F, Manos MM, Munoz N, Sherman M, Jansen AM, Peto J,
Schiffman M H, Moreno V, Kurman R, Shah KV, International Biological Study on Cervical Cancer (IBSCC) Study Group (1995), “Prevalence of human papillomavirus in cervical cancer: a worldwide perspective”,J. Natl.
Cancer Inst, 87, pp. 796–802.
18.Bosch FX, Munoz N (2002), “The viral etiology of cervical cancer”, Virus
Res. 89, pp. 183-190.
19.Bosch FX, Lorincz A, Munoz N, Meijer CJ, Shah KV (2002), “The causal
relation between human papillomavirus and cervical cancer”, J. Clin.
20.Chan AT, Teo PM, Huang DP (2004), Pathogenesis and treatment of nasopharyngeal carcinoma. Semin Oncology, 31(6), pp. 794-801. (adsbygoogle = window.adsbygoogle || []).push({});
21.Chan SY, Bernard HU, Ratterree M, Birkebak TA, Faras AJ, Ostrow RS
(1997), “Genomic diversity and evolution of papillomaviruses in rhesus monkeys”, J. Virol. 71, pp. 4938- 4943.
22.Chan SY, Delius H, Halpern AL, Bernard HU (1995), “Analysis of genomic
sequences of 95 papillomavirus types: uniting typing, phylogeny, and taxonomy”, J. Virol. 69, pp. 3074-3083.
23.Chang ET, Adami HO, (2006), The enigmatic epidemiology of
nasopharyngeal carcinoma. Cancer Epidemiol Biomarkers Prev, 15(10),
pp. 1765-77, Review.
24.Chou J, Lin YC, Kim J, You L, Xu Z, He B, Jablons DM (2008),.