Genomika Populasi Naga Komodo: Ulasan Makalah Iannucci dkk. (2021)

📖 14 menit baca~3086 kata

Daftar Isi

• Ringkasan Makalah
• Latar Belakang: Mengapa Genomika Populasi Penting
• Metode: Pengambilan Sampel dan Sekuensing Genom Utuh
• Temuan Utama: Heterozigositas Lintas Subpopulasi
• Koefisien Perkawinan Sedarah dan Substruktur Populasi
• Estimasi Ukuran Populasi Efektif (Ne)
• Perbandingan dengan Hasil Mikrosatelit Ciofi 1999
• Implikasi bagi Manajemen Konservasi
• Keterbatasan dan Arah Penelitian Mendatang
• Fakta Cepat
• Mitos vs Fakta
• Poin Penting
• Pertanyaan yang Sering Diajukan
• Sumber & Bacaan Lanjutan

Ringkasan Makalah

Pada tahun 2021, sebuah tim penelitian internasional yang dipimpin Alessio Iannucci dari Universitas Florence menerbitkan studi tingkat populasi paling komprehensif secara genomik tentang naga Komodo (Varanus komodoensis) hingga saat ini. Makalah tersebut — Iannucci, A., Benazzo, A., Natali, C., Arida, E., Zein, M., Jessop, T.S., Bertorelle, G., & Ciofi, C. (2021). "Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing." Molecular Ecology 30(23): 6309–6324. DOI: 10.1111/mec.16121 — menggunakan sekuensing genom utuh dari 24 naga Komodo liar yang diambil dari seluruh lima pulau utama dalam jangkauan spesies.

Temuan utamanya cukup mengkhawatirkan. Di seluruh populasi yang diambil sampelnya, keragaman genomik V. komodoensis ternyata lebih rendah dari nilai yang dilaporkan untuk spesies reptil lain yang dianggap rentan atau terancam kritis. Runs of homozygosity yang meningkat — tanda genomik perkawinan sedarah baru-baru ini — terdeteksi di berbagai populasi pulau, terutama di Flores Utara dan Nusa Kode. Pemodelan demografi mengungkap bahwa ukuran populasi efektif (Ne) spesies ini telah anjlok dalam 150.000 tahun terakhir, hanya menyisakan beberapa ratus individu pada masa geologi yang relatif baru. Dan di mana studi mikrosatelit sebelumnya mengidentifikasi dua kluster genetik, analisis seluruh genom menemukan tiga garis keturunan genomik berbeda.

Secara keseluruhan, temuan ini menggambarkan spesies di bawah tekanan genetik berlapis: bukan sekadar kecil dalam hitungan sensus, tetapi miskin secara genetik dengan cara yang melemahkan potensi adaptif, kebugaran reproduksi, dan ketahanan jangka panjang. Makalah ini diterbitkan pada tahun yang sama IUCN meningkatkan status Daftar Merah komodo dari Rentan menjadi Terancam Punah.

Latar Belakang: Mengapa Genomika Populasi Penting

Selama sebagian besar abad kedua puluh, manajemen konservasi spesies yang terancam punah bergantung pada dua jenis informasi: jumlah sensus dan observasi lapangan. Berapa banyak hewan yang tersisa? Di mana mereka tinggal? Apakah mereka bereproduksi? Data ini tetap penting, tetapi tidak lengkap. Populasi tiga ribu hewan mungkin membawa bahan genetik yang jauh lebih sedikit dari yang dikira jumlah sensusnya, jika ketiga ribu individu itu semuanya berasal dari bottleneck baru-baru ini, jika isolasi geografis telah mencegah persilangan, atau jika proporsi tinggi perkawinan terjadi antar kerabat dekat.

Genomika populasi — analisis variasi skala-genom di berbagai individu dari berbagai lokasi geografis — mengatasi kesenjangan ini secara langsung. Dengan memeriksa polimorfisme nukleotida tunggal (SNP) yang tersebar di seluruh genom, peneliti dapat memperkirakan: seberapa besar keragaman genetik yang dipertahankan berbagai populasi; apakah populasi berbeda telah kawin bebas atau telah terisolasi secara reproduksi; bagaimana ukuran populasi efektif berubah selama ribuan hingga jutaan tahun; dan apakah individu tertentu membawa tanda perkawinan sedarah baru-baru ini berupa bentangan DNA homozigot panjang yang dikenal sebagai runs of homozygosity (ROH).

Metode: Pengambilan Sampel dan Sekuensing Genom Utuh

Studi ini mengambil sampel 24 naga Komodo dari enam populasi geografis di lima pulau Indonesia: Komodo (n=4), Rinca (n=4), Nusa Kode (n=4), Gili Motang (n=4), Flores Barat (n=4), dan Flores Utara (n=4). Sekuensing dilakukan pada platform Illumina dengan strategi berlapis: enam sampel disekuensing dengan kedalaman sekitar 25× dan 18 sisanya dengan sekitar 10×. Pembacaan diselaraskan ke genom referensi Varanus komodoensis yang diterbitkan Lind dkk. (2019).

Setelah penyaringan kualitas ketat, tim mempertahankan 608.471 SNP yang tersebar di autosom. Dataset ini menjadi dasar semua analisis berikutnya: analisis komponen utama (PCA) dan pemodelan ADMIXTURE untuk struktur populasi; perhitungan heterozigositas menggunakan estimator Watterson (theta-W); deteksi ROH untuk koefisien perkawinan sedarah tingkat individu (F_ROH); pairwise F_ST untuk diferensiasi antar populasi; dan rekonstruksi sejarah demografi menggunakan pemodelan koalesen sekuensial-Markovian berganda (MSMC) pada sampel dengan kedalaman lebih tinggi.

Temuan Utama: Heterozigositas Lintas Subpopulasi

Ukuran keragaman genetik yang paling langsung adalah heterozigositas nukleotida — proporsi lokus genomik di mana individu membawa dua alel berbeda. Studi Iannucci dkk. (2021) memperkirakan heterozigositas per individu menggunakan estimator theta-W Watterson. Nilai berkisar dari sekitar 7,62 × 10⁻⁵ hingga 1,31 × 10⁻⁴ di seluruh individu, dengan perbedaan nyata antar populasi.

Gili Motang — pulau terkecil yang diambil sampelnya — memiliki rata-rata theta-W terendah sekitar 7,37 × 10⁻⁵, konsisten dengan efek populasi pendiri kecil dan hanyutan genetik dalam isolasi. Flores Barat menunjukkan rata-rata tertinggi sekitar 9,81 × 10⁻⁵. Nilai-nilai ini bermakna dalam konteks komparatif: keragaman genomik komodo berada di bawah nilai yang dilaporkan untuk spesies yang diklasifikasikan sebagai rentan atau bahkan terancam kritis.

Nilai pairwise F_ST antara populasi lebih lanjut mengkonfirmasi kedalaman diferensiasi genetik. Flores Utara menunjukkan divergensi tertinggi dari semua populasi lain, konsisten dengan periode isolasi geografis dan genetiknya yang lebih lama. Pola-pola ini menunjukkan bahwa fragmentasi genetik telah maju pada tingkat yang biasanya memicu kekhawatiran manajemen dalam populasi vertebrata mana pun.

Koefisien Perkawinan Sedarah dan Substruktur Populasi

Sementara heterozigositas mengukur keragaman di tingkat populasi, runs of homozygosity (ROH) memberikan ukuran tingkat individu dari perkawinan sedarah. ROH adalah bentangan panjang genom di mana kedua salinan kromosom identik — tanda telltale bahwa kedua orang tua individu tersebut berbagi nenek moyang bersama baru-baru ini. Proporsi autosom individu yang ditutupi segmen ROH lebih panjang dari 1 Mb disebut F_ROH, dan ini melacak perkawinan sedarah yang terealisasi dengan presisi lebih tinggi dari estimasi pedigree.

Studi Iannucci dkk. (2021) menemukan F_ROH yang meningkat di berbagai populasi pulau, dengan nilai tertinggi terkonsentrasi di subpopulasi yang paling terpencil dan paling sedikit terhubung. Rata-rata F_ROH menurut populasi (untuk ROH ≥ 1 Mb) tertinggi untuk Flores Utara sekitar 12%, diikuti Nusa Kode sekitar 11%, lalu Komodo sekitar 10%. Flores Barat, Rinca, dan Gili Motang menunjukkan nilai rata-rata F_ROH di bawah 10%.

Yang sangat mencolok adalah individu Flores Utara yang ditandai NF1, yang membawa F_ROH sekitar 30%. Nilai ekstrem seperti ini pada individu liar menunjukkan perkawinan yang sangat baru antara kerabat dekat, atau sejarah panjang perkawinan dalam populasi yang sangat kecil di mana semua individu saling berkerabat melalui berbagai jalur silsilah.

Analisis ADMIXTURE mendukung model tiga kluster sebagai deskripsi paling parsimoni tentang struktur genomik. Kluster 1 sebagian besar berkorespondensi dengan individu Pulau Komodo; Kluster 2 dengan Flores Utara; dan Kluster 3 dengan leluhur yang tersebar luas di Rinca, Nusa Kode, Gili Motang, dan Flores Barat. Analisis komponen utama (PCA) menceritakan kisah yang sama dalam dua dimensi, mencerminkan geografi fisik Kepulauan Sunda Kecil.

Estimasi Ukuran Populasi Efektif (Ne)

Rekonstruksi demografi berbasis MSMC menelusuri dinamika Ne komodo selama sekitar satu juta tahun terakhir, memanfaatkan enam genom dengan kedalaman tinggi. Hasilnya mengungkap sejarah demografi yang luas dan konsisten dengan spesies yang dahulu jauh lebih banyak dan tersebar luas. Sekitar satu juta hingga 500.000 tahun lalu, Ne tampaknya sangat besar — puluhan ribu atau lebih tinggi — konsisten dengan spesies yang menempati jangkauan lebih luas selama periode interglasial yang lebih hangat.

Fitur paling mencolok dari lintasan demografi adalah penurunan cepat dan parah yang dimulai sekitar 150.000 tahun lalu dan berlanjut melalui Holosen hingga sekitar 3.000 tahun sebelum sekarang. Pada akhir lintasan ini, estimasi Ne di semua populasi telah menyatu pada nilai dalam ratusan rendah — jauh di bawah apa yang dianggap ukuran efektif genetik minimum yang layak oleh ahli genetika populasi.

Analisis RCCR mengindikasikan bahwa aliran gen antara sebagian besar pasangan populasi secara efektif berhenti selama periode yang sama: populasi Flores Utara terisolasi sekitar 15.000 tahun lalu; Pulau Komodo terisolasi sekitar 5.000 tahun lalu; dan kelompok pulau tengah (Rinca, Nusa Kode, Gili Motang) mempertahankan beberapa konektivitas hingga sekitar 1.000 tahun lalu. Waktu isolasi ini erat mengikuti sejarah naiknya permukaan laut di kawasan Sunda Kecil.

Perbandingan dengan Hasil Mikrosatelit Ciofi 1999

Studi 2021 tidak muncul dalam kevakuman. Baseline komparatif dasar untuk genetika populasi komodo adalah sepasang makalah yang diterbitkan Ciofi dan rekan-rekan pada pergantian milenium. Ciofi, Beaumont, Swingland, & Bruford (1999) menganalisis variasi alel pada lokus mikrosatelit dan melaporkan bahwa Pulau Komodo menunjukkan divergensi genetik terbesar dan harus dikelola secara terpisah.

Secara garis besar, studi seluruh-genom 2021 mengkonfirmasi temuan 1999: populasi pulau terdiferensiasi secara genetik, Ne rendah, dan Pulau Komodo menempati posisi genomik berbeda. Namun, resolusi yang lebih tinggi dari data SNP seluruh-genom menghasilkan revisi penting. Yang paling signifikan, di mana analisis mikrosatelit mengenali dua kluster genetik utama, analisis ADMIXTURE dari 608.471 SNP mengungkap tiga kluster — dengan Flores Utara muncul sebagai unit genomik berbeda yang terpisah dari Pulau Komodo dan kelompok pulau tengah. Perbedaan ini tidak terlihat dari data mikrosatelit, kemungkinan karena 10 lokus yang digunakan dalam studi 1999 kurang memiliki kekuatan statistik untuk menyelesaikan diferensiasi tambahan tersebut.

Implikasi bagi Manajemen Konservasi

Para penulis studi 2021 eksplisit dalam menerjemahkan temuan genomik mereka ke dalam rekomendasi konservasi. Pada tingkat yang paling mendasar, identifikasi tiga kluster genomik berbeda menuntut pengakuan Pulau Komodo, Flores Utara, dan kelompok pulau tengah (Rinca, Nusa Kode, Gili Motang, Flores Barat) sebagai unit manajemen terpisah. Secara praktis, ini berarti program translokasi atau pembiakan penangkaran tidak boleh mencampurkan individu dari garis keturunan yang berbeda secara genomik ini tanpa pertimbangan cermat tentang konsekuensinya — termasuk risiko outbreeding depression.

Populasi Flores Utara memerlukan perhatian khusus dari berbagai sudut pandang. Ini adalah yang paling terisolasi secara geografis, paling berbeda secara genetik, membawa nilai F_ROH tertinggi, dan — yang sangat krusial — hampir seluruhnya berada di luar batas kawasan lindung Taman Nasional Komodo. Tidak seperti komodo di Pulau Komodo dan Rinca yang mendapat perlindungan hukum dan patroli ranger, komodo Flores Utara menghuni lanskap pesisir yang didominasi manusia, di mana kehilangan habitat, berkurangnya mangsa, dan persekusi langsung mengancam kelangsungannya. Bukti genomik 2021 memberikan argumen ilmiah yang kuat untuk memperluas perlindungan formal ke populasi ini.

Populasi Gili Motang juga perlu mendapat perhatian. Pulau kecil ini sudah diidentifikasi oleh studi demografi sebelumnya (terutama Purwandana dkk., 2014) sebagai berpotensi menurun. Dikombinasikan dengan bukti genomik heterozigositas rendah dan isolasi yang tinggi, Gili Motang mewakili populasi dalam bahaya ganda: rapuh secara demografis dan miskin secara genetik. Para penulis merekomendasikan Gili Motang sebagai kandidat untuk augmentasi atau penyelamatan genetik.

Keterbatasan dan Arah Penelitian Mendatang

Studi Iannucci dkk. (2021) memiliki keterbatasan yang diakui para penulisnya. Yang paling signifikan adalah ukuran sampel: dua puluh empat individu dari enam populasi — empat per populasi — adalah minimum yang diperlukan untuk menghasilkan statistik tingkat populasi yang bermakna, tetapi tidak cukup untuk menangkap seluruh rentang variasi genomik tingkat individu dalam populasi mana pun.

Kedua, kedalaman sekuensing tidak homogen. Hanya enam sampel disekuensing pada kedalaman 25×; sisanya sekitar 10×. Sekuensing kedalaman rendah memperkenalkan kesalahan stokastik dalam pemanggilan genotipe yang dapat menggelembungkan atau mengempiskan estimasi ROH. Para penulis menerapkan ambang penyaringan konservatif untuk mengurangi hal ini, tetapi beberapa ketidakpastian tetap ada.

Penelitian mendatang harus memprioritaskan setidaknya tiga arah. Pertama, pengambilan sampel genomik populasi Flores Utara harus diperluas untuk mengkarakterisasi rentang penuh tingkat perkawinan sedarah. Kedua, analisis genomik fungsional akan memungkinkan peneliti untuk bergerak dari gambaran perkawinan sedarah yang murni deskriptif ke penilaian konsekuensi kebugarannya. Ketiga, integrasi data genomik dengan data lapangan demografis akan memungkinkan estimasi Ne ditempatkan dalam konteks tren ukuran populasi aktual.

Fakta Cepat

Parameter	Nilai / Rincian
Diterbitkan	2021 (online Agustus 2021; cetak Desember 2021)
Jurnal	Molecular Ecology, Vol. 30, Edisi 23
Halaman	6309–6324
DOI	10.1111/mec.16121
Penulis utama	Alessio Iannucci (Universitas Florence, Italia)
Metode	Sekuensing genom utuh (Illumina); ADMIXTURE; inferensi demografi MSMC2; analisis ROH
Ukuran sampel	24 individu liar dari 6 populasi di 5 pulau
SNP yang dipertahankan	608.471 (setelah penyaringan kualitas)
Kluster genomik (K)	3 (Pulau Komodo; Flores Utara; Rinca + Nusa Kode + Gili Motang + Flores Barat)
F_ROH tertinggi	Flores Utara (~12%); individu maksimum: NF1 (~30%)
Tanggal isolasi Flores Utara	~15.000 tahun lalu (analisis RCCR)
Tanggal isolasi Pulau Komodo	~5.000 tahun lalu (analisis RCCR)

Mitos vs Fakta

Asumsi Umum	Apa yang Ditunjukkan Genomika
Populasi ~3.000 hewan cukup besar untuk mempertahankan kesehatan genetik.	Jumlah sensus tidak relevan jika Ne sangat kecil. Estimasi Ne Holosen akhir hanya beberapa ratus menunjukkan bottleneck genetik dalam meskipun total populasi lebih besar.
Taman Nasional Komodo melindungi populasi yang paling penting secara genetik.	Flores Utara — yang sebagian besar berada di luar taman — adalah populasi yang paling berbeda secara genomik dan membawa koefisien perkawinan sedarah tertinggi, menjadikannya prioritas konservasi tertinggi.
Pulau Komodo dan Rinca pada dasarnya adalah stok genetik yang sama.	Analisis ADMIXTURE dan PCA menyelesaikan Pulau Komodo sebagai kluster genomik berbeda, terpisah dari kelompok yang berpusat di Rinca. Keduanya harus dikelola sebagai unit konservasi yang berbeda.
Studi mikrosatelit 1999 mengidentifikasi semua unit konservasi penting.	Data seluruh-genom mengungkap kluster ketiga (Flores Utara) yang tidak dapat diselesaikan oleh analisis mikrosatelit. Kerangka manajemen yang hanya didasarkan pada data 1999 tidak lengkap.
Translokasi individu antar pulau adalah alat konservasi yang mudah untuk komodo.	Mengingat diferensiasi genomik yang dalam antar kluster, translokasi tanpa pengawasan berisiko menimbulkan outbreeding depression. Program augmentasi apa pun memerlukan skrining genomik cermat.

Poin Penting

• Tiga kluster genomik, bukan dua. Spesies ini memiliki setidaknya tiga garis keturunan penting secara evolusioner yang memerlukan pengelolaan terpisah: Pulau Komodo, Flores Utara, dan kluster pulau tengah. Perencanaan konservasi harus memperlakukan ini sebagai entitas berbeda.
• Flores Utara adalah kasus darurat. Populasi yang paling terisolasi secara genetik dan paling banyak perkawinan sedarah sebagian besar hidup di luar batas perlindungan. Memperluas perlindungan formal ke habitat Flores Utara adalah prioritas konservasi yang dapat dikuantifikasi.
• Penurunan Ne adalah tren jangka panjang, bukan krisis baru-baru ini. Bottleneck dimulai sekitar 150.000 tahun lalu, artinya spesies ini telah miskin secara genetik sebelum kedatangan manusia modern di kawasan ini. Ancaman saat ini mempercepat kerentanan yang sudah ada sebelumnya.
• Gili Motang memerlukan pemantauan aktif. Heterozigositas rendah, ukuran pulau kecil, fragilitas demografis yang didokumentasikan dalam studi lapangan sebelumnya, dan isolasi tinggi menjadikan Gili Motang kandidat untuk augmentasi genetik jika tren populasinya menurun.
• Reproduktibilitas adalah kekuatan. Semua data mentah dan skrip diarsipkan secara publik (PRJNA738464; GitHub). Analisis ulang independen dan peneliti masa depan dapat membangun langsung di atas fondasi ini.

Pertanyaan yang Sering Diajukan

Mengapa Iannucci dkk. (2021) memilih sekuensing genom utuh alih-alih mikrosatelit?

Sekuensing genom utuh menghasilkan ratusan ribu SNP yang tersebar di seluruh genom, memungkinkan estimasi jauh lebih tepat tentang struktur populasi, perkawinan sedarah, dan sejarah demografi dibanding 10–15 lokus mikrosatelit yang tersedia pada akhir 1990-an. Data SNP skala-genom juga memungkinkan analisis runs-of-homozygosity (ROH) yang mengidentifikasi perkawinan sedarah individu dengan sensitivitas jauh lebih tinggi. Studi 2021 menghasilkan 608.471 SNP yang difilter dari 24 individu — beberapa orde besarnya lebih banyak data dibanding studi genetika komodo sebelumnya manapun.

Populasi komodo mana yang menunjukkan perkawinan sedarah tertinggi dalam studi 2021?

Individu dari Flores Utara menunjukkan rata-rata F_ROH tertinggi sekitar 12%, dengan satu individu (NF1) membawa koefisien perkawinan sedarah genomik sebesar 30% — artinya hampir sepertiga genomnya autozigot. Nusa Kode berada di belakang dengan rata-rata sekitar 11%. Nilai-nilai ini tinggi dibanding populasi reptil yang dipantau dan menunjukkan perkawinan sedarah substansial baru-baru ini, kemungkinan akibat ukuran populasi kecil dan aliran gen terbatas ke subpopulasi perifer ini.

Apa yang dimaksud 'ukuran populasi efektif' dan mengapa estimasi Ne dalam studi ini mengkhawatirkan?

Ukuran populasi efektif (Ne) adalah ukuran populasi ideal yang kawin secara acak dan kehilangan keragaman genetik pada laju yang sama dengan populasi nyata yang dikaji. Hampir selalu lebih kecil dari jumlah sensus. Studi 2021 memperkirakan bahwa Ne leluhur telah turun drastis selama 150.000 tahun terakhir, mencapai hanya beberapa ratus individu sekitar 3.000 tahun lalu — lintasan yang mengindikasikan spesies yang berada di ambang viabilitas genetik jangka panjang.

Bagaimana studi 2021 mendefinisikan ulang unit konservasi komodo?

Sebelum 2021, manajemen konservasi sebagian besar memperlakukan populasi Taman Nasional Komodo sebagai satu unit genetik. Analisis ADMIXTURE dan koalesen studi 2021 mengungkap bahwa Pulau Komodo dan pesisir Flores Utara merupakan kluster genomik berbeda dengan sejarah evolusi terpisah dan waktu penghentian aliran gen independen — Flores Utara terisolasi sekitar 15.000 tahun lalu, Pulau Komodo sekitar 5.000 tahun lalu. Ini mendukung kuat pengelolaan masing-masing sebagai unit konservasi terpisah, dengan urgensi khusus untuk populasi Flores Utara yang tidak terlindungi.

Apakah temuan genomik 2021 memengaruhi status Daftar Merah IUCN komodo?

IUCN meningkatkan status Varanus komodoensis dari Rentan menjadi Terancam Punah pada 2021, mengutip jangkauan terbatas, subpopulasi kecil yang terisolasi, dan kehilangan habitat yang diproyeksikan akibat perubahan iklim. Bukti genomik yang diterbitkan Iannucci dkk. pada tahun yang sama menguatkan peningkatan status tersebut dengan menyediakan bukti kuantitatif independen tentang keragaman genomik rendah, perkawinan sedarah meningkat, dan penurunan jangka panjang ukuran populasi efektif. Meskipun penilaian IUCN terutama didorong oleh data sensus dan pemodelan iklim, studi genomik memperkuat konsensus ilmiah bahwa spesies ini menghadapi risiko kepunahan yang lebih parah dari yang sebelumnya diakui.

Sumber & Bacaan Lanjutan

1. Iannucci, A., Benazzo, A., Natali, C., Arida, E., Zein, M., Jessop, T.S., Bertorelle, G., & Ciofi, C. (2021). "Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing." Molecular Ecology, 30(23), 6309–6324. https://doi.org/10.1111/mec.16121
2. Ciofi, C., Beaumont, M.A., Swingland, I.R., & Bruford, M.W. (1999). "Genetic divergence and units for conservation in the Komodo dragon Varanus komodoensis." Proceedings of the Royal Society Series B, 266(1435), 2269–2274. https://doi.org/10.1098/rspb.1999.0918
3. Lind, A.L. et al. (2019). "Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards." Nature Ecology & Evolution, 3, 1241–1252. https://doi.org/10.1038/s41559-019-0945-8
4. Purwandana, D., Ariefiandy, A., Imansyah, M.J., Rudiharto, H., Seno, A., Ciofi, C., Fordham, D.A., & Jessop, T.S. (2014). "Demographic status of Komodo dragon populations in Komodo National Park." Biological Conservation, 171, 29–35. https://doi.org/10.1016/j.biocon.2014.01.007
5. Auffenberg, W. (1981). The Behavioral Ecology of the Komodo Monitor. University Presses of Florida. Monografi lapangan dasar tentang ekologi dan sejarah alam komodo.
6. IUCN SSC Monitor Lizard Specialist Group. (2021). Varanus komodoensis. The IUCN Red List of Threatened Species 2021. https://doi.org/10.2305/IUCN.UK.1996.RLTS.T22884A9396736.en