Analisi genomica della domesticazione e dell'evoluzione post-conquista spagnola del lama e dell'alpacaGenomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca

Titolo originale

«Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca»

Titolo in italiano (traduzione redazionale): «Analisi genomica della domesticazione e dell'evoluzione post-conquista spagnola del lama e dell'alpaca».

Original title

«Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca»

Italian title (editorial translation): «Analisi genomica della domesticazione e dell'evoluzione post-conquista spagnola del lama e dell'alpaca».

Autori

Fan Ruiwen · Gu Zhongru · Guang Xuanmin · Marín Juan Carlos · Varas Valeria · González Benito A. · Wheeler Jane C. · Hu Yafei · Li Erli · Sun Xiaohui · Yang Xukui · Zhang Chi · Gao Wenjun · He Junping · Munch Kasper · Corbett-Detig Russel · Barbato Mario · Pan Shengkai · Zhan Xiangjiang · Bruford Michael W. · Dong Changsheng.

Elenco completo, ventuno autori, nell'ordine e nella forma del record Crossref (interrogato il 31/08/2026). Il record Crossref non riporta ORCID per nessuno dei ventuno autori, e questa scheda non ne inventa.

Autori corrispondenti (tre, come indicato nel frontespizio del deposito PMC): Zhan Xiangjiang, Bruford Michael W., Dong Changsheng. Cinque autori sono dichiarati come contributori in parti uguali: Fan, Gu, Zhan, Bruford, Dong.

Affiliazioni, come dichiarate nell'articolo (quattordici sedi, sei paesi): College of Animal Science and Veterinary Medicine, Shanxi Agricultural University, Taigu, Cina · CAS Key Lab of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Pechino · Cardiff University – Institute of Zoology Joint Laboratory for Biocomplexity Research, CAS, Pechino · University of Chinese Academy of Sciences, Pechino · BGI Genomics, Shenzhen · Universidad del Bío-Bío, Chillán, Cile · Universidad Austral de Chile, Valdivia · Universidad de Chile, Santiago · CONOPA — Instituto de Investigación y Desarrollo de Camélidos Sudamericanos, Pachacamac, Lima, Perù · Bioinformatics Research Centre, Aarhus University, Danimarca · Genomics Institute, UC Santa Cruz, Stati Uniti · Department of Animal Science, Food and Technology (DIANA), Università Cattolica del Sacro Cuore, Piacenza, Italia · Center for Excellence in Animal Evolution and Genetics, CAS, Kunming · School of Biosciences and Sustainable Places Institute, Cardiff University, Galles.

Due nomi vanno segnalati a un lettore italiano. Il primo è Jane C. Wheeler, del CONOPA di Lima: è l'archeozoologa che per prima ha sostenuto, su base morfologica e poi molecolare, che l'alpaca discende dalla vigogna e non dal guanaco — cioè la tesi che questo lavoro conferma con i genomi. Il secondo è Mario Barbato, dell'Università Cattolica di Piacenza: c'è un italiano fra i ventuno autori del genoma del lama, e vale la pena saperlo quando si dice che in Italia sui camelidi non lavora nessuno.

Authors

Fan Ruiwen · Gu Zhongru · Guang Xuanmin · Marín Juan Carlos · Varas Valeria · González Benito A. · Wheeler Jane C. · Hu Yafei · Li Erli · Sun Xiaohui · Yang Xukui · Zhang Chi · Gao Wenjun · He Junping · Munch Kasper · Corbett-Detig Russel · Barbato Mario · Pan Shengkai · Zhan Xiangjiang · Bruford Michael W. · Dong Changsheng.

Complete list, twenty-one authors, in the order and form of the Crossref record (queried on 31/08/2026). The Crossref record reports no ORCID for any of the twenty-one authors, and this factsheet does not invent one.

Corresponding authors (three, as indicated on the frontispiece of the PMC deposit): Zhan Xiangjiang, Bruford Michael W., Dong Changsheng. Five authors are declared as contributing in equal parts: Fan, Gu, Zhan, Bruford, Dong.

Affiliations, as declared in the article (fourteen sites, six countries): College of Animal Science and Veterinary Medicine, Shanxi Agricultural University, Taigu, China · CAS Key Lab of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing · Cardiff University – Institute of Zoology Joint Laboratory for Biocomplexity Research, CAS, Beijing · University of Chinese Academy of Sciences, Beijing · BGI Genomics, Shenzhen · Universidad del Bío-Bío, Chillán, Chile · Universidad Austral de Chile, Valdivia · Universidad de Chile, Santiago · CONOPA — Instituto de Investigación y Desarrollo de Camélidos Sudamericanos, Pachacamac, Lima, Peru · Bioinformatics Research Centre, Aarhus University, Denmark · Genomics Institute, UC Santa Cruz, United States · Department of Animal Science, Food and Technology (DIANA), Università Cattolica del Sacro Cuore, Piacenza, Italy · Center for Excellence in Animal Evolution and Genetics, CAS, Kunming · School of Biosciences and Sustainable Places Institute, Cardiff University, Wales.

Two names deserve a note for an Italian reader. The first is Jane C. Wheeler, of CONOPA in Lima: she is the archaeozoologist who first argued, on morphological and later molecular grounds, that the alpaca descends from the vicuña and not the guanaco — the very thesis this work confirms with genomes. The second is Mario Barbato, of the Università Cattolica di Piacenza: there is an Italian among the twenty-one authors of the llama genome, and it is worth knowing when it is said that in Italy no one works on camelids.

Rivista, anno, volume, fascicolo, pagine, identificatori

Genome Biology, 2020, volume 21, fascicolo 1, numero d'articolo 159. 159 è il numero d'articolo, non una pagina. Pubblicazione elettronica: 2 luglio 2020; collection date 2020. Cronologia editoriale dichiarata: ricevuto l'8 ottobre 2019, accettato il 21 giugno 2020 — otto mesi e mezzo di revisione, che sono l'esatto contrario del ritmo della scheda precedente e altrettanto informativi.

• DOI 10.1186/s13059-020-02080-6
• PMID 32616020 · PMCID PMC7331169
• Testo completo libero: sì, depositato in PubMed Central ed Europe PMC, licenza CC BY 4.0 (letta su Crossref e sul deposito PMC il 31/08/2026).
• Tipo di pubblicazione secondo Europe PMC: studio comparativo, ricerca con sostegno pubblico non statunitense, articolo di ricerca, articolo su rivista.
• La storia della revisione paritaria è pubblicata come Additional file 2, ed è una scelta editoriale della rivista che questa collana considera un elemento di qualità: chi vuole sapere che cosa i revisori hanno contestato può leggerlo.
• Accesso aperto dell'articolo: sì. Unpaywall, interrogato il 03/09/2026: accesso aperto pieno presso l'editore, versione pubblicata presso l'editore, licenza CC BY.

Stato della rivista (verificato il 31/08/2026): Genome Biology, ISSN elettronico 1474-760X (ISSN a stampa 1474-7596), NLM ID 100960660, editore BioMed Central / Springer Nature. PubMed: sì. DOAJ: sì. OpenAlex: 9.154 lavori indicizzati, rivista interamente open access. È una delle testate di riferimento della genomica: che un lavoro sui camelidi sudamericani ci sia entrato dopo otto mesi e mezzo di revisione è di per sé un'informazione.

Licenza dell'articolo originale: CC BY 4.0. Questa scheda può riprodurre: dati, tabelle e figure con attribuzione. In concreto riproduce solo dati numerici, riorganizzati.

Journal, year, volume, issue, pages, identifiers

Genome Biology, 2020, volume 21, issue 1, article number 159. 159 is the article number, not a page. Electronic publication: 2 July 2020; collection date 2020. Declared editorial timeline: received 8 October 2019, accepted 21 June 2020 — eight and a half months of review, which is the exact opposite of the pace in the previous factsheet, and equally informative.

• DOI 10.1186/s13059-020-02080-6
• PMID 32616020 · PMCID PMC7331169
• Free full text: yes, deposited in PubMed Central and Europe PMC, license CC BY 4.0 (checked on Crossref and on the PMC deposit on 31/08/2026).
• Publication type per Europe PMC: comparative study, non-US-funded research, research article, journal article.
• The peer review history is published as Additional file 2, and it is an editorial choice of the journal that this collection regards as a quality element: anyone who wants to know what the reviewers challenged can read it.
• Open access of the article: yes. Unpaywall, queried on 03/09/2026: full open access at the publisher, published version at the publisher, license CC BY.

Journal status (verified 31/08/2026): Genome Biology, electronic ISSN 1474-760X (print ISSN 1474-7596), NLM ID 100960660, publisher BioMed Central / Springer Nature. PubMed: yes. DOAJ: yes. OpenAlex: 9,154 indexed works, fully open-access journal. It is one of the reference journals of genomics: that a work on South American camelids got in after eight and a half months of review is in itself informative.

License of the original article: CC BY 4.0. This factsheet may reproduce: data, tables and figures with attribution. Concretely it reproduces numerical data only, reorganized.

Specie e soggetti

Quattro specie, due generi, e la distinzione fra loro è letteralmente l'oggetto dello studio:

⚠️ Punto che quasi tutte le divulgazioni sbagliano: in questo lavoro il genoma dell'alpaca non è nuovo. I tre genomi assemblati de novo qui sono lama, guanaco e vigogna; l'alpaca entra con il genoma di riferimento pubblicato da altri. Chi cita questo articolo come «il genoma dell'alpaca» lo cita male.

Numerosità — e qui va detta una cosa scomoda, perché l'articolo la riporta in tre modi diversi.

• L'abstract dice: «resequence seven individuals of all four species» — sette individui per ciascuna delle quattro specie.
• La sezione Metodi dice: campioni di sangue e tessuto ottenuti da otto guanachi selvatici, otto vigogne, otto lama domestici e sette alpaca.
• La legenda della Figura 2 indica N = 8 individui sia per l'alpaca sia per il lama negli assi delle proporzioni di ascendenza; la sezione sulla storia demografica dice di aver usato «6 (alpaca e lama) o 8 (vigogna e guanaco)» genomi riesequenziati per l'analisi MSMC.
• L'albero filogenetico è costruito su 32 camelidi sudamericani in tutto (dichiarato nella legenda della Figura 1c per l'analisi ADMIXTURE).

Questa scheda non sceglie per gli autori quale numero sia quello giusto. Riporta i tre e segnala che l'ordine di grandezza è sette-otto individui per specie, che è la cifra da cui dipende ogni intervallo di confidenza dell'articolo. La ripartizione esatta per analisi: [DA VERIFICARE SUL TESTO PIENO, tabelle supplementari S8 e S9].

Provenienza dei campioni, come dichiarata: guanachi e vigogne campionati su tutto l'areale di distribuzione di ciascuna specie; per gli assemblaggi de novo, un guanaco maschio adulto di Putre (Cile) e una vigogna femmina adulta del Parco Nazionale Lauca (Cile), scelti come rappresentanti delle sottospecie settentrionali Lama guanicoe cacsilensis e Vicugna vicugna mensalis; tessuti di lama da Putre e di alpaca da Temuco (Cile). Per i dati di espressione, campioni di cute di 17 alpaca dalla Cina.

Nessun animale europeo. Nessuna alpaca italiana, nessun lama europeo: le proporzioni di introgressione misurate qui sono di animali sudamericani e cinesi, e la trasferibilità alle mandrie europee è un'assunzione ragionevole ma non misurata.

Species and subjects

Four species, two genera, and the distinction between them is literally the object of the study:

⚠️ A point almost every popularization gets wrong: in this work the alpaca genome is not new. The three genomes assembled de novo here are llama, guanaco and vicuña; the alpaca enters with the reference genome published by others. Anyone who cites this article as "the alpaca genome" is citing it wrongly.

Sample size — and here something uncomfortable must be said, because the article reports it in three different ways.

• The abstract says: «resequence seven individuals of all four species» — seven individuals for each of the four species.
• The Methods section says: blood and tissue samples obtained from eight wild guanacos, eight vicuñas, eight domestic llamas and seven alpacas.
• The legend of Figure 2 indicates N = 8 individuals for both alpaca and llama on the ancestry-proportion axes; the section on demographic history says it used «6 (alpaca and llama) or 8 (vicuña and guanaco)» resequenced genomes for the MSMC analysis.
• The phylogenetic tree is built on 32 South American camelids in total (declared in the legend of Figure 1c for the ADMIXTURE analysis).

This factsheet does not choose for the authors which number is the correct one. It reports the three and notes that the order of magnitude is seven to eight individuals per species, which is the figure on which every confidence interval in the article depends. The exact breakdown per analysis: [TO BE CHECKED AGAINST THE FULL TEXT, supplementary tables S8 and S9].

Sample provenance, as declared: guanacos and vicuñas sampled across the whole distribution range of each species; for the de novo assemblies, one adult male guanaco from Putre (Chile) and one adult female vicuña from Lauca National Park (Chile), chosen as representatives of the northern subspecies Lama guanicoe cacsilensis and Vicugna vicugna mensalis; llama tissue from Putre and alpaca tissue from Temuco (Chile). For expression data, skin samples from 17 alpacas from China.

No European animal. No Italian alpaca, no European llama: the introgression proportions measured here are from South American and Chinese animals, and transferability to European herds is a reasonable assumption but not a measured one.

Tipo di studio e livello di evidenza

Studio genomico comparativo di popolazione. Non è un esperimento clinico, non ha gruppi assegnati, non ha esito clinico. È il confronto fra genomi di quattro specie per rispondere a quattro domande dichiarate: origine evolutiva, entità e direzione dell'introgressione, storia demografica, geni sotto selezione.

Livello di evidenza: come per la scheda CRL-E-02, è una categoria che le gerarchie cliniche non coprono. È studio sperimentale con approvazione etica e numerosità dichiarata nel vocabolario del §1.3 del capitolo 14 della KT — con la riserva sulla numerosità appena esposta.

La forza di questo lavoro non sta nella numerosità: sta nella convergenza di metodi indipendenti. La stessa conclusione — un'introgressione massiccia e unidirezionale dal lama all'alpaca — è raggiunta da cinque strade separate: albero neighbor-joining, ADMIXTURE, Treemix, statistica ABBA-BABA (D e f_d) e Local Ancestry Inference. Quando cinque metodi con assunzioni diverse danno la stessa risposta, la risposta non dipende dalle assunzioni di uno di essi. È il criterio con cui va giudicato uno studio genomico, e non è la numerosità.

Standard di reporting pertinente: ARRIVE 2.0 per la componente in vivo (prelievo di sangue e tessuto su animali selvatici e domestici) e per la parte sperimentale in vitro. Dichiarato dagli autori: no. È dichiarata l'approvazione etica — «All animals were handled and treated according to the alpaca guidelines approved by the Animal Ethics Committee [2017(55)] at Shanxi Agricultural University» — ma non la conformità a una checklist. Come sempre in questa collana, la griglia si usa retroattivamente per dire cosa manca.

Study type and level of evidence

Comparative population genomic study. It is not a clinical trial, it has no assigned groups, it has no clinical outcome. It is the comparison between the genomes of four species to answer four stated questions: evolutionary origin, extent and direction of introgression, demographic history, genes under selection.

Level of evidence: as with factsheet CRL-E-02, this is a category that clinical hierarchies do not cover. It is an experimental study with ethical approval and declared sample size in the vocabulary of §1.3 of chapter 14 of the KT — with the reservation on sample size just set out.

The strength of this work does not lie in sample size: it lies in the convergence of independent methods. The same conclusion — massive, unidirectional introgression from llama to alpaca — is reached by five separate routes: neighbor-joining tree, ADMIXTURE, Treemix, ABBA-BABA statistic (D and f_d) and Local Ancestry Inference. When five methods with different assumptions give the same answer, the answer does not depend on the assumptions of any one of them. This is the criterion by which a genomic study should be judged, not sample size.

Relevant reporting standard: ARRIVE 2.0 for the in vivo component (blood and tissue sampling on wild and domestic animals) and for the in vitro experimental part. Declared by the authors: no. Ethical approval is declared — «All animals were handled and treated according to the alpaca guidelines approved by the Animal Ethics Committee [2017(55)] at Shanxi Agricultural University» — but not compliance with a checklist. As always in this collection, the grid is used retroactively to say what is missing.

In breve

Questo è il lavoro che ha risolto, con i genomi interi, una questione discussa per due secoli: da chi discendono il lama e l'alpaca. La risposta conferma l'ipotesi che era già la più accreditata — il lama viene dal guanaco, l'alpaca viene dalla vigogna — ma nel farlo scopre qualcosa che nessuno aveva quantificato: il genoma di ogni alpaca analizzata contiene circa il 36% di ascendenza guanaco/lama, contro un 5% circa di ascendenza vigogna/alpaca nel genoma del lama. L'incrocio, insomma, c'è stato, è stato massiccio e squilibrato, e ha lasciato tracce che si possono datare: l'analisi della lunghezza dei segmenti introgressi colloca l'evento a circa 121 generazioni fa per l'alpaca e 115 per il lama, cioè — a seconda della durata di generazione che si assume — fra i 345 e i 727 anni fa. In quella finestra, nelle Ande, è successa una cosa sola: la conquista spagnola, cominciata nel 1529, con la caduta di Cuzco nel 1532.

Gli autori collegano i due fatti in modo esplicito e argomentato: dopo la conquista, con la scomparsa quasi totale delle mandrie, la morte di gran parte dei pastori e la fine della trasmissione orale delle pratiche di allevamento, le barriere fra le due specie domestiche cedettero. Nello stesso periodo, e con metodi indipendenti, si misura un'accelerazione del declino della dimensione efficace di popolazione a circa 110 generazioni fa.

Il resto dell'articolo è la mappa di che cosa quell'incrocio ha spostato: geni della qualità e del colore del vello, recettori olfattivi, un blocco di geni legato alla pressione arteriosa e all'adattamento all'ipossia, e — sul cromosoma X — il segno che la selezione ha fatto resistenza.

In brief

This is the work that resolved, with whole genomes, a question debated for two centuries: who do the llama and the alpaca descend from. The answer confirms the hypothesis that already had the most support — the llama comes from the guanaco, the alpaca comes from the vicuña — but in doing so it discovers something no one had quantified: the genome of every alpaca analyzed contains about 36% guanaco/llama ancestry, against about 5% vicuña/alpaca ancestry in the llama genome. Crossing, in short, did happen, was massive and unbalanced, and left traces that can be dated: analysis of the length of introgressed segments places the event at about 121 generations ago for the alpaca and 115 for the llama, that is — depending on the assumed generation length — between 345 and 727 years ago. In that window, in the Andes, one thing happened: the Spanish conquest, beginning in 1529, with the fall of Cuzco in 1532.

The authors link the two facts explicitly and with argument: after the conquest, with the almost total disappearance of herds, the death of most herders and the end of the oral transmission of husbandry practices, the barriers between the two domestic species gave way. In the same period, and with independent methods, an acceleration of the decline in effective population size around 110 generations ago is measured.

The rest of the article is the map of what that crossing moved: fleece quality and color genes, olfactory receptors, a gene block linked to blood pressure and hypoxia adaptation, and — on the X chromosome — the sign that selection pushed back.

Perché questo studio conta

Per una ragione che riguarda direttamente chi ha animali in stalla, e non solo chi fa filogenesi.

Perché toglie di mezzo la parola «purezza». In una specie in cui ogni individuo analizzato porta oltre il 30% di genoma di un altro genere, l'idea di un'alpaca «pura» come categoria biologica non ha oggetto. Non è un'opinione dell'associazione: è la conseguenza aritmetica di quel 36%, misurato con cinque metodi. Le distinzioni che contano in allevamento — tipo di vello, finezza, morfologia, registro genealogico — restano tutte in piedi; quella che cade è l'idea che corrispondano a linee genomiche separate.

E perché spiega dove sta la fibra. Fra i geni introgressi dal lama all'alpaca c'è FGF5, il gene del vello lungo; fra quelli passati dall'alpaca al lama ci sono TAGLN, legato allo sviluppo del follicolo nella pecora merino, e EDN3, che gli autori collegano al fenotipo bianco con occhi azzurri presente in entrambe le specie. I caratteri su cui si fonda il valore commerciale di questi animali stanno, in buona parte, nelle regioni scambiate fra le due specie. Chi seleziona, seleziona lì.

Why this study matters

For a reason that concerns directly anyone with animals in the barn, not only phylogeneticists.

Because it removes the word "purity" from the table. In a species in which every individual analyzed carries over 30% of the genome of another genus, the idea of a "pure" alpaca as a biological category has no referent. This is not an opinion of the association: it is the arithmetic consequence of that 36%, measured with five methods. The distinctions that matter in breeding — fleece type, fineness, morphology, pedigree register — all still stand; what falls is the idea that they correspond to separate genomic lineages.

And because it explains where the fiber lies. Among the genes introgressed from llama to alpaca is FGF5, the long-fleece gene; among those passed from alpaca to llama are TAGLN, linked to follicle development in the Merino sheep, and EDN3, which the authors link to the white phenotype with blue eyes present in both species. The traits on which the commercial value of these animals is founded lie, in large part, in the regions exchanged between the two species. Whoever selects, selects there.

Obiettivo

Dichiarato dagli autori alla fine dell'introduzione, in quattro punti: chiarire, per il lama e l'alpaca, (1) le origini evolutive; (2) livello, origine e distribuzione dell'introgressione in entrambe le specie; (3) la storia demografica dalla domesticazione a oggi, compresa la conquista spagnola; (4) i geni selezionati durante la domesticazione e potenzialmente per effetto dell'incrocio, e il loro significato per gli adattamenti delle due specie.

Domanda in forma PICO: non applicabile. È uno studio di genomica delle popolazioni; costruirgli addosso una PICO sarebbe un travestimento.

Objective

Declared by the authors at the end of the introduction, in four points: to clarify, for the llama and the alpaca, (1) their evolutionary origins; (2) the level, origin and distribution of introgression in both species; (3) the demographic history from domestication to the present, including the Spanish conquest; (4) the genes selected during domestication and potentially as an effect of crossing, and their significance for the adaptations of the two species.

Question in PICO form: not applicable. This is a population genomics study; forcing a PICO onto it would be a disguise.

Disegno dello studio

Tre passaggi, in questo ordine.

1. Assemblaggio de novo dei genomi di guanaco, vigogna e lama, con annotazione e trascrittomica, confrontati con il genoma di riferimento dell'alpaca già pubblicato.
2. Riesequenziamento di sette-otto individui per specie (vedi sopra), chiamata delle varianti e analisi di ascendenza con cinque metodi indipendenti.
3. Analisi funzionale su un sottoinsieme di geni candidati: espressione in cute di alpaca Suri e Huacaya, e coltura di melanociti di alpaca e di pecora in normossia e ipossia.

Nessuna assegnazione, nessun cieco, nessun follow-up — e, a differenza della scheda E-01, qui non sono mancanze: non esiste un disegno randomizzato per la domanda «da chi discende l'alpaca».

Study design

Three steps, in this order.

1. De novo assembly of the guanaco, vicuña and llama genomes, with annotation and transcriptomics, compared against the already published alpaca reference genome.
2. Resequencing of seven-eight individuals per species (see above), variant calling and ancestry analysis with five independent methods.
3. Functional analysis on a subset of candidate genes: expression in Suri and Huacaya alpaca skin, and culture of alpaca and sheep melanocytes under normoxia and hypoxia.

No allocation, no blinding, no follow-up — and, unlike factsheet E-01, here these are not shortcomings: there is no randomized design for the question "who does the alpaca descend from".

Metodi / Protocollo

Sequenziamento e assemblaggio. Piattaforme Illumina HiSeq 2000 e 2500; librerie paired-end con inserti di circa 170, 500 e 800 bp e librerie mate pair con inserti di 2, 5, 10, 20 e 40 kb. 335, 285 e 262 Gb di dati di alta qualità rispettivamente per guanaco, vigogna e lama. Assemblaggio con SOAPdenovo v2.04.

Qualità degli assemblaggi, con i numeri dichiarati:

Nota tecnica utile, e che vale in generale: l'assemblaggio del lama ha un N50 degli scaffold di 3,2 Mb, cioè più di quattro volte peggiore di quello del guanaco (14,6 Mb). I tre genomi non hanno la stessa qualità, e questo va tenuto presente quando si confrontano risultati che dipendono dalla contiguità. Il controllo BUSCO ha allineato oltre il 93,7% dei geni conservati di riferimento in tutti e tre gli assemblaggi, e le assemblee coprono oltre l'82% dei geni unici di ciascuna specie secondo l'allineamento degli unigeni da RNA-seq: sono valori buoni, e la percentuale di ripetizioni — fra il 27% e il 33% — è più bassa che nella maggior parte dei mammiferi, cosa che gli autori attribuiscono onestamente ai limiti del sequenziamento a letture corte, non a una peculiarità biologica.

Riesequenziamento: 35-49 Gb di dati grezzi per individuo, profondità media 16× dopo il filtraggio.

Variabilità genetica misurata (SNP e tasso di eterozigosi), per specie:

⚠️ Da leggere con attenzione: l'alpaca ha il tasso di eterozigosi più alto delle quattro specie (0,0030), più alto perfino di quello della vigogna selvatica. Non è un segno di buona salute genetica: è in buona parte l'effetto meccanico dell'introgressione, cioè della presenza nello stesso genoma di segmenti provenienti da due generi diversi. È un ottimo esempio di perché un singolo indice di diversità, letto da solo, può dire l'opposto della verità.

Analisi di ascendenza — i cinque metodi.

• Albero neighbor-joining su 5.901.447 SNP: due grandi linee con 100% di bootstrap, vigogna + alpaca da una parte (genere Vicugna), guanaco + lama dall'altra (genere Lama).
• ADMIXTURE: raggruppamento ottimale a k = 2 (validazione incrociata CV = 0,593 per k = 2, contro 0,620, 0,657 e 0,744 per k = 3, 4 e 5).
• Treemix, a conferma del flusso genico unidirezionale.
• Statistica ABBA-BABA: D per ((alpaca, vigogna), lama, battriano) = 0,565 ± 0,006 (z-test bilaterale, P < 0,01); D per ((guanaco, lama), vigogna, battriano) = −0,185 ± 0,006 (P < 0,01). La proporzione di introgressione stimata con la statistica f_d è ~38,9% dal lama all'alpaca e ~4,5% in direzione opposta.
• Local Ancestry Inference (LAI): ~94,5% del genoma del lama di origine guanaco/lama e ~36% del genoma dell'alpaca di origine guanaco/lama, con il 64% restante di origine vigogna/alpaca.

Datazione dell'evento: analisi della lunghezza dei segmenti di ascendenza introgressa (più il segmento è corto, più antico è l'incrocio, perché la ricombinazione lo ha spezzettato).

Storia demografica: PSMC per il tempo profondo, MSMC per il tempo intermedio, SNeP per le ultime ~200 generazioni, più l'analisi NeS della pendenza del declino.

Firme di selezione: combinazione di outlier di F_ST e XP-EHH (omozigosi estesa degli aplotipi), confrontando ciascuna specie domestica con il proprio antenato selvatico.

Validazione funzionale: qPCR e Western blot su cute di alpaca Suri e Huacaya, maschi e femmine; melanociti di alpaca e di pecora coltivati in normossia (21% O₂) e ipossia (13% O₂), con tre ripetizioni dell'esperimento (N = 3) e quattro repliche per le reazioni di qPCR.

Methods / Protocol

Sequencing and assembly. Illumina HiSeq 2000 and 2500 platforms; paired-end libraries with inserts of about 170, 500 and 800 bp and mate pair libraries with inserts of 2, 5, 10, 20 and 40 kb. 335, 285 and 262 Gb of high-quality data respectively for guanaco, vicuña and llama. Assembly with SOAPdenovo v2.04.

Assembly quality, with declared figures:

Useful technical note, of general validity: the llama assembly has a scaffold N50 of 3.2 Mb, that is more than four times worse than the guanaco's (14.6 Mb). The three genomes do not have the same quality, and this must be kept in mind when comparing results that depend on contiguity. The BUSCO check aligned over 93.7% of conserved reference genes in all three assemblies, and the assemblies cover over 82% of unigenes unique to each species according to RNA-seq unigene alignment: these are good values, and the repeat percentage — between 27% and 33% — is lower than in most mammals, which the authors honestly attribute to the limits of short-read sequencing, not to a biological peculiarity.

Resequencing: 35-49 Gb of raw data per individual, mean depth 16× after filtering.

Genetic variability measured (SNPs and heterozygosity rate), by species:

⚠️ To be read with attention: the alpaca has the highest heterozygosity rate of the four species (0.0030), higher even than that of the wild vicuña. This is not a sign of good genetic health: it is largely the mechanical effect of introgression, that is, the presence within the same genome of segments from two different genera. It is an excellent example of why a single diversity index, read on its own, can say the opposite of the truth.

Ancestry analysis — the five methods.

• Neighbor-joining tree on 5,901,447 SNPs: two major lineages with 100% bootstrap, vicuña + alpaca on one side (genus Vicugna), guanaco + llama on the other (genus Lama).
• ADMIXTURE: optimal clustering at k = 2 (cross-validation CV = 0.593 for k = 2, against 0.620, 0.657 and 0.744 for k = 3, 4 and 5).
• Treemix, confirming unidirectional gene flow.
• ABBA-BABA statistic: D for ((alpaca, vicuña), llama, Bactrian) = 0.565 ± 0.006 (two-tailed z-test, P < 0.01); D for ((guanaco, llama), vicuña, Bactrian) = −0.185 ± 0.006 (P < 0.01). The introgression proportion estimated with the f_d statistic is ~38.9% from llama to alpaca and ~4.5% in the opposite direction.
• Local Ancestry Inference (LAI): ~94.5% of the llama genome of guanaco/llama origin and ~36% of the alpaca genome of guanaco/llama origin, with the remaining 64% of vicuña/alpaca origin.

Dating of the event: analysis of the length of introgressed ancestry segments (the shorter the segment, the older the crossing, because recombination has broken it up).

Demographic history: PSMC for deep time, MSMC for intermediate time, SNeP for the last ~200 generations, plus NeS analysis of the decline slope.

Selection signatures: combination of F_ST outliers and XP-EHH (extended haplotype homozygosity), comparing each domestic species with its own wild ancestor.

Functional validation: qPCR and Western blot on Suri and Huacaya alpaca skin, males and females; alpaca and sheep melanocytes cultured under normoxia (21% O₂) and hypoxia (13% O₂), with three replicates of the experiment (N = 3) and four replicates for the qPCR reactions.

Risultati principali

A · Le origini: confermate, e con una precisazione

L'albero costruito su tutto il genoma colloca l'alpaca dentro il genere Vicugna con 100% di bootstrap, ma non dentro nessuna delle due sottospecie di vigogna: forma un gruppo a parte, con rami più lunghi. Quando si rimuovono dal genoma dell'alpaca tutte le regioni prive di ascendenza vigogna significativa, l'alpaca si sposta e va a raggrupparsi con Vicugna vicugna mensalis, la sottospecie settentrionale.

È il risultato più elegante dell'articolo, ed è anche un avvertimento metodologico: l'introgressione deforma l'albero filogenetico. Un'analisi condotta senza toglierla avrebbe concluso che l'alpaca è una linea a sé — che è, guarda caso, la conclusione a cui erano arrivati alcuni studi precedenti su pochi marcatori.

Per il lama: tutti i lama si raggruppano con i guanachi Lama guanicoe cacsilensis del Cile settentrionale, il che sostiene la domesticazione del lama a partire da quel pool genico. Gli autori segnalano inoltre, con onestà, che due campioni etichettati come cacsilensis (guanaco 22 e 24, raccolti a Paposo e Ovalle, rispettivamente 7° e 11° di latitudine più a sud degli altri) sono verosimilmente etichettati male e appartengono all'altra sottospecie.

Che cosa gli autori dichiarano di NON aver affrontato — ed è la frase che ogni divulgatore salta: «we did not attempt to address how many times each species was domesticated nor the likely precise origins of domestication». Quante volte l'alpaca sia stata domesticata, e dove esattamente, questo studio non lo dice, e per saperlo servirebbero più popolazioni attuali (soprattutto peruviane) e il DNA antico — che è precisamente il mestiere della scheda CRL-E-05.

B · L'introgressione: quanto, in che direzione, e su quali cromosomi

L'introgressione è dieci volte più forte in una direzione che nell'altra, ed è pervasiva e notevolmente costante fra individui: tutte le stime, su tutti gli alpaca sequenziati, superano il 30%. È questa costanza — non la media — che porta gli autori a concludere per un evento unico e consistente, non per una deriva di incroci sparsi.

Il cromosoma X racconta la storia opposta, ed è il dato più fine dell'articolo.

Sull'X è passato molto meno. Gli autori interpretano il dato come un costo adattativo dell'incrocio: la selezione ha eliminato i segmenti eterospecifici sul cromosoma sessuale più efficacemente che sugli autosomi — un fenomeno osservato anche fra bovini taurini e indicini, dove però gli antenati erano sottospecie, non specie di generi diversi, e il costo atteso era quindi molto minore. Segnalano anche la spiegazione alternativa (antichi sweep selettivi che abbassano la divergenza e rendono l'introgressione non rilevabile), senza sceglierne una: è il modo corretto di procedere.

Sull'X, otto regioni si distinguono: quattro ad alta introgressione (79 geni) e quattro a introgressione quasi nulla (116 geni). In queste otto regioni ci sono 53 geni con effetti noti deleteri o letali legati all'X nell'uomo o nel topo — 25 nelle regioni a bassa introgressione e 28 in quelle ad alta. Fra gli esempi citati, PTCHD1, le cui mutazioni deleterie sono descritte in pazienti maschi con disturbo dello spettro autistico legato all'X e/o disabilità intellettiva.

⚠️ Un avvertimento che questa scheda considera obbligatorio. «Geni con fenotipi deleteri legati all'X nell'uomo» non significa che quei geni causino qualcosa di analogo in alpaca e lama. Significa che quei geni, nella specie umana, sono annotati come tali — ed è così che gli autori li usano: come indicatore che la regione è funzionalmente sensibile. Non esiste in questo articolo alcuna osservazione clinica su alpaca. Chi trasformasse questo elenco in un sospetto diagnostico starebbe facendo un salto che gli autori non fanno.

C · La data, e come è stata ottenuta

Il numero di generazioni è misurato; gli anni no. Convertirlo in anni richiede la durata della generazione, che per lama e alpaca non è mai stata pubblicata. Gli autori lo dicono espressamente e mostrano il calcolo per quattro ipotesi, invece di sceglierne una in silenzio:

(Per il guanaco selvatico le stime pubblicate di durata della generazione variano fra 3 e 11 anni, con media di 7 anni; gli autori assumono che le specie domestiche abbiano generazioni più brevi delle selvatiche, come è comune sotto domesticazione.)

A questo punto il ragionamento si inverte, e va capito bene perché è la parte più contestabile dell'articolo. Gli autori osservano che l'unico evento noto in quella finestra è la conquista spagnola (inizio 1529, caduta di Cuzco 1532), e da lì deducono all'indietro che la durata della generazione di lama e alpaca sia di circa 4 anni. È un'inferenza dichiarata, plausibile e trasparente — ma è circolare se la si usa poi come prova indipendente della data. La datazione dell'evento e la stima della durata di generazione non sono due risultati indipendenti: sono lo stesso risultato letto in due direzioni.

Il sostegno indipendente esiste ed è altrove, ed è la ragione per cui l'ipotesi regge comunque: l'analisi NeS della pendenza del declino demografico identifica un'accelerazione del declino a circa 110 generazioni fa in entrambe le specie domestiche, ottenuta con un metodo che non usa la lunghezza dei segmenti introgressi. Due orologi diversi segnano la stessa ora.

D · La demografia: un declino che non si è fermato

• PSMC: espansione forte a partire da un milione di anni fa per vigogna e alpaca; per guanaco e lama l'espansione comincia solo dopo circa 400.000 anni fa, preceduta da un declino.
• Tutte e quattro le specie mostrano un forte declino da 100.000-200.000 anni fa fino a 25.000 anni fa, sovrapposto agli ultimi due grandi periodi glaciali.
• MSMC: declino concorde in tutte e quattro le specie fra 100.000 e 10.000 anni fa.
• SNeP (ultime ~200 generazioni): declino ulteriore e continuo, il più marcato nell'alpaca.

⚠️ Le due cifre di N_e dell'articolo si riferiscono a momenti diversi, e vanno lette insieme. La sezione Risultati riporta valori di dimensione efficace di popolazione (N_e) da 2.500 (lama) a 3.500 (alpaca considerata solo nelle porzioni di ascendenza vigogna), ma aggiunge che somigliano ai valori di mille anni fa nelle simulazioni MSMC delle due specie (testo pieno, Risultati, demografia): vanno quindi letti come la cifra di un passato di circa duecento generazioni fa, non come quella di oggi. La sezione Discussione porta lo stesso declino fino a «una dimensione efficace attuale inferiore a 1.000» (testo pieno, Discussione). La N_e attuale secondo questo studio è quindi inferiore a 1.000, senza una stima puntuale nel testo; i valori esatti per specie stanno nella Figura 2d e nei materiali supplementari, che qui non sono stati letti. La direzione è la stessa: la dimensione efficace di popolazione di lama e alpaca scende da duecento generazioni a oggi.

Perché questo interessa un allevatore e non solo un genetista. Una dimensione efficace di popolazione nell'ordine delle migliaia — o del migliaio — significa che la varianza genetica su cui si può selezionare è limitata e si sta restringendo, e che la consanguineità è un rischio strutturale della specie, non un incidente del singolo allevamento. È lo stesso problema che la scheda CRL-E-03 vede dal lato clinico quando riferisce difetti congeniti attribuiti a un pool genico ristretto.

E · I geni: dove sta la fibra, e dove sta l'ipossia

Geni introgressi identificati con entrambi i metodi (LAI e f_d), con la direzione dello scambio, come riportati dagli autori:

Gli autori riferiscono inoltre l'introgressione dal lama all'alpaca del cluster HoxD (HoxD3, HoxD4, HoxD10, HoxD12, HoxD13), di EVX2 e di OLA1, geni le cui mutazioni sono note per alterare lo sviluppo scheletrico, e ipotizzano che questo possa aver contribuito alla taglia maggiore dell'alpaca rispetto alla vigogna — differenza evidente, per esempio, nell'altezza al garrese. È un'ipotesi degli autori, non una dimostrazione.

Il blocco ANTXR2 / PRDM8 / FGF5 / C4orf22 merita due righe a parte, perché è il punto in cui genomica ed ecologia si toccano. È un blocco sintenico al locus umano HSA4q21, associato negli studi di associazione genome-wide alla regolazione della pressione arteriosa in popolazioni asiatiche ed europee. Gli autori misurano lungo la regione un crollo progressivo della diversità: θπ = 0,17 nel guanaco, 0,13 nel lama, 0,02 nell'alpaca — cioè l'alpaca è, in quella regione, quasi priva di SNP. Notano infine che Bacillus anthracis è noto essere presente nel suolo andino al di sotto dei 4.000 m di altitudine e assente al di sopra: un aplotipo che porta insieme un recettore dell'antrace, un gene del vello e una regione della pressione arteriosa, e che passa dal lama (specie di quote più basse) all'alpaca, è un'ipotesi di introgressione adattativa che si può raccontare. Gli autori la presentano come tale, e questa scheda non la promuove a fatto.

La verifica funzionale, con i suoi numeri. Le sequenze codificanti di FGF5 nella cute di alpaca Suri e Huacaya, maschi e femmine, sono risultate identiche a eccezione di due sostituzioni T→G alle posizioni 50 e 346. L'espressione dell'mRNA di FGF5 è risultata significativamente più alta nell'alpaca Suri, in entrambi i sessi (ANOVA a due vie, P < 0,001), risultato replicato al Western blot (P < 0,001). Sul versante ipossia, in melanociti di alpaca e di pecora coltivati a 13% di O₂ contro 21%, FGF5, ANTXR2 e C4orf22 sono risultati più espressi in ipossia, mentre PRDM8 è risultato significativamente meno espresso — insieme a SOX6 e a HIF1α (P < 0,05 e P < 0,01 rispettivamente; N = 3 esperimenti, 4 repliche per reazione).

⚠️ Su questo esperimento la cautela è d'obbligo: sono melanociti in coltura, di due specie, con tre ripetizioni. È una verifica di plausibilità del meccanismo, non una misura su animali in quota. Gli autori la presentano correttamente come un meccanismo di regolazione trascrizionale potenziale.

F · Le firme di domesticazione, e perché vanno lette con prudenza

La riga da guardare è l'ultima. Due metodi che cercano la stessa cosa — le tracce della selezione durante la domesticazione — concordano su 4 regioni su circa cento nell'alpaca e su 15 su circa centotrenta nel lama. Gli autori lo scrivono («overlap between the two methods was limited») e non lo nascondono. Significa che le liste di «geni della domesticazione» di questo articolo vanno prese come candidati da verificare, non come inventari. È l'esatto parallelo del punto sollevato nella scheda CRL-E-02 sul modo in cui si sceglie un gene candidato.

Le categorie funzionali riferite dagli autori: per l'alpaca, caratteri del vello, fecondità e sviluppo dei caratteri sessuali secondari, colore del mantello, ipossia, qualità della carne, olfatto; per il lama, controllo della pressione arteriosa, memoria spaziale e comportamento sociale, colore di cute e mantello, caratteri di performance e regolazione dell'appetito.

Un dato collaterale ma memorabile: il numero di recettori olfattivi funzionali dei camelidi sudamericani è risultato circa la metà di quello di bovino e cavallo (che ne hanno 566-602).

Interpretazione degli autori (attribuita, non fusa con i dati)

Gli autori concludono che l'alpaca discende dalla vigogna e il lama dal guanaco, e considerano risolta la questione delle origini filogenetiche. Su tutto il resto usano il condizionale, e va riportato com'è:

• il livello di introgressione nell'alpaca è «fra i più alti mai rilevati in un animale domestico»;
• la spiegazione «più plausibile» dell'evento contemporaneo nelle due specie è l'arrivo dei conquistadores, con il crollo delle mandrie, la mortalità di massa dei pastori e la rottura delle pratiche tradizionali di allevamento trasmesse solo oralmente; gli autori citano il fatto documentato che gli scrittori spagnoli non distinguevano il lama dall'alpaca nei loro scritti;
• «se questo evento spiega i segmenti introgressi», allora la durata di generazione delle due specie è di circa quattro anni;
• se l'introgressione sia stata deliberata o accidentale, «non è noto e merita ulteriore studio, soprattutto con i documenti storici»;
• l'introgressione non va letta solo in negativo: alcuni segmenti hanno aperto nuove strade alla selezione — colore, caratteri della fibra, adattamento all'alta quota — e nelle specie domestiche merita «un'interpretazione più sfumata» di quella che ha di solito in conservazione.

Main results

A · The origins: confirmed, with a caveat

The tree built on the whole genome places the alpaca within the genus Vicugna with 100% bootstrap, but not within either of the two vicuña subspecies: it forms a separate group, with longer branches. When all regions of the alpaca genome lacking significant vicuña ancestry are removed, the alpaca shifts and clusters with Vicugna vicugna mensalis, the northern subspecies.

This is the article's most elegant result, and it is also a methodological warning: introgression distorts the phylogenetic tree. An analysis conducted without removing it would have concluded that the alpaca is a lineage of its own — which, as it happens, is the conclusion that some earlier studies based on few markers had reached.

For the llama: all llamas cluster with the Lama guanicoe cacsilensis guanacos of northern Chile, supporting llama domestication from that gene pool. The authors also honestly note that two samples labeled cacsilensis (guanacos 22 and 24, collected in Paposo and Ovalle, respectively 7° and 11° of latitude further south than the others) are likely mislabeled and belong to the other subspecies.

What the authors declare they did NOT address — the sentence every popularization skips: «we did not attempt to address how many times each species was domesticated nor the likely precise origins of domestication». How many times the alpaca was domesticated, and exactly where, this study does not say, and to find out one would need more present-day populations (especially Peruvian ones) and ancient DNA — which is precisely the business of factsheet CRL-E-05.

B · Introgression: how much, in which direction, and on which chromosomes

Introgression is ten times stronger in one direction than the other, and it is pervasive and remarkably consistent across individuals: all estimates, across all sequenced alpacas, exceed 30%. It is this consistency — not the average — that leads the authors to conclude in favor of a single, consistent event, not a scatter of sporadic crosses.

The X chromosome tells the opposite story, and it is the article's finest datum.

Much less passed through on the X. The authors interpret the datum as an adaptive cost of crossing: selection eliminated heterospecific segments on the sex chromosome more effectively than on the autosomes — a phenomenon also observed between taurine and indicine cattle, where however the ancestors were subspecies, not species of different genera, and the expected cost was therefore much lower. They also note the alternative explanation (ancient selective sweeps lowering divergence and making introgression undetectable), without choosing between them: this is the correct way to proceed.

On the X, eight regions stand out: four with high introgression (79 genes) and four with almost no introgression (116 genes). In these eight regions there are 53 genes with known deleterious or lethal X-linked effects in humans or mice — 25 in the low-introgression regions and 28 in the high-introgression ones. Among the examples cited, PTCHD1, whose deleterious mutations are described in male patients with X-linked autism spectrum disorder and/or intellectual disability.

⚠️ A warning this factsheet considers mandatory. "Genes with deleterious X-linked phenotypes in humans" does not mean that these genes cause anything analogous in alpaca and llama. It means that these genes, in the human species, are annotated as such — and that is how the authors use them: as an indicator that the region is functionally sensitive. There is no clinical observation on alpacas in this article. Anyone turning this list into a diagnostic suspicion would be making a leap the authors do not make.

C · The date, and how it was obtained

The number of generations is measured; the years are not. Converting it into years requires the generation length, which for llama and alpaca has never been published. The authors say so explicitly and show the calculation for four hypotheses, instead of silently choosing one:

(For the wild guanaco, published estimates of generation length vary between 3 and 11 years, with a mean of 7 years; the authors assume that domestic species have shorter generations than wild ones, as is common under domestication.)

At this point the reasoning reverses, and it must be understood clearly why this is the most contestable part of the article. The authors observe that the only known event in that window is the Spanish conquest (beginning 1529, fall of Cuzco 1532), and from that they deduce backward that the generation length of llama and alpaca is about 4 years. This is a declared, plausible and transparent inference — but it is circular if it is then used as independent proof of the date. The dating of the event and the estimate of generation length are not two independent results: they are the same result read in two directions.

Independent support does exist, and it lies elsewhere, which is why the hypothesis holds up regardless: the NeS analysis of the demographic decline slope identifies an acceleration of the decline about 110 generations ago in both domestic species, obtained with a method that does not use the length of introgressed segments. Two different clocks show the same hour.

D · Demography: a decline that has not stopped

• PSMC: strong expansion starting one million years ago for vicuña and alpaca; for guanaco and llama the expansion only begins after about 400,000 years ago, preceded by a decline.
• All four species show a strong decline from 100,000-200,000 years ago to 25,000 years ago, overlapping with the last two major glacial periods.
• MSMC: concordant decline in all four species between 100,000 and 10,000 years ago.
• SNeP (last ~200 generations): further, continuing decline, most marked in the alpaca.

⚠️ The article's two N_e figures refer to different moments, and should be read together. The Results section reports effective population size (N_e) values from 2,500 (llama) to 3,500 (alpaca, considering only the vicuña-ancestry portions), but adds that they closely resemble the values of 1,000 years ago in the MSMC simulations of the two species (full text, Results, demographic history): they should therefore be read as the figure of a past some two hundred generations back, not of today. The Discussion section carries the same decline down to «a current effective size of less than 1000» (full text, Discussion). The current N_e according to this study is therefore below 1,000, with no point estimate in the text; the exact per-species values are in Figure 2d and the supplementary materials, which were not read here. The direction is the same: the effective population size of llama and alpaca has been falling from two hundred generations ago to the present.

Why this matters to a breeder and not only to a geneticist. An effective population size on the order of thousands — or of a thousand — means that the genetic variance available for selection is limited and shrinking, and that inbreeding is a structural risk of the species, not an accident of the individual herd. It is the same problem that factsheet CRL-E-03 sees from the clinical side when it reports congenital defects attributed to a restricted gene pool.

E · The genes: where the fiber lies, and where hypoxia lies

Introgressed genes identified with both methods (LAI and f_d), with the direction of the exchange, as reported by the authors:

The authors further report the introgression from llama to alpaca of the HoxD cluster (HoxD3, HoxD4, HoxD10, HoxD12, HoxD13), of EVX2 and of OLA1, genes whose mutations are known to alter skeletal development, and hypothesize that this may have contributed to the alpaca's larger size relative to the vicuña — a difference evident, for example, in withers height. This is a hypothesis of the authors, not a demonstration.

The ANTXR2 / PRDM8 / FGF5 / C4orf22 block deserves a separate note, because it is the point where genomics and ecology meet. It is a block syntenic to the human locus HSA4q21, associated in genome-wide association studies with blood pressure regulation in Asian and European populations. The authors measure a progressive collapse of diversity along the region: θπ = 0.17 in guanaco, 0.13 in llama, 0.02 in alpaca — that is, the alpaca is, in that region, almost devoid of SNPs. They finally note that Bacillus anthracis is known to be present in Andean soil below 4,000 m of altitude and absent above it: a haplotype that carries together an anthrax receptor, a fleece gene and a blood-pressure region, and that passes from the llama (a lower-altitude species) to the alpaca, is a tellable hypothesis of adaptive introgression. The authors present it as such, and this factsheet does not promote it to fact.

Functional verification, with its numbers. The coding sequences of FGF5 in Suri and Huacaya alpaca skin, males and females, were found identical except for two T→G substitutions at positions 50 and 346. FGF5 mRNA expression was found significantly higher in Suri alpaca, in both sexes (two-way ANOVA, P < 0.001), a result replicated by Western blot (P < 0.001). On the hypoxia side, in alpaca and sheep melanocytes cultured at 13% O₂ versus 21%, FGF5, ANTXR2 and C4orf22 were found more expressed under hypoxia, while PRDM8 was found significantly less expressed — together with SOX6 and HIF1α (P < 0.05 and P < 0.01 respectively; N = 3 experiments, 4 replicates per reaction).

⚠️ Caution is mandatory on this experiment: these are melanocytes in culture, of two species, with three replicates. It is a plausibility check of the mechanism, not a measurement on animals at altitude. The authors correctly present it as a potential transcriptional regulation mechanism.

F · The domestication signatures, and why they must be read with caution

The line to look at is the last one. Two methods searching for the same thing — the traces of selection during domestication — agree on 4 regions out of about a hundred in the alpaca and on 15 out of about a hundred and thirty in the llama. The authors write this («overlap between the two methods was limited») and do not hide it. It means that this article's lists of "domestication genes" should be taken as candidates to be verified, not as inventories. It is the exact parallel of the point raised in factsheet CRL-E-02 about how a candidate gene is chosen.

The functional categories reported by the authors: for the alpaca, fleece traits, fertility and development of secondary sexual characteristics, coat color, hypoxia, meat quality, olfaction; for the llama, blood pressure control, spatial memory and social behavior, skin and coat color, performance traits and appetite regulation.

A collateral but memorable datum: the number of functional olfactory receptors in South American camelids was found to be about half that of cattle and horse (which have 566-602).

Authors' interpretation (attributed, not merged with the data)

The authors conclude that the alpaca descends from the vicuña and the llama from the guanaco, and regard the question of phylogenetic origins as settled. On everything else they use the conditional, and it must be reported as such:

• the level of introgression in the alpaca is «among the highest ever detected in a domestic animal»;
• the «most plausible» explanation for the event being contemporaneous in the two species is the arrival of the conquistadores, with the collapse of herds, mass mortality of herders and the breakdown of traditional husbandry practices transmitted only orally; the authors cite the documented fact that Spanish writers did not distinguish llama from alpaca in their writings;
• «if this event explains the introgressed segments», then the generation length of the two species is about four years;
• whether the introgression was deliberate or accidental, «is not known and deserves further study, especially with historical documents»;
• introgression should not be read only negatively: some segments opened new avenues for selection — color, fiber traits, high-altitude adaptation — and in domestic species it deserves «a more nuanced interpretation» than it usually receives in conservation.

Implicazioni pratiche

1. La categoria «alpaca pura» non ha un correlato genomico. Ogni alpaca analizzata porta oltre il 30% di ascendenza guanaco/lama. Chi vende, compra o registra animali può continuare a usare registri, tipi e standard morfologici — ma non può fondarli su una purezza di specie che il genoma non mostra.
2. La fibra sta nelle regioni scambiate. FGF5 (vello lungo) è arrivato all'alpaca dal lama; TAGLN, legato allo sviluppo del follicolo nella merino, è passato dall'alpaca al lama. Chi seleziona per la fibra lavora, senza saperlo, su un'eredità interspecifica.
3. Il bianco con occhi azzurri ha un gene candidato, ed è un problema di benessere. Gli autori collegano il fenotipo a EDN3, introgresso dall'alpaca al lama, e ricordano che nelle alpaca e nei lama bianchi con occhi azzurri la sordità è frequente. Per chi seleziona il mantello bianco, e per chi valuta un animale che «non risponde», è un'informazione clinica.
4. La dimensione efficace di popolazione è piccola e in calo. La consanguineità in questa specie non è un rischio del singolo allevamento: è una condizione di partenza. Vale per l'Europa, dove la base di riproduttori è ancora più stretta (vedi il dato LAREU citato nella scheda CRL-E-01: 7,7 animali per proprietario in media).
5. Un dato di terminologia che fa perdere letteratura. L'alpaca è stata riclassificata da Lama pacos a Vicugna pacos proprio in conseguenza di questa linea di ricerca. Chi interroga le banche dati solo con il nome nuovo non trova la letteratura anteriore — ed è lo stesso avvertimento già dato nella scheda CRL-E-03.

Practical implications

1. The category "pure alpaca" has no genomic correlate. Every alpaca analyzed carries over 30% guanaco/llama ancestry. Those who sell, buy or register animals may continue to use registries, types and morphological standards — but cannot ground them in a species purity the genome does not show.
2. The fiber lies in the exchanged regions. FGF5 (long fleece) arrived in the alpaca from the llama; TAGLN, linked to follicle development in Merino, passed from alpaca to llama. Those selecting for fiber are working, without knowing it, on interspecific inheritance.
3. White with blue eyes has a candidate gene, and it is a welfare issue. The authors link the phenotype to EDN3, introgressed from alpaca to llama, and note that in white, blue-eyed alpacas and llamas deafness is frequent. For those selecting for white coat, and for those evaluating an animal that "does not respond," this is clinical information.
4. Effective population size is small and declining. Inbreeding in this species is not a risk of the individual herd: it is a starting condition. This holds for Europe, where the breeding base is even narrower (see the LAREU figure cited in factsheet CRL-E-01: 7.7 animals per owner on average).
5. A terminology fact that loses literature. The alpaca was reclassified from Lama pacos to Vicugna pacos precisely as a consequence of this line of research. Anyone searching databases only with the new name will not find the earlier literature — the same warning already given in factsheet CRL-E-03.

Limiti

Limiti dichiarati dagli autori (ed è un elenco vero, non di forma): non hanno affrontato quante volte ciascuna specie sia stata domesticata né le origini precise della domesticazione, per cui servirebbero più popolazioni attuali — soprattutto peruviane — e il DNA antico; è possibile che segmenti introgressi residui non rilevabili con LAI abbiano influenzato la forma dell'albero e spieghino la separazione dell'alpaca da V. v. mensalis; la popolazione di mensalis effettivamente domesticata potrebbe non essere rappresentata dai loro campioni; alcuni SNP delle famiglie dei recettori olfattivi potrebbero essere artefatti di allineamento errato; la durata di generazione delle due specie domestiche non è disponibile in letteratura.

Limiti aggiunti dalla redazione:

• La numerosità è riportata in tre modi diversi nello stesso articolo (sette per specie nell'abstract; otto, otto, otto e sette nei metodi; N = 8 in Figura 2 e sei o otto nell'analisi MSMC). Nessuna di queste cifre è giustificata da un calcolo. Sette-otto individui per specie sono pochi per stimare proporzioni di ascendenza a livello di popolazione, e la costanza fra individui — che gli autori usano come argomento — è misurata su quei sette-otto.
• Nessun animale europeo. Le proporzioni di introgressione valgono per animali sudamericani; che un'alpaca italiana abbia lo stesso 36% è estremamente probabile ma non misurato in questo studio.
• La datazione dipende da un parametro che lo studio stesso ricava dalla datazione. Il ragionamento è dichiarato e trasparente, ma va detto: la data «conquista spagnola» e la durata di generazione «quattro anni» non sono due prove indipendenti.
• Le tre assemblee non hanno la stessa qualità (N50 degli scaffold da 14,6 Mb del guanaco a 3,2 Mb del lama), e il confronto fra genomi di contiguità diversa può introdurre asimmetrie che l'articolo non quantifica.
• La sovrapposizione fra i due metodi di firma di selezione è bassissima (4 regioni su ~100 nell'alpaca), il che rende le liste di geni della domesticazione dei candidati e non un risultato consolidato.
• Le funzioni geniche sono quasi tutte annotate su altre specie — uomo, topo, bovino, pecora merino, bovino Hanwoo. Nessuna è stata dimostrata in Vicugna pacos in questo lavoro, con l'eccezione dell'espressione differenziale di FGF5 fra Suri e Huacaya.
• La validazione funzionale è in vitro, su melanociti, con N = 3.
• La N_e attuale non ha una stima puntuale: 2.500-3.500 è il valore di circa mille anni fa (Risultati), «meno di 1.000» quello di oggi (Discussione), come spiegato sopra.
• Nessuna misura fenotipica sugli animali sequenziati: non un micron di diametro di fibra, non un peso di vello. È la stessa lacuna della scheda CRL-E-02, ed è ricorrente in questo campo.

Finanziamento: dichiarato ed esteso — progetti di ricerca agro-scientifica di interesse pubblico della Cina (n. 201303119), Strategic Priority Program dell'Accademia Cinese delle Scienze (XDB31000000), National Natural Science Foundation of China (31821001), Cardiff University e CAS President's International Fellowship Initiative, Young Sanjin Scholars, Shanxi Agricultural University (CXTD201201), FONDECYT/CONICYT del Cile (1140785), progetto «Sanjiangyuan National Park Animal Genome Program» e Second Tibetan Plateau Scientific Expedition (2019QZKK0501). Tutti finanziatori pubblici; nessuno sponsor industriale.

Conflitti di interesse: dichiarati assenti.

Disponibilità dei dati: dichiarata con identificatori verificabili — BioProject NCBI PRJNA427832, PRJNA421373, PRJNA427644; progetto whole-genome shotgun dei camelidi sudamericani PRJNA612032; dati RNA-seq dell'alpaca PRJNA636766.

🔗 Un collegamento verificabile fra due schede di questo fascicolo

PRJNA612032 è lo stesso identificatore di progetto da cui la scheda CRL-E-02 dichiara che provengono le sei vigogne e i sei guanachi usati come controlli selvatici dallo studio di Camerino su TRPV3. In altre parole: i genomi selvatici che hanno permesso al gruppo italiano di filtrare le varianti candidate del fenotipo Suri sono quelli depositati da questo lavoro. È un esempio concreto di che cosa significa depositare i dati grezzi con un identificatore — e di perché il §7 del capitolo 14 lo considera un criterio di qualità e non una formalità amministrativa.

Limitations

Limitations declared by the authors (and it is a genuine list, not a formal one): they did not address how many times each species was domesticated nor the precise origins of domestication, for which more present-day populations — especially Peruvian ones — and ancient DNA would be needed; it is possible that undetectable residual introgressed segments with LAI influenced the shape of the tree and explain the separation of the alpaca from V. v. mensalis; the mensalis population actually domesticated may not be represented by their samples; some SNPs of the olfactory receptor families may be artifacts of misalignment; the generation length of the two domestic species is not available in the literature.

Limitations added by the editorial team:

• Sample size is reported in three different ways in the same article (seven per species in the abstract; eight, eight, eight and seven in the methods; N = 8 in Figure 2 and six or eight in the MSMC analysis). None of these figures is justified by a calculation. Seven-eight individuals per species are few for estimating population-level ancestry proportions, and the consistency across individuals — which the authors use as an argument — is measured on those seven-eight.
• No European animal. The introgression proportions hold for South American animals; whether an Italian alpaca has the same 36% is extremely likely but not measured in this study.
• The dating depends on a parameter the study itself derives from the dating. The reasoning is declared and transparent, but it must be said: the "Spanish conquest" date and the "four years" generation length are not two independent proofs.
• The three assemblies do not have the same quality (scaffold N50 from 14.6 Mb for the guanaco to 3.2 Mb for the llama), and comparing genomes of different contiguity can introduce asymmetries the article does not quantify.
• The overlap between the two selection-signature methods is very low (4 regions out of ~100 in the alpaca), which makes the domestication-gene lists candidates and not a consolidated result.
• Gene functions are almost all annotated in other species — human, mouse, cattle, Merino sheep, Hanwoo cattle. None was demonstrated in Vicugna pacos in this work, with the exception of the differential expression of FGF5 between Suri and Huacaya.
• Functional validation is in vitro, on melanocytes, with N = 3.
• The current N_e has no point estimate: 2,500-3,500 is the value of about a thousand years ago (Results), "below 1,000" the value of today (Discussion), as explained above.
• No phenotypic measurement on the sequenced animals: not a micron of fiber diameter, not a fleece weight. It is the same gap as factsheet CRL-E-02, and it recurs throughout this field.

Funding: declared and extensive — public-interest agro-scientific research projects of China (no. 201303119), Strategic Priority Program of the Chinese Academy of Sciences (XDB31000000), National Natural Science Foundation of China (31821001), Cardiff University and CAS President's International Fellowship Initiative, Young Sanjin Scholars, Shanxi Agricultural University (CXTD201201), FONDECYT/CONICYT of Chile (1140785), the "Sanjiangyuan National Park Animal Genome Program" project and the Second Tibetan Plateau Scientific Expedition (2019QZKK0501). All public funders; no industry sponsor.

Conflicts of interest: declared none.

Data availability: declared with verifiable identifiers — NCBI BioProject PRJNA427832, PRJNA421373, PRJNA427644; South American camelid whole-genome shotgun project PRJNA612032; alpaca RNA-seq data PRJNA636766.

🔗 A verifiable link between two factsheets in this issue

PRJNA612032 is the same project identifier from which factsheet CRL-E-02 declares that the six vicuñas and six guanacos used as wild controls by the Camerino group's study on TRPV3 come from. In other words: the wild genomes that allowed the Italian group to filter candidate variants for the Suri phenotype are the ones deposited by this work. It is a concrete example of what it means to deposit raw data with an identifier — and of why §7 of chapter 14 regards it as a quality criterion and not an administrative formality.

Cosa manca rispetto allo standard di reporting (ARRIVE 2.0, applicato retroattivamente)

• La numerosità non è giustificata da alcun calcolo, e — cosa più grave — non è riportata in modo coerente nelle diverse sezioni dell'articolo.
• Mancano le caratteristiche degli animali nel testo principale: età, sesso, condizione, provenienza individuale sono rinviate alle tabelle supplementari.
• Mancano criteri di inclusione ed esclusione dei campioni.
• Manca qualunque misura fenotipica sugli animali sequenziati.
• Manca la dichiarazione di conformità a una checklist (l'approvazione etica c'è, la checklist no).
• Presente e ben fatto, invece, e va detto con la stessa nettezza: approvazione etica con numero di protocollo, deposito completo dei dati grezzi con cinque identificatori di progetto, descrizione dettagliata di software e versioni, convergenza dichiarata di cinque metodi indipendenti, limiti discussi in modo sostanziale, finanziamento e conflitti dichiarati, e storia della revisione paritaria pubblicata come file supplementare.

◆ COMMENTO DELL'ASSOCIAZIONE

Camelid Research Studies APS — redazione CRL · 31/08/2026

Il 36% è la proporzione di genoma di ascendenza guanaco/lama in ogni alpaca analizzata: un'alpaca moderna, geneticamente, è per più di un terzo l'altro genere.

Il 36% è la proporzione di genoma di ascendenza guanaco/lama in ogni alpaca analizzata. Detto in una frase che si può ripetere a voce: un'alpaca moderna, geneticamente, è per più di un terzo l'altro genere. Non è una curiosità tassonomica. È la ragione per cui la parola «purezza», in questa specie, non ha un referente biologico; è la ragione per cui l'alpaca ha l'eterozigosi più alta delle quattro specie senza per questo stare bene; ed è la ragione per cui i geni che interessano un allevatore — FGF5 per la lunghezza della fibra, TAGLN per il follicolo, EDN3 per il bianco con occhi azzurri — stanno nelle regioni scambiate.

Il 121,17 è invece il numero da cui si impara come si legge un articolo. È il tempo dell'incrocio in generazioni, con un intervallo di confidenza al 95% stretto (113,99-130,04). È misurato. Gli anni, invece, non sono misurati: dipendono dalla durata di una generazione di alpaca, che nessuno ha mai pubblicato. Gli autori mostrano quattro conversioni possibili — da 364 a 727 anni fa — e poi fanno una cosa che va vista bene: osservano che nella finestra c'è la conquista spagnola e da lì deducono la durata della generazione. È un ragionamento legittimo e dichiarato; ma la data e la durata della generazione non sono due prove indipendenti, e chi ripete «l'incrocio risale alla conquista spagnola, che è la prova che le alpaca fanno una generazione ogni quattro anni» sta chiudendo un cerchio su sé stesso. Quello che regge davvero l'ipotesi è un secondo orologio, che gli autori mettono in campo e che quasi nessuno cita: l'accelerazione del declino demografico a 110 generazioni fa, stimata con un metodo che non guarda i segmenti introgressi. Due metodi indipendenti che segnano la stessa ora sono un argomento; un metodo che si tara su sé stesso non lo è.

Perché questa scheda ci sembra necessaria per un veterinario italiano, e non solo per un genetista. Perché l'introgressione spiega, in un colpo solo, tre cose che in ambulatorio si vedono e non si spiegano: che alpaca e lama si incrociano e fanno figli fertili (e quindi che ogni «divieto naturale» è culturale, non biologico); che la base genetica su cui si può selezionare in Europa è strettissima e la consanguineità è un rischio strutturale; e che bianco più occhi azzurri è associato a sordità in entrambe le specie, con un gene candidato che ha un nome. Chi valuta un animale che «non risponde al richiamo» dovrebbe avere questa riga in testa prima di attribuirlo al carattere.

La cosa che vogliamo dire come associazione, però, è un'altra, e riguarda i dati. Questo articolo ha depositato tutto: cinque identificatori di progetto NCBI, software con le versioni, storia della revisione pubblicata. Uno di quei depositi — PRJNA612032 — è esattamente ciò che quattro anni dopo ha permesso a un gruppo dell'Università di Camerino, senza finanziamento dedicato, di avere sei vigogne e sei guanachi come controlli selvatici per cercare il gene del Suri (scheda CRL-E-02). Un laboratorio italiano ha potuto fare genomica comparata su specie selvatiche andine senza mandare nessuno sulle Ande, perché ventuno autori fra Cina, Cile, Perù, Danimarca, Stati Uniti, Italia e Galles avevano messo i dati grezzi in un archivio pubblico. La scienza aperta, qui, non è un principio: è la ragione per cui esiste la scheda E-02.

Che cosa non concludiamo. Non concludiamo che l'incrocio sia avvenuto nel 1532: la datazione è in generazioni, e la conversione in anni dipende da un parametro non misurato. Non concludiamo che il 36% valga per le alpaca europee: nessun animale europeo è stato sequenziato in questo studio, e per quanto sia probabile, «probabile» non è «misurato». Non concludiamo che i geni elencati causino i caratteri a cui sono associati: quasi tutte le funzioni sono annotate su uomo, topo, bovino o pecora, e l'unica dimostrazione fatta qui in Vicugna pacos è l'espressione differenziale di FGF5 fra Suri e Huacaya. E sulla dimensione efficace attuale delle due specie sappiamo solo che è inferiore a 1.000: l'articolo non dà una stima puntuale, e il 2.500-3.500 dei Risultati è il valore di circa mille anni fa.

Una nota di raccordo che vale il fascicolo intero. La scheda CRL-E-02 propone TRPV3 come gene candidato del fenotipo Suri; questa scheda mostra FGF5 — introgresso dal lama — significativamente più espresso nella cute delle alpaca Suri (P < 0,001). Non sono due affermazioni in conflitto, e chi le legge come tali sbaglia: il lavoro di Camerino cita esplicitamente FGF5 fra i precedenti, ricordando che vi era già stato descritto un codone di stop prematuro (p.Arg167*) e che studi genome-wide vi avevano trovato segnali di selezione. La lunghezza della fibra dell'alpaca è, con ogni probabilità, un carattere a più geni, e questo fascicolo ne mette due sul tavolo con due metodi diversi. Chi vende un test per «il gene del Suri» deve spiegare quale dei due, e perché.

Valutazione a stelle

Asse	Voto	Motivazione
A · Qualità scientifica	★★★☆☆	Studio sperimentale con approvazione etica e numero di protocollo dichiarati, numerosità dichiarata, deposito completo dei dati grezzi con cinque identificatori, software e versioni descritti, convergenza di cinque metodi indipendenti sulla conclusione principale, limiti discussi in modo sostanziale, finanziamento e conflitti dichiarati, storia della revisione paritaria pubblicata. La quarta stella richiederebbe coerenza con l'Essential 10 di ARRIVE 2.0, e mancano: giustificazione della numerosità, caratteristiche degli animali nel testo, criteri di esclusione — e soprattutto la numerosità è riportata in tre modi diversi nello stesso articolo (7 nell'abstract, 8/8/8/7 nei metodi, 6 o 8 nell'analisi MSMC). Con lo stesso metro applicato alla scheda E-02, il voto è lo stesso.
B · Utilità clinica	★★☆☆☆	Suggerisce ipotesi da verificare — geni candidati per fibra, colore e sordità, e un quadro di consanguineità strutturale — ma non offre oggi alcuna procedura applicabile: nessun test disponibile, nessuna raccomandazione riproduttiva sostenibile su questa base. L'unico elemento immediatamente utile in ambulatorio (bianco + occhi azzurri → sospetto di sordità) era già noto per altra via. Vincolo incrociato rispettato con margine. Nessuna dose, via o protocollo è coinvolto.
C · Valore didattico	★★★★★	È il riferimento di prima scelta sull'origine di lama e alpaca: chi deve imparare quella cosa parte da qui, e la letteratura successiva ci torna. Insegna inoltre quattro cose di metodo che valgono ben oltre i camelidi: che l'introgressione deforma l'albero filogenetico (l'alpaca cambia posto quando si rimuovono i segmenti introgressi); che un tempo misurato in generazioni non è un tempo misurato in anni; che due metodi che cercano la stessa selezione possono concordare su quattro regioni su cento; e che un tasso di eterozigosi alto può essere l'effetto di un incrocio, non un segno di salute. Questa collana prevede esplicitamente che C possa valere cinque con A a tre: è questo il caso.
D · Rilevanza per la specie	★★★★☆	Studio su tutte e quattro le specie di camelidi sudamericani, con campionamento delle due specie selvatiche esteso all'intero areale, e con un risultato che riguarda ogni alpaca esistente, comprese quelle italiane. Non arriva a cinque stelle perché nessun animale europeo è stato sequenziato, perché la numerosità per specie è di sette-otto individui e perché i dati non sono direttamente trasferibili a una decisione di allevamento europea: dicono da dove viene la specie, non che cosa farne domani.

Nessun voto complessivo, per costruzione. A = 3 e C = 5 non sono una contraddizione: un lavoro con sette animali per specie può essere il riferimento didattico di un intero campo, e questo lo è.

Bottom line. Si può concludere che il lama discende dal guanaco e l'alpaca dalla vigogna (sottospecie Vicugna vicugna mensalis, una volta rimossi i segmenti introgressi), con 100% di bootstrap su 5.901.447 SNP; che il genoma dell'alpaca contiene circa il 36% di ascendenza guanaco/lama (38,9% con la statistica f_d) contro circa il 5% di ascendenza vigogna/alpaca nel lama, misurati con cinque metodi indipendenti su sette-otto individui per specie; che l'introgressione è significativamente più bassa sul cromosoma X che sugli autosomi (P < 0,001); e che i segmenti scambiati contengono geni del vello (FGF5, TAGLN), del colore (EDN3, ASIP, CORIN) e una regione sintenica al locus umano della pressione arteriosa. Non si può concludere che l'evento risalga a una data del calendario: la stima è di 121,17 generazioni (IC 95% 113,99-130,04) per l'alpaca, e la sua conversione in anni dipende da una durata di generazione che nessuno ha misurato e che gli autori deducono, in modo dichiarato ma circolare, dalla data della conquista stessa.

What is missing relative to the reporting standard (ARRIVE 2.0, applied retroactively)

• Sample size is not justified by any calculation, and — more seriously — is not reported consistently across the different sections of the article.
• Animal characteristics are missing from the main text: age, sex, condition, individual provenance are relegated to supplementary tables.
• Inclusion and exclusion criteria for the samples are missing.
• Any phenotypic measurement on the sequenced animals is missing.
• A declaration of compliance with a checklist is missing (ethical approval is present, the checklist is not).
• Present and well done, however, and it must be said with equal clarity: ethical approval with protocol number, complete deposit of raw data with five project identifiers, detailed description of software and versions, declared convergence of five independent methods, substantially discussed limitations, declared funding and conflicts, and published peer review history as a supplementary file.

◆ ASSOCIATION COMMENT

Camelid Research Studies APS — CRL editorial team · 31/08/2026

The 36% is the proportion of guanaco/llama-ancestry genome in every alpaca analyzed: a modern alpaca, genetically, is more than a third the other genus.

The 36% is the proportion of guanaco/llama-ancestry genome in every alpaca analyzed. Said in a sentence that can be repeated out loud: a modern alpaca, genetically, is more than a third the other genus. This is not a taxonomic curiosity. It is the reason why the word "purity," in this species, has no biological referent; it is the reason why the alpaca has the highest heterozygosity of the four species without this meaning it is doing well; and it is the reason why the genes that matter to a breeder — FGF5 for fiber length, TAGLN for the follicle, EDN3 for white with blue eyes — lie in the exchanged regions.

The 121.17, instead, is the number from which one learns how to read an article. It is the time of the crossing in generations, with a tight 95% confidence interval (113.99-130.04). It is measured. The years, instead, are not measured: they depend on the length of an alpaca generation, which no one has ever published. The authors show four possible conversions — from 364 to 727 years ago — and then do something that must be looked at closely: they observe that the Spanish conquest falls within the window and from that deduce the generation length. This is a legitimate and declared piece of reasoning; but the date and the generation length are not two independent proofs, and anyone who repeats "the crossing dates to the Spanish conquest, which is the proof that alpacas have a generation every four years" is closing a circle on themselves. What truly supports the hypothesis is a second clock, which the authors bring into play and which almost no one cites: the acceleration of the demographic decline 110 generations ago, estimated with a method that does not look at introgressed segments. Two independent methods showing the same hour are an argument; a method calibrated on itself is not.

Why we think this factsheet is necessary for an Italian veterinarian, and not only for a geneticist. Because introgression explains, in one stroke, three things seen in the clinic and not otherwise explained: that alpaca and llama cross and produce fertile offspring (and therefore that any "natural prohibition" is cultural, not biological); that the genetic base available for selection in Europe is extremely narrow and inbreeding is a structural risk; and that white plus blue eyes is associated with deafness in both species, with a candidate gene that has a name. Anyone evaluating an animal that "does not respond to being called" should have this line in mind before attributing it to temperament.

What we want to say as an association, however, is something else, and it concerns the data. This article deposited everything: five NCBI project identifiers, software with versions, published review history. One of those deposits — PRJNA612032 — is exactly what, four years later, allowed a group at the University of Camerino, without dedicated funding, to have six vicuñas and six guanacos as wild controls to search for the Suri gene (factsheet CRL-E-02). An Italian laboratory was able to do comparative genomics on wild Andean species without sending anyone to the Andes, because twenty-one authors across China, Chile, Peru, Denmark, the United States, Italy and Wales had put the raw data in a public archive. Open science, here, is not a principle: it is the reason factsheet E-02 exists.

What we do not conclude. We do not conclude that the crossing occurred in 1532: the dating is in generations, and the conversion to years depends on an unmeasured parameter. We do not conclude that the 36% holds for European alpacas: no European animal was sequenced in this study, and however probable it may be, "probable" is not "measured." We do not conclude that the listed genes cause the traits they are associated with: almost all functions are annotated in human, mouse, cattle or sheep, and the only demonstration made here in Vicugna pacos is the differential expression of FGF5 between Suri and Huacaya. On the current effective size of the two species we know only that it is below 1,000: the article gives no point estimate, and the 2,500-3,500 of the Results is the value of about a thousand years ago.

A cross-reference note worth the whole issue. Factsheet CRL-E-02 proposes TRPV3 as the candidate gene of the Suri phenotype; this factsheet shows FGF5 — introgressed from the llama — significantly more expressed in Suri alpaca skin (P < 0.001). These are not two conflicting statements, and anyone reading them as such is mistaken: the Camerino group's work explicitly cites FGF5 among the precedents, recalling that a premature stop codon (p.Arg167*) had already been described there and that genome-wide studies had found selection signals there. Alpaca fiber length is, in all likelihood, a polygenic trait, and this issue puts two of those genes on the table with two different methods. Whoever sells a test for "the Suri gene" must explain which of the two, and why.

Star rating

Axis	Score	Rationale
A · Scientific quality	★★★☆☆	Experimental study with declared ethical approval and protocol number, declared sample size, complete deposit of raw data with five identifiers, software and versions described, convergence of five independent methods on the main conclusion, substantially discussed limitations, declared funding and conflicts, published peer review history. The fourth star would require compliance with ARRIVE 2.0's Essential 10, and missing are: sample size justification, animal characteristics in the text, exclusion criteria — and above all sample size is reported in three different ways in the same article (7 in the abstract, 8/8/8/7 in the methods, 6 or 8 in the MSMC analysis). Applying the same yardstick as factsheet E-02, the score is the same.
B · Clinical utility	★★☆☆☆	Suggests hypotheses to verify — candidate genes for fiber, color and deafness, and a picture of structural inbreeding — but offers no applicable procedure today: no available test, no sustainable breeding recommendation on this basis. The only element immediately useful in the clinic (white + blue eyes → suspicion of deafness) was already known through other means. Cross-axis constraint respected with margin. No dose, route or protocol is involved.
C · Educational value	★★★★★	This is the first-choice reference on the origin of llama and alpaca: anyone who needs to learn that starts here, and later literature returns to it. It also teaches four methodological lessons that hold well beyond camelids: that introgression distorts the phylogenetic tree (the alpaca changes place once introgressed segments are removed); that a time measured in generations is not a time measured in years; that two methods searching for the same selection can agree on four regions out of a hundred; and that a high heterozygosity rate can be the effect of a crossing, not a sign of health. This collection explicitly allows C to score five with A at three: this is that case.
D · Relevance to the species	★★★★☆	Study of all four South American camelid species, with sampling of the two wild species covering the entire range, and with a result that concerns every existing alpaca, including Italian ones. It does not reach five stars because no European animal was sequenced, because the sample size per species is seven-eight individuals, and because the data are not directly transferable to a European breeding decision: they say where the species comes from, not what to do with it tomorrow.

No overall score, by design. A = 3 and C = 5 are not a contradiction: a work with seven animals per species can be the educational reference for an entire field, and this one is.

Bottom line. It can be concluded that the llama descends from the guanaco and the alpaca from the vicuña (subspecies Vicugna vicugna mensalis, once introgressed segments are removed), with 100% bootstrap on 5,901,447 SNPs; that the alpaca genome contains about 36% guanaco/llama ancestry (38.9% with the f_d statistic) against about 5% vicuña/alpaca ancestry in the llama, measured with five independent methods on seven-eight individuals per species; that introgression is significantly lower on the X chromosome than on the autosomes (P < 0.001); and that the exchanged segments contain fleece genes (FGF5, TAGLN), color genes (EDN3, ASIP, CORIN) and a region syntenic to the human blood-pressure locus. It cannot be concluded that the event dates to a calendar date: the estimate is 121.17 generations (95% CI 113.99-130.04) for the alpaca, and its conversion to years depends on a generation length that no one has measured and that the authors deduce, in a declared but circular way, from the date of the conquest itself.

Messaggi da portare a casa

1. Il 36% del genoma di ogni alpaca analizzata è di ascendenza guanaco/lama; nel lama, solo il 5% circa è di ascendenza vigogna/alpaca. Cinque metodi indipendenti, sette-otto individui per specie.
2. Il lama viene dal guanaco, l'alpaca viene dalla vigogna (V. v. mensalis). Ma l'alpaca si colloca correttamente nell'albero solo dopo aver rimosso i segmenti introgressi: l'introgressione deforma la filogenesi.
3. L'evento è datato a 121,17 generazioni fa (IC 95% 113,99-130,04) per l'alpaca e 115,13 per il lama. Generazioni, non anni. In anni: fra 345 e 727, a seconda dell'ipotesi.
4. Sul cromosoma X l'introgressione è molto più bassa (D 0,278 contro 0,560 sugli autosomi, P < 0,001): la selezione ha fatto resistenza dove il costo era più alto.
5. FGF5, il gene del vello lungo, è arrivato all'alpaca dal lama, ed è più espresso nella cute dell'alpaca Suri che dell'Huacaya (P < 0,001). La fibra sta nelle regioni scambiate.
6. Bianco con occhi azzurri: gene candidato EDN3, introgresso dall'alpaca al lama; gli autori ricordano che questi individui sono spesso sordi in entrambe le specie.
7. La dimensione efficace di popolazione di lama e alpaca è in declino da 200 generazioni. La consanguineità è una condizione della specie, non un incidente dell'allevamento.

Take-home messages

1. 36% of the genome of every alpaca analyzed is of guanaco/llama ancestry; in the llama, only about 5% is of vicuña/alpaca ancestry. Five independent methods, seven-eight individuals per species.
2. The llama comes from the guanaco, the alpaca comes from the vicuña (V. v. mensalis). But the alpaca is correctly placed in the tree only after removing the introgressed segments: introgression distorts phylogeny.
3. The event is dated to 121.17 generations ago (95% CI 113.99-130.04) for the alpaca and 115.13 for the llama. Generations, not years. In years: between 345 and 727, depending on the assumption.
4. On the X chromosome introgression is much lower (D 0.278 versus 0.560 on the autosomes, P < 0.001): selection pushed back where the cost was highest.
5. FGF5, the long-fleece gene, arrived in the alpaca from the llama, and is more expressed in Suri alpaca skin than Huacaya (P < 0.001). The fiber lies in the exchanged regions.
6. White with blue eyes: candidate gene EDN3, introgressed from alpaca to llama; the authors note that these individuals are often deaf in both species.
7. The effective population size of llama and alpaca has been declining for 200 generations. Inbreeding is a condition of the species, not an accident of breeding.

Citazione completa

Fan R, Gu Z, Guang X, Marín JC, Varas V, González BA, Wheeler JC, Hu Y, Li E, Sun X, Yang X, Zhang C, Gao W, He J, Munch K, Corbett-Detig R, Barbato M, Pan S, Zhan X, Bruford MW, Dong C. Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca. Genome Biol. 2020;21(1):159. DOI 10.1186/s13059-020-02080-6. PMID 32616020. PMCID PMC7331169. Testo completo libero: sì (CC BY 4.0). Storia della revisione paritaria pubblicata come file supplementare.

Full citation

Fan R, Gu Z, Guang X, Marín JC, Varas V, González BA, Wheeler JC, Hu Y, Li E, Sun X, Yang X, Zhang C, Gao W, He J, Munch K, Corbett-Detig R, Barbato M, Pan S, Zhan X, Bruford MW, Dong C. Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca. Genome Biol. 2020;21(1):159. DOI 10.1186/s13059-020-02080-6. PMID 32616020. PMCID PMC7331169. Free full text: yes (CC BY 4.0). Peer review history published as a supplementary file.

Schede CRL correlate

CRL-E-02 (il gene candidato del Suri secondo il gruppo di Camerino — e i controlli selvatici che vengono dal deposito PRJNA612032 di questo stesso studio) · CRL-E-01 (che cosa misurano, in micron e millimetri, i caratteri di vello che qui compaiono come geni) · CRL-E-05 (che cosa c'era prima: i camelidi dei primi pastori dell'Atacama, letti con il DNA antico che gli autori di questa scheda dichiarano necessario) · CRL-E-03 (il pool genico ristretto visto dal lato dei difetti congeniti e della clinica).

Scheda CRL-E-04 · versione 1.0 · 31/08/2026 · sintetizza una fonte sola, non un corpo di letteratura. Non è stato cercato, e quindi non è escluso, che esistano studi genomici successivi al luglio 2020 che rivedano al rialzo o al ribasso la stima del 36% di introgressione o la datazione in generazioni; non è stata condotta una ricerca sistematica delle citazioni in avanti, che è la prima cosa da fare prima di aggiornare questa scheda. Le tabelle supplementari dell'articolo (S1-S20) non sono state consultate: tutti i valori qui riportati provengono dal testo principale del deposito PMC, e i punti che dipendono dal materiale supplementare sono marcati [DA VERIFICARE SUL TESTO PIENO].

Related CRL factsheets

CRL-E-02 (the candidate gene of the Suri phenotype according to the Camerino group — and the wild controls that come from this same study's PRJNA612032 deposit) · CRL-E-01 (what the fleece traits that appear here as genes measure, in microns and millimeters) · CRL-E-05 (what came before: the camelids of the first Atacama herders, read with the ancient DNA that this factsheet's authors declare necessary) · CRL-E-03 (the restricted gene pool seen from the side of congenital defects and the clinic).

Factsheet CRL-E-04 · version 1.0 · 31/08/2026 · synthesizes a single source, not a body of literature. It has not been searched, and therefore it is not excluded, whether genomic studies subsequent to July 2020 exist that revise upward or downward the 36% introgression estimate or the dating in generations; a systematic search of forward citations has not been conducted, which is the first thing to do before updating this factsheet. The article's supplementary tables (S1-S20) have not been consulted: all values reported here come from the main text of the PMC deposit, and the points that depend on supplementary material are marked [TO BE CHECKED AGAINST THE FULL TEXT].