LECTURA 2. 2021 - Campbell-Staton-et-al.-Science - En.es
LECTURA 2. 2021 - Campbell-Staton-et-al.-Science - En.es
LECTURA 2. 2021 - Campbell-Staton-et-al.-Science - En.es
com
INVESTIGAR
La caza furtiva de marfil y la rápida evolución de la completos a partir de muestras de sangre de 18 elefantes
hembra (norte =7 con colmillos, 11 sin colmillos). Asignamos
falta de colmillos en los elefantes africanos lecturas de secuencia al genoma del elefante de sabana
africana anotado (Loxafr3.0) y generamos alineaciones con una
Shane C. Campbell-Staton1,2,3*†,Brian J. Arnold4,5†,Dominique Gonçalves6,7, Petter Granli8, cobertura de ~ 30 × para 13 muestras y una cobertura de 14 ×
Joyce Poole8, Ryan A. Largo9, Robert M Pringle1 para 5 muestras (materiales complementarios). Utilizando las
muestras de cobertura de 30× (norte =6 con colmillos, 7 sin
Comprender las consecuencias evolutivas de la explotación de la vida silvestre es cada vez más importante a medida que colmillos), calculamos el índice de TajimaD (18)en todo el
la recolección se vuelve más eficiente. Examinamos los impactos de la caza furtiva de marfil durante la Guerra Civil de genoma en ventanas de 10 kb no superpuestas. Ambos grupos
Mozambique (1977 a 1992) en la evolución de los elefantes de la sabana africana (Loxodonta africana)en el Parque mostraron un ligero exceso de variantes raras, indicadas con
Nacional Gorongosa. La caza furtiva resultó en una fuerte selección que favoreció la falta de colmillos en medio de una negativoDvalores (sin colmillos: −0,27, colmillos: −0,2). Sin
rápida disminución de la población. Los datos de la encuesta revelaron patrones de herencia de colmillo consistentes con embargo, las muestras con colmillos tenían significativamente
un rasgo letal masculino dominante ligado al cromosoma X. Las exploraciones del genoma completo implicaron a dos menos variantes raras que las muestras sin colmillos (la
genes candidatos con funciones conocidas en el desarrollo de los dientes de los mamíferos (AMELXyMEP1a), incluyendo la muestra de dos muestras de Welch).tprueba:pag <0.0001) (Fig.
formación de esmalte, dentina, cemento y periodonto. Uno de estos lugares (AMELX)se asocia con un síndrome 1D y materiales complementarios), lo que es consistente con
masculino letal dominante ligado al cromosoma X en humanos que disminuye el crecimiento de los incisivos laterales una contracción más severa de la población de individuos con
superiores (homólogo a los colmillos de elefante). Este estudio proporciona evidencia de una selección rápida mediada colmillos.
por la caza furtiva para la pérdida de un rasgo anatómico prominente en una especie clave.
Para evaluar la respuesta evolutiva a la selección,
cuantificamos la frecuencia de fenotipos de colmillo
D de Tajima
150K
Frecuencia
2500
100K
50K 0.4
2000
(norte=54) (norte=108) (norte=91) 0 − 0,26
Tamaño de la poblacion
18,5% mozambiqueño 50,9% 33% 0,0 0,1 0,2 0,3 0,4 0,5 0,6
0.3
Densidad
1500 Proporción de
Civil con colmillos sin colmillos
C
Guerra 0.2
1000
9,3% 68,5% (1977-1992) 8,3% 40,7% 7,7% 58,4%
registro (probabilidad)
0.1
- 8 -7 -6 -5 -4 -3 -2
500
0.0
0
1970 1980 1990 2000 2010 −2 −1 0 1 2 3
Año 3 4 5 6 7 8 D de Tajima
Sesgo de supervivencia sin colmillos
Fig. 1. Cambios demográficos durante la Guerra Civil de Mozambique y evidencia de selección sin selección para la falta de colmillos. Las líneas verticales sólidas indican las
por falta de colmillos en el Parque Nacional Gorongosa. (A)Cambio en el tamaño de la población proporciones observadas de hembras sin colmillos entre los supervivientes de la
y morfología de los colmillos. La línea representa el tamaño mínimo de la población de elefantes guerra y la primera generación de posguerra (nacidas entre 1995 y 2004). La línea
en el Parque Nacional Gorongosa por año a partir de censos aéreos (30);el segmento discontinua muestra el cuantil superior del 1% de la distribución simulada.
discontinuo indica un período para el cual no hay datos censales sólidos disponibles. Los (C)Logaritmo (probabilidad) del aumento de la supervivencia de las hembras sin colmillos
gráficos circulares muestran la proporción de hembras de dos colmillos (azul oscuro), de un solo observado en Gorongosa, dado un rango de probabilidades de supervivencia de las hembras sin
X−los portadores muestran el fenotipo (es decir, En conjunto, el 87,1% de las asociaciones caza furtiva intensiva. Analizamos los genomas para
penetrancia completa). fenotípicas madre-hijo fueron consistentes con un identificar los loci candidatos más fuertes para la
Encontramos que el 91,3% de las hijas nacidas de modelo de herencia dominante ligado al cromosoma selección en un rasgo letal masculino con sesgo
madres con dos colmillos portaban dos colmillos (norte = X de un solo locus, y el sesgo sexual asociado con la sexual. Primero realizamos un escaneo de todo el
21 con dos colmillos, 1 sin colmillos, 1 con un colmillo) (Fig. falta de colmillos estuvo dentro del 1% de lo esperado genoma en busca de loci que muestren evidencia de
2A). Por el contrario, las hijas de madres sin colmillos bajo letalidad masculina completa. La variabilidad una fuerte selección reciente, específica del fenotipo
mostraron proporciones aproximadamente iguales de inexplicable en el rasgo, incluida la presencia de sin colmillos, y nos enfocamos en genomas de
fenotipos con colmillos y sin colmillos.norte =19 con dos ausencia de colmillos unilateral (si es genética), individuos sin colmillos (norte =11). Aunque la falta de
colmillos (40,9%), 21 sin colmillos (44,7%), prueba de sugiere que las interacciones epistáticas entre al colmillos no se ha fijado en la población de
igualdad de proporciones de dos muestras con corrección menos dos loci pueden influir en la expresión de la Gorongosa (Fig. 1), esperaríamos fuertes señales de
de continuidad,PAG =0,42] (Fig. 2A). Las madres de ambos morfología de los colmillos. Se sabe que las relaciones selección reciente al analizar los genomas del
morfos de colmillo se observaron con hijas que mostraban genotipo-fenotipo asociadas con la variación en la subgrupo sin colmillos.
un fenotipo intermedio de un colmillo [madres de dos morfogénesis dental son muy variables, y se han Usando un análisis de ventana deslizante de
colmillos:norte =1 (4,3%), madres sin colmillos:norte =7 identificado mutaciones distintivas como específicas ventanas de 10 kb con un tamaño de paso de 2 kb,
(14,9%), prueba de igualdad de proporcionesPAG =0,37]. de la población o incluso de la familia.24). Además, el buscamos dos firmas que esperábamos bajo la
Las madres sin colmillos también mostraron una patrón epigenético y la inactivación del cromosoma X hipótesis de selección en un rasgo letal masculino
proporción de sexos de descendencia sesgada (≤5 años de en mosaico pueden dar como resultado la expresión dominante ligado a X. Usando muestras de cobertura
edad,norte =67); El 65,7% eran mujeres, lo que difiere fenotípica variable de mutaciones idénticas asociadas de 30x (6 con colmillos, 7 sin colmillos), primero
significativamente de la hipótesis nula de proporciones de con la agenesia dental, incluso entre gemelos buscamos regiones genómicas que mostraran un
sexos iguales (prueba binomial exacta,PAG = 0.027) y es monocigóticos.25).Sin embargo, los datos de exceso de heterocigosidad dentro de los alelos
estadísticamente indistinguible del 66.7% de sesgo Gorongosa respaldan la hipótesis de que el fenotipo privados de la morfología sin colmillos y
femenino esperado bajo letalidad masculina completa ( sin colmillos está controlado por al menos un locus cuantificamos la heterocigosidad como desviaciones
PAG =0,90) (Fig. 2B). No encontramos evidencia de sesgo letal masculino dominante ligado al X de gran efecto, del equilibrio de Hardy-Weinberg. Luego buscamos
sexual entre las hembras con colmillos (54,2 % de crías con posibles loci modificadores adicionales que loci que mostraran mutaciones altamente
hembras,norte = afectan la expresión fenotípica del rasgo. correlacionadas entre muestras sin colmillos [medido
48,PAG =0.67) (Fig. 2B), y la investigación previa no ha como desequilibrio de ligamiento (LD)]. Cuantificamos
mostrado un nacimiento sesgado por el sexo general Usamos datos del genoma completo para LD mediante el uso de vectores de polimorfismo de
en elefantes africanos (23),lo que sugiere que el sesgo identificar el locus de efecto principal putativo que un solo nucleótido (SNP) (metroLD), tal como se
observado en la proporción de sexos de la sustenta la herencia dominante ligada al X implementa en el RAiSD (26) software.
descendencia se correlaciona con la expresión del hipotética del fenotipo sin colmillos y buscamos Posteriormente, comparamos todos los colmillos (norte
fenotipo sin colmillos. firmas de selección asociadas con recientes =7) y sin colmillos (norte =11) genomas para buscar
µLD
contiguos. Cinco de estas ventanas superponen 305
una región de ~100 kb que contiene la isoforma X: 8
X diente: 5
de amelogenina ligada al X (AMELX) (Fig. 2D), que 13072 4155
codifica una proteína de la matriz extracelular 97 157
MEP1a
A 0.30
0.25
0.20
0.15
primero
0.10
0.05
0
0 20 40 60 80 100 120
B 0.005
0
Δπ
- 0.005
- 0.01
0 20 40 60 80 100 120
C 25
20
15
µLD
10
5
0
0 20 40 60 80 100 120
Posición del cromosoma 1 (andamio 0) (Mb)
Fig. 3. Locus candidato autosómico para la carencia de colmillos en elefantes de equilibrio cuantificado mediante el uso de vectores SNP (metroLD) entre muestras sin colmillos.
Gorongosa. (AaC)Andamio 0 dentro del cromosoma 1 del elefante africano, que destaca Valores atípicos extremos en cada distribución (pag <0.001) se indican con puntos azul oscuro. Los
un locus candidato para la falta de colmillos,MEP1a. (A) Divergencia genética (FS T) entre puntos de color azul claro indican valores no extremos. Las líneas verticales continuas representan
morfos de dos colmillos y sin colmillos. (B) Diversidad relativa (Dp = la posición deMEP1aa lo largo del cromosoma 1 del genoma del elefante. Las líneas discontinuas
asocian con varios trastornos odontogénicos, incluida posición anatómica de (a) esmalte,
la displasia dentinaria, que provoca malformación de d (b) cemento, (c) dentina (marfil), (d)
la raíz del diente y pérdida prematura del diente.34). periodonto y (e) raíz del colmillo.
dentina (35).Estas tres ventanas están anidadas que se vean afectadas por el gen
b
dentro de un intervalo genómico extendido (~6 Mb) candidatoAMELX.Se propone que
que muestra un LD elevado (pag <0.001) (Fig. 3C y fig. los círculos azul claro se vean
a
S2), lo que sugiere una selección positiva reciente en afectados por el gen candidato
distintas del diente de los mamíferos, incluidos el dental (interior negro de la sección
de cualquier estudio que haya establecido vidades sino también para diseñar estrategias 39. E. Abe,Paquidermo22,46–47 (1996).
firmemente una frecuencia de machos sin colmillos para mitigarlas. 40. JS Brashareset al., Ciencia345,376–378
(2014).
más allá de lo que podría explicarse plausiblemente 41. JBS Haldane,J. Genet.44,296–304 (1942).
por lesiones raras o errores del observador (texto 42. KH Andersen, K. Brander,proc. nacional Academia ciencia EE.UU106,
REFERENCIAS Y NOTAS
complementario y tabla S3), pero no podemos 11657–11660 (2009).
1. SK Wassery col., Sci. Adv.4,eaat0625 (2018). 43. M. Heino, B. Díaz Pauli, U. Dieckmann,año Rev. Ecol. Evol.
descartar la posibilidad de mecanismos genéticos 2. S. Wasseret al., Ciencia327,1331–1332 (2010). sist.46,461–480 (2015).
alternativos y/o interacciones genotipo-ambiente . 3. D. Biggs, F. Courchamp, R. Martin, HP Possingham,Ciencias 44. FW Allendorf, JJ Duro,proc. nacional Academia ciencia EE.UU106
Además, los fenotipos intermedios de un solo colmillo 339,1038–1039 (2013). (suplemento 1), 9987–9994 (2009).
4. M. Cardilloet al., Ciencia309,1239–1241 (2005). 45. JS Vertedero,J. Zool.158,293–310 (1969).
comúnmente ocurren en grupos familiares que
5. Onda WJy col., Sci. Adv.1,e1400103 (2015). 46. WH Anciano,Zool. africana158,293–310 (1969).
también incluyen hembras sin colmillos bilaterales.17, 47. AHK Beuchner, HC Dawkins,Ecología42,752–766
6. JJ Midgley, D. Balfour, GI Kerley,S.Afr. J. Ciencia.101,213–215
20, 37). Aunque la evidencia de Gorongosa es (2005). (1961).
consistente con un rasgo letal masculino dominante 7. Médico de cabecera Asnery col., Proc. nacional Academia ciencia EE.UU106,4947– 48. Ala LD, bus IO,Salvaje Monogr.19,3-92
4952 (2009). (1970).
ligado al cromosoma X, los patrones de expresión y 49. RM Pringle,Ecología89,26–33 (2008).
8. MS Waldram, Bonos WJ, Acciones WD,ecosistemas11,
heredabilidad de los colmillos en todo el continente 101-112 (2008). 50. TC Coverdaleet al., Ecología97,3219–3230
pueden ser el resultado de la variación geográfica en (2016).
9. CT Darimonty col., Proc. nacional Academia ciencia EE.UU106,952–954
51. EP Palkovacs, MT Kinnison, C. Correa, CM Dalton,
LD entreAMELX y loci letales masculinos adyacentes, (2009).
AP Hendry,Evol. aplicación5,183–191 (2012).
10. V. Fugère, AP Hendry,proc. nacional Academia ciencia EE.UU115,
loci adicionales en otras partes del genoma de la 52. KB Bienet al., Nat. común11,4155 (2020).
10070–10075 (2018).
especie, variación individual en los patrones de 11. JH Daskin, RM Pringle,Naturaleza553,328–332
EXPRESIONES DE GRATITUD
inactivación del cromosoma X o algún mecanismo (2018).
Agradecemos a A. Gadani, R. Branco, M. Stalmans, P. Tonecas,
genético completamente diferente. Se necesitan más 12. I. Douglas-Hamilton,Orix21,11-24 (1987).
L. Van Wyk, M. Pingo, G. Carr, al personal del Parque Nacional
13. RL Beyerset al., PLOS UNO6,e27129 (2011).
estudios para establecer el número exacto y la Gorongosa y a la República de Mozambique.Fondos:Esta
14. V. Varun, T. Ewing, “La maldición del marfil: la militarización y
identidad de las variantes causales que codifican la investigación fue apoyada por el apoyo financiero a SCC-S. de la
puede causar cambios evolutivos abruptos y persistentes (1995). datos de encuestas sobre crías. JP y PG crearon, mantuvieron y completaron
21. G. Steenkamp, SM Ferreira, MN Bester,JS Af. Veterinario. Asoc. la base de datos de identificación y avistamientos de elefantes. SCC-S. y BJA
en animales longevos, incluso en medio de una
78,75-80 (2007). analizó todos los datos de encuestas, simulación y genómicos. Todos los
disminución extrema de la población. En Gorongosa, la 22. D. Germain, A. Mehta, M. Beck, G. Sunder-Plassmann, en autores participaron en la interpretación y redacción de los datos.Conflicto de
recuperación tanto de la abundancia de elefantes como de Enfermedad de Fabry: Perspectivas de 5 años de FOS,M. Atul, M. intereses:Varios autores tienen relaciones actuales o anteriores con el
Beck, G. Sunder-Plassmann, eds. (Oxford PharmaGenesis, 2006). Proyecto Gorongosa (GP), una organización sin fines de lucro 501(c)(3) de los
la morfología de los colmillos ancestrales puede ser crucial
EE. UU. que apoya la conservación y el alcance comunitario en el Parque
para la restauración del ecosistema. Los colmillos de Nacional Gorongosa, donde se llevó a cabo la investigación. DG es empleado
23. PC Lee, V. Fishlock, CE Webber, CJ Moss,Comportamiento Ecol.
elefante son herramientas multipropósito que se utilizan Sociobiol.70,337–345 (2016). de GP (2017 al presente); PG y JP son directores de ElephantVoices, una
organización sin fines de lucro 501(c)(3) de EE. UU. que recibió una tarifa de
para excavar alimentos y minerales subterráneos (45, 46)y 24. N. Chhabra, M. Goswami, A. Chhabra,Medicina. Patol Oral. Cirugía
consultoría por la investigación realizada en el Parque Nacional Gorongosa;
Oral. bucal19,e112–e119 (2014).
arrancar y pelar la corteza, lo que puede matar árboles RMP fue miembro no remunerado de un comité asesor de GP de 2010 a 2020.
25. M. Varela, MJ Trujillo-Tiebas, P. García-Camba,EUR. Arco.
(47, 48).Estos comportamientos pueden catalizar pediatra Mella.12,318–322 (2011).
Disponibilidad de datos y materiales:Los datos de secuencia están disponibles
a través del archivo de lectura breve dentro de BioProject PRJNA 761769. Los
transiciones de bosque a pastizal a gran escala (45) y crear 26. N. Alachiotis, P. Pavlidis,común Biol.1,79
datos y el código asociados con este estudio están almacenados y disponibles
hábitat para otras especies a escala local (49, 50).En (2018).
públicamente en Github (https://github.com/brian-arnold/
27. KE Holsinger, BS Weir,Nat. Rev. Genet.10,639–650
consecuencia, un aumento en la falta de colmillos en toda AfricanElephantGenomics). La base de datos de quién es quién y su paradero
(2009).
del elefante de Gorongosa está archivada en www.elephantvoices. org/
la población puede tener impactos aguas abajo, como una 28. M. Nei, JC Miller,Genética125,873–879 (1990). gorongosaelephants/. Esta base de datos está protegida con contraseña
bioturbación reducida, cambios en la composición de las 29. HL Viswanatánet al., J. Periodontol.74,1423–1431 como precaución contra la caza furtiva, pero PG y JP pondrán a disposición los
especies de plantas, heterogeneidad espacial reducida y (2003). datos utilizados en este estudio previa solicitud razonable de acuerdo con los
30. CW Gibsony col., J. Biol. química276,31871–31875 acuerdos de confidencialidad. Todos los datos de esta base de datos que son
una mayor cubierta de árboles, cualquiera de los cuales
(2001). necesarios para los análisis realizados en este estudio están disponibles
podría afectar una miríada de otras propiedades del 31. Gerente General Hobsonet al., Am. J.Med. Gineta. UNA.149A,1698-1705 abiertamente en el enlace proporcionado de GitHub.
ecosistema. En otros lugares, la evolución de las especies (2009).
que realizan funciones ecológicas clave ha ejercido efectos 32. X. Duany col., J. Dent. Res.98,870–878 (2019).
33. SK Prakashet al., Hum. mol. Gineta.11,3237-3248
potentes sobre la estructura de la red alimentaria, la
(2002).
composición de la comunidad y el transporte de MATERIALES COMPLEMENTARIOS
34. S. Malik, S. Gupta, V. Wadhwan, GP Suhasini,J. Oral
science.org/doi/10.1126/science.abe7389
nutrientes.51, 52).La restauración de estas funciones Maxilofacial. Patol.19,110 (2015).
Materiales y métodos
puede requerir escalas de tiempo 35. P. Arnoldy col., Cell Tissue Res.367,351–358
Texto complementario
(2017).
desproporcionadamente más largas que el evento de higos. S1 y S2
36. SR, J. Sarukhan, A. Whyte, R. Hassan, “Evaluación de ecosistemas
selección inicial (44)y, por lo tanto, puede limitar el ritmo Tablas S1 a S8
del milenio: ecosistemas y bienestar humano” (Island Press,
Referencias (53–72)
de los esfuerzos de reconstrucción. Comprender la 2006).
Lista de verificación de reproducibilidad de MDAR
37. CJ Moss,Memorias de elefantes: trece años en la vida de una
dinámica de la rápida evolución en el Antropoceno es, por
familia de elefantes (William Morrow, 1988).
lo tanto, esencial, no solo para revelar los impactos 38. IJ Whyte, A. Hall-Martin,Paquidermo55,31–40 11 de septiembre de 2020; aceptado el 3 de septiembre de 2021
La caza furtiva de marfil y la rápida evolución de la falta de colmillos en los elefantes africanos
Materiales y métodos
Texto complementario
higos. S1 y S2
Referencias
Tablas S1 a S8
Lista de verificación de reproducibilidad de MDAR
Materiales y métodos
Encuestas de población de Gorongosa -ElephantVoices (JP y PG) completó diez estudios de campo de ~3-4 semanas en
Gorongosa (2011–2019). ElephantVoices creó y llenó la base de datos de Quién es quién y paradero de Gorongosa
Elephant (http://www.elephantvoices.org/gorongosaelephants/ ), una base de datos en línea que alberga un registro de
búsqueda de elefantes conocidos individualmente ("Quién es quién") y registros geoespaciales de avistamientos de
grupos de elefantes observados, así como aquellos capturados en cámaras de seguimiento ("Paradero"). Los usuarios
pueden, por ejemplo, utilizar la base de datos para identificar un elefante por sus características físicas distintivas o para
buscar el número de individuos sin colmillos en la población por familia, edad o sexo. Esta base de datos está protegida
con contraseña como medida de precaución contra los cazadores furtivos que podrían utilizar los datos geoespaciales.
Identificación y registro de elefantes individuales -Usamos la identificación individual como base para monitorear la
población de elefantes. Los elefantes se caracterizaron y registraron en la base de datos de Gorongosa siguiendo métodos
desarrollados por Amboseli Trust for Elephants (53) y utilizados por ElephantVoices durante más de cuatro décadas. Para
identificar individualmente a los elefantes, usamos sexo, edad, precisión de edad (dentro de 10 años, 5 años, 2 años, 1 año, 6
meses, 1 mes, 1 semana, 1 día), familia, precisión familiar (4 = conocido, 3 =buena idea, 2=suposición fundamentada,
1=desconocida), madre, precisión de la madre (4=conocida, 3=buena idea, 2=suposición fundamentada, 1=desconocida) y
características fisonómicas que incluyen la forma y configuración del colmillo, la forma de la oreja, muescas, rasgaduras y
agujeros, cola y otras marcas corporales.
Utilizamos ocho categorías de edad: 0A (0-4,9 años), 0B (5-9,9 años), 1A (10-14,9 años), 1B (15-
19,9 años), 2 (20-24,5 años), 3 (25-34,5), 4 (35-49,9 años), 5 (50+ años) y se definió como "adulto" a los ≥15
años. Nos enfocamos en la identificación y registro de individuos adultos hasta 2017, cuando comenzamos a
registrar los elefantes inmaduros de familias conocidas. A fines de 2019, habíamos identificado 350 elefantes
adultos, de los cuales 198 eran hembras adultas.
Antes de registrar un nuevo elefante, buscamos cuidadosamente en la base de datos cualquier individuo
coincidente para evitar la duplicación. Asignamos a cada individuo registrado un código de identificación alfanumérico
único que comienza con "gf" para mujeres y "gm" para hombres. Subimos hasta seis fotografías a la tarjeta de
identificación digital, agregando el nombre del fotógrafo y el año tomado como metadatos.
Localización y registro de elefantes -Tenemos registros fotográficos de madres conocidas y sus crías que se
remontan a 2011. Recorrimos las pistas principales para buscar grupos de elefantes y usamos la aplicación de
teléfono celular GPS-Trk para llevar un registro de nuestra ruta. Gorongosa está muy arbolado y sus elefantes son
muy defensivos. Los grupos familiares en particular suelen reaccionar ante la aproximación de vehículos con
vigilancia, evasión, confrontación o acoso y ataque. Por lo tanto, cuando avistamos un elefante o un grupo de
elefantes, intentamos inmediatamente fotografiar tantos individuos como fuera posible. Usamos una Canon 6D con
una lente 200-400 y un extensor 1.4 para tomar fotografías de identificación desde una distancia considerable para
no molestar a los elefantes.
Los registros geoespaciales de grupos de elefantes ("avistamientos") se recopilaron mediante Gorongosa EleApp,
una aplicación de Android personalizada diseñada por ElephantVoices que registra directamente la fecha, la hora y la
ubicación del GPS. Ingresamos manualmente el observador, el nombre del lugar (si se conoce), el tipo de grupo (grupo de
machos, grupo familiar, grupo familiar con machos asociados), número de individuos, precisión del conteo (exacta, buena
estimación, conjetura), apellido si se conocen, y los machos y hembras individuales reconocidos. Una vez conectado a WiFi,
los avistamientos registrados en EleApp se cargaron en la base de datos en línea. Las notas de campo generalmente se
registraban en un teléfono móvil y luego se transcribían en el registro en línea. Las notas de campo incluyeron la reacción
de los elefantes al vehículo y otros comportamientos, y,
cuando sea posible, un censo del grupo, incluidas las hembras adultas presentes y la edad y sexo de sus crías. Los
elefantes que no pudimos reconocer en el campo se identificaron más tarde usando la base de datos y sus identificaciones
se agregaron al registro de avistamientos. Se registraron nuevos individuos como se describe anteriormente. Usamos
Adobe Bridge para escribir palabras clave en cada fotografía con el número de registro de todos los elefantes
fotografiados.
Proporciones históricas del fenotipo del colmillo:Para estimar la frecuencia de falta de colmillos en Gorongosa antes de
la guerra civil, JP clasificó el estado de los colmillos de elefantes individuales a partir de catorce videos históricos grabados
dentro del parque antes del conflicto (antes de 1972; Table_S1.xlsx). En todas las imágenes de video, se identificaron 149
elefantes adultos individuales (89 machos, 54 hembras). Las personas cuya categoría de colmillo no se pudo determinar
de manera confiable (debido a la posición del cuerpo o al tiempo en la pantalla) se clasificaron como "sin clasificar". Para
estimar la frecuencia de ausencia de colmillos después de la guerra civil, utilizamos individuos registrados en la base de
datos de elefantes de Gorongosa (Tabla_S2.xlsx). Esta base de datos incluía 217 mujeres registradas y 174 hombres
registrados.
Observaciones contemporáneas de madre e hijo -Para los fines de este estudio, buscamos en nuestras notas de campo, 5
TB de imágenes de video y 30 000 fotografías tomadas durante las observaciones de 280 grupos de elefantes para identificar
crías cuyas madres se conocían con un alto grado de precisión. Ingresamos a cada elefante inmaduro al que pudimos
asignar un puntaje de precisión de la madre de ≥2 en una hoja de cálculo. La hoja de cálculo incluía la siguiente información:
nombre de la madre y código de identificación, configuración de los colmillos de la madre (2 colmillos, 1 colmillo, sin
colmillos), precisión de la madre (utilizando las categorías descritas anteriormente), notas de campo, fecha de registro de la
cría por primera vez, edad estimada de la becerro en el primer registro, código de identificación del becerro, fecha de
nacimiento estimada del becerro, puntaje de precisión de la fecha de nacimiento, sexo del becerro y configuración de
colmillos del becerro (desconocido, 2 colmillos, 1 colmillo, sin colmillos). Los elefantes inmaduros de hasta 8 años de edad se
encuentran a menos de 5 metros de sus madres el 80% del tiempo,54). A medida que los terneros envejecen, tienden a
encontrarse a mayor distancia de sus madres, particularmente en el caso de los machos jóvenes. Debido a esta
diferenciación entre los sexos, limitamos nuestro análisis a los terneros que fueron< 5 años de edad en el primer registro.
Los terneros que se registraron por primera vez como bebés (se estima que son< 1 año de edad) y fueron vistos en estrecha
asociación con una hembra adulta, junto con crías o crías que fueron observadas, filmadas o fotografiadas amamantando,
se les asignó una puntuación de precisión de madre de 4. Las crías y crías observadas en estrecha asociación con,
orientando hacia, y/o siguiendo persistentemente a una hembra adulta en particular, se les asignó un puntaje de precisión
de la madre de 3. Para el análisis, incluimos solo terneros con puntajes de precisión de la madre de 3 o 4. De los 270 terneros
ingresados en la hoja de cálculo, 130 eran de sexo conocido, tenían una precisión de la madre de 3 o 4, y eran< 5 años de
edad en el momento del primer registro. Los análisis estadísticos en el texto principal excluyen individuos con colmillos
únicos (se informan números y porcentajes) porque los individuos con colmillos únicos pueden ocurrir debido a la rotura de
colmillos y, por lo tanto, representan una categoría ambigua para comprender la heredabilidad de la ausencia de colmillos
con base genética. Estos datos de avistamientos madre-hijo se incluyen en Table_S4.xlsx madre-hijo, y todo el código R
utilizado para filtrar y analizar estos datos está archivado en GitHub (y estará disponible públicamente tras su publicación).
Selección de hembras sin colmillos en Gorongosa-Con estimaciones del tamaño de la población antes y después de la guerra, y
asumiendo de manera conservadora una proporción de sexos igual de machos y hembras, utilizamos la distribución de Wallenius
(paquete BiasedUrn R) para modelar las probabilidades de que las hembras sin colmillos sobrevivan a la guerra. Usando un valor de
probabilidad de 1 correspondiente a ninguna selección, primero simulamos la proporción esperada de hembras sin colmillos
después de la guerra (Figura 2a). Luego calculamos la probabilidad, o la probabilidad
de observar nuestros datos dado un valor de probabilidades, para probabilidades que van del 1 al 10 en intervalos de 0.001. Un
valor de probabilidad de 10 significa que las hembras sin colmillos tienen 10 veces más probabilidades de sobrevivir a la guerra que
las hembras con colmillos. Esta curva de probabilidad aproximada (Figura 2b) se usó para inferir la estimación de máxima
verosimilitud del valor de probabilidad, junto con el intervalo de confianza del 95 % de esta estimación utilizando una caída en 1,92
unidades logarítmicas de verosimilitud a partir de la estimación de máxima verosimilitud. Todo el código R utilizado para estos
cálculos está archivado en GitHub (y estará disponible públicamente tras su publicación).
Captura de animales y muestreo genético -Durante mayo-agosto de 2018, inmovilizamos químicamente a 18 elefantes
hembra (norte=7 colmillos,norte=11 sin colmillos) por inyección remota (dardos) desde un helicóptero utilizando una
combinación de oxalato de tiafentanilo y azaperona. La dosificación se basó en el tamaño aproximado y la edad del
individuo, que evaluamos visualmente en el momento de la captura. Supervisamos cuidadosamente a cada elefante durante
el manejo y medimos la frecuencia cardíaca (normal: 25–30 lpm), la frecuencia respiratoria (normal: 4–6 respiraciones/
minuto) y la temperatura rectal (normal: 36–37oC) a intervalos de 5 min. También cubrimos los ojos de cada individuo para
minimizar la exposición a estímulos externos. Mientras estaba inmovilizado, a cada elefante se le colocó un collar GPS de
iridio (Savannah Tracking, Nairobi, Kenia) para rastrear los movimientos como parte de un estudio relacionado. Obtuvimos
una muestra de sangre (10 mL) de cada individuo por venopunción de una de las principales venas de la oreja. Cuando se
completó la recolección de muestras, antagonizamos el thiafentanil con naltrexona y observamos a cada elefante desde una
distancia segura hasta que recuperó el equilibrio y se alejó del área. Todos los procedimientos de manipulación de animales
fueron aprobados por el Comité de uso y cuidado de animales de la Universidad de Idaho (protocolo n.° 2015-39) y se
realizaron de acuerdo con las pautas establecidas por la Sociedad Estadounidense de Mamólogos (55).
Alineación de secuencias y llamada de variantes-Las lecturas se procesaron mediante un flujo de trabajo de Snakemake (
https://github.com/harvardinformatics/shortRead_mapping_variantCalling ) para crear un archivo VCF que contenga SNP de
alta calidad para análisis posteriores. Específicamente, las lecturas se filtraron primero con fastp (56) utilizando la opción
"detect_adapter_for_pe" para eliminar las lecturas de baja calidad y recortar las secuencias de adaptador restantes.
Posteriormente, las lecturas se asignaron al ensamblaje del genoma Loxafr3.0 para Loxodonta africanausando el algoritmo
BWA-MEM (57) (junto con la opción -M); Los duplicados de PCR se marcaron para análisis posteriores utilizando Picard (
http://broadinstitute.github.io/picard ). En Table_S5.xlsx se proporcionan varios resúmenes del proceso de mapeo y filtrado
de lectura, incluidas las fracciones de lecturas de alta calidad, alineaciones de lectura de alta calidad y duplicados de PCR,
junto con las profundidades de secuenciación medias de todo el genoma. Secuenciamos 13 muestras a mayores
profundidades (promedio ~30X) y las 5 muestras restantes a profundidades intermedias (promedio ~10-15X) debido a los
bajos rendimientos de ADN.
Luego usamos lecturas alineadas para detectar SNP e indels pequeños usando el flujo de trabajo
GATK (v. 4.1.8) HaplotypeCaller. Los sitios SNP de baja calidad se filtraron utilizando la expresión "QD < 2,0 ||
FS > 60,0 || SOR > 3,0 || MQ < 40,0 || MQRankSum < -12,5 || ReadPosRankSum < -8,0 ", y no filtramos sitios
con exceso de heterocigosidad para evitar el descarte de sitios potencialmente relevantes para este estudio.
Después del filtrado, teníamos 11 931 133 SNP. Por último, anotamos variantes en el VCF final usando snpEff (
58) (base de datos Loxafr3.99). Aunque asignamos lecturas a la referencia Loxafr3.0,
utilizó información del ensamblaje LoxAfr4 para agrupar andamios en cromosomas para análisis
particulares.
Control de calidad -Utilizamos vcftools (59) para calcular la relación (estadística Ajk no ajustada) de los
individuos entre sí para verificar que las muestras no se originaron en la misma familia
(Tabla_S6.xlsx). Dos individuos, 0045B y 2986A, parecían ser hermanos (Ajk » 0.5), por lo que solo usamos el
individuo 0045B para análisis posteriores para evitar que la relación sesgue los resultados.
Además del filtro SNP por sitio aplicado anteriormente, aplicamos filtrado adicional por
individual para análisis que dependían de genotipos de alta calidad: Tajima's D, LD (usando vectores SNP en
RAiSD), heterocigosidad y diversidad relativa (pagssin colmillos-pagsdos colmillos), descrito abajo. Para estos análisis, solo usamos las
13 muestras (7 sin colmillos, 6 con colmillos) que tenían profundidades altas (~30X; Table_S5.xlsx) y solo consideramos sitios
variables en los que cada uno de estos individuos tenía al menos 8X de profundidad pero no más de 80X profundidad. No
consideramos ningún sitio al que le faltaran datos. En particular, nuestras restricciones de profundidad de secuenciación
por individuo y la intolerancia de los datos faltantes no descartaron una cantidad significativa de sitios; de los ~13,5 millones
de SNP que cumplieron con un umbral de profundidad mínima inferior de 4X por persona, el 97 % se retuvo cuando se
utilizó nuestra profundidad mínima superior de 8X por persona (Tabla_S7.xlsx). Los únicos análisis que utilizaron todas las
muestras, de profundidad alta e intermedia, involucraron ANGSD (60) para evaluar la diferenciación entre muestras sin
colmillos y con colmillos (FS TyDXY). Este programa utiliza las probabilidades de genotipo para dar cuenta de la incertidumbre
dentro de las regiones de menor profundidad de secuenciación. Para obtener una lista de qué muestras se usaron para qué
análisis, consulte Table_S8.xlsx.
Para definir ventanas genómicas para análisis, buscamos bloques de SNP en LD relativamente alto. Utilizamos vcftools (59
) para cuantificar la DL de genotipo por sitio (medida comor2) para nuestras muestras de alta profundidad y descubrió que los
bloques de 10 kb normalmente teníanr2valores > 0.5 (Figura S1). Dado que las ventanas genómicas para muchos de nuestros
análisis se definieron en función de los SNP, elegimos tamaños de ventana SNP tales que la longitud media de la ventana, en pares
de bases, fuera de aproximadamente 10 kb o menos (Table_S8.xlsx).
Análisis genómico de las diferencias demográficas entre fenotipos de colmillo–Utilizamos vcftools (59) para calcular el TajimaD
en ventanas de 10 kb que no se superponen para grupos con colmillos y sin colmillos. Utilizando todas las ventanas de 10 kb para las
que había estimaciones disponibles, comparamos la diferencia media enDentre grupos usando undos muestras de Welcht-prueba.
Para asegurar que las diferencias observadas en la mediaDno se vieron afectados significativamente por las diferencias en el tamaño
de la muestra entre los grupos (7 sin colmillos, 6 con colmillos) o la elección de individuos específicos, volvimos a ejecutar el análisis
siete veces más, reduciendo iterativamente el grupo sin colmillos a seis individuos mediante la eliminación de un individuo diferente.
Los resultados de los análisis de muestreo reducido siguieron siendo estadísticamente significativos en cada caso. Por lo tanto,
informamos la comparación original entre los grupos de morfología de colmillo utilizando todas las muestras de alta profundidad.
Escaneos de selección y desequilibrio de ligamiento -Para detectar firmas de selección dentro de elefantes sin
colmillos, usamos RAiSD (72), que calcula varias métricas sensibles a los barridos. Calculamos estas métricas en
ventanas deslizantes de 50 SNP a lo largo del genoma (con valores predeterminados para otros parámetros) y solo
usamos la métrica que cuantifica LD usando vectores SNP (µviejo). Las puntuaciones individuales por ventana se
construyen de manera que los valores altamente positivos sean consistentes con la selección reciente (es decir, LD
alto).
Identificación de loci diferenciados -Utilizamos un análisis de ventana deslizante a lo largo del genoma para cuantificar la
diferenciación genética entre muestras con colmillos y sin colmillos. Este análisis utilizó dos estadísticas genéticas de población:FS Ty
DXY.FS Tutiliza diferencias normalizadas en las frecuencias alélicas entre grupos de muestras. TiempoFS Tes una forma popular de
cuantificar la diferenciación, sus valores dependen en gran medida de los niveles de diversidad, de modo que se espera que las
regiones genómicas con baja diversidad dentro de los grupos tengan valores más altos (60, 61). Por lo tanto, para asegurarnos de
que detectamos regiones genómicas que eran verdaderamente
diferenciado entre muestras sin colmillos y con colmillos, y no solo regiones de baja diversidad, también usamosDXY
(62), que es independiente de la diversidad de secuencias dentro de los grupos. Esta estadística se calcula mediante
el número promedio de diferencias por pares entre secuencias de dos grupos, ignorando todas las comparaciones
entre secuencias dentro de grupos.
Para calcular estas estadísticas, utilizamos ANGSD (60), que utiliza las probabilidades de genotipo para dar cuenta
de la incertidumbre dentro de las regiones de menor profundidad de secuenciación. Por lo tanto, utilizamos las 18 muestras
para estos análisis, ya que ANGSD tiene en cuenta la incertidumbre del genotipo para las muestras de profundidad
intermedia (Tabla_S5.xlsx). Tuvimos ANGSD calcularFS Ten ventanas de 10kb a lo largo del genoma, y también lo usamos
para calcular frecuencias de alelos menores dentro de muestras con colmillos y sin colmillos, excluyendo sitios que tenían
una frecuencia de alelos menores (en las 18 muestras) < 10%, y también excluyendo bases de baja calidad (BQ < 20) y
lecturas mapeadas de baja calidad (MQ < 20). Usando ventanas deslizantes de 50 SNP (con un tamaño de paso de 10 SNP),
calculamosDXYcomo∑($1 (́1 − $2) + $2 (́1 − $1))dóndepags1 ypags2 representan frecuencias alélicas en grupos sin colmillos y
con colmillos, respectivamente. Valores deDXYfueron normalizados por el número de SNP en cada ventana (63,64). Estas
ventanas 50-SNP usadas para calcularDXYtenía una longitud similar, en pares de bases, a otros análisis (Tabla_S8.xlsx).
Análisis de heterocigosidad -Para estudiar los niveles de heterocigosidad a lo largo del genoma, usamos cyvcf2 (26) para analizar
nuestro archivo VCF y calcular las desviaciones locales de la heterocigosidad observada a partir de las proporciones.
Para estos análisis de heterocigosidad, solo consideramos los SNP que eran comunes en las muestras sin colmillos (15-85
% de frecuencia) pero que estaban ausentes en las muestras con colmillos. Si bien esperamos que una mutación letal masculina
ligada al cromosoma X sea heterocigota en las mujeres, es posible que los SNP vinculados a esta mutación no lo sean, especialmente
si estos SNP vinculados tenían frecuencias intermedias cuando surgió la mutación causal por primera vez y estaban presentes tanto
en mutantes como en los de tipo salvaje. antecedentes genéticos. Sin embargo, los SNP vinculados que tenían frecuencias muy
bajas, como los exclusivos del trasfondo genético en el que surgió la mutación causal, pueden enriquecerse para los genotipos
heterocigotos después de que el alelo mutante letal masculino ligado al X experimentara una selección positiva y aumentara la
frecuencia de estos SNP vinculados. Por lo tanto, para apuntar a este subconjunto de SNP vinculados que también pueden
enriquecerse para genotipos heterocigotos, solo consideramos los SNP que eran comunes en muestras sin colmillos (15-85% de
frecuencia) pero ausentes en muestras con colmillos. Dado que este procedimiento de creación de subconjuntos aumenta la
distancia entre SNP de los SNP que se analizan en última instancia, utilizamos ventanas de 10 SNP (con un tamaño de paso de 2 SNP),
que correspondían a ventanas genómicas de aproximadamente 10 kb de tamaño. Este tamaño de ventana SNP más pequeño se
eligió para hacer que el tamaño de estas ventanas (en pares de bases) sea más comparable con los otros análisis (Table_S8.xlsx).
Análisis de la diversidad genética -La diversidad genética se calculó para individuos con colmillos y sin colmillos por
separado en ventanas de 10 kb sin superposición utilizando vcftools. La diversidad relativa se calculó como la
diferencia en la diversidad observada entre el grupo sin colmillos y el grupo con colmillos (pagssin colmillos-pagscon colmillos). Si el
locus candidato respalda el fenotipo sin colmillos, esperaríamos una diversidad genética significativamente menor, en
comparación con el promedio de todo el genoma, en la región candidata debido a la fuerte selección reciente dentro de los
individuos sin colmillos. Alternativamente, la diversidad dentro de los individuos con colmillos no debe ser
significativamente más bajo que el promedio de todo el genoma. Significarpagsel genoma completo fue casi idéntico
para ambos grupos, con una media de 0,001666 para muestras con colmillos (desv. constante = 0,00295) y 0,001668 para
muestras sin colmillos (desv. constante = 0,00287). Dentro de las 3 ventanas candidatas asociadas con MEP1a (scaffold_0:
posición 117320001 – 117330000, posición 117330001 – 117340000 y posición
117340001 – 117350000),pagslas estimaciones fueron 0,009, 0,009 y 0,011, respectivamente, para el grupo con colmillos, y 0,0002,
0,0004 y 0,0004, respectivamente, para el grupo sin colmillos. Por lo tanto, mientras observamos una disminución de la diversidad
dentro de las muestras sin colmillos, las muestras con colmillos también tenían una diversidad elevada dentro de las muestras.
región candidata. Esto es impulsado por un conjunto de 246 variantes que no son de referencia que están todas fijadas dentro
de muestras sin colmillos pero presentes en frecuencias intermedias en muestras con colmillos (ver “Región candidata en el
cromosoma 1" abajo).
Análisis de SNP diferenciados-Además de nuestro análisis basado en ventanas para detectar regiones diferenciadas
entre muestras sin colmillos y con colmillos (usandoFS TyDXY,arriba), usamos cyvcf2 para escanear nuestro VCF en
busca de SNP individuales altamente diferenciados e indeles pequeños con diferencias de frecuencia alélica entre
muestras sin colmillos y con colmillos que eran ≥ 50%. Usamos solo aquellos individuos que tenían una alta
profundidad de secuenciación. Cuantificamos la diferenciación de SNP utilizando la diferencia total de frecuencia de
alelos en todos los alelos que no son de referencia. Por ejemplo, si las muestras sin colmillos tuvieran múltiples alelos
que no son de referencia que no se encuentran en las muestras con colmillos, cada uno contribuiría a la diferencia
general de frecuencias alélicas. Si bien la mayoría de los SNP son bialélicos en nuestras muestras (95,6 %), utilizamos
este procedimiento para acomodar los sitios multialélicos menos frecuentes que pueden ser biológicamente
relevantes. Curiosamente, a pesar de las altas profundidades de secuenciación y relativamente pocas muestras por
grupo,
Análisis de grandes variantes estructurales-Para identificar variantes estructurales más grandes no detectables por GATK4,
utilizamossv-callers(66),un flujo de trabajo de Snakemake que ejecuta cuatro llamadores de variantes estructurales: Manta (v1.1.0),
DELLY (v0.7.7), LUMPY (v0.2.13) y GRIDSS (v1.3.4). El flujo de trabajo de sv-callers utiliza los llamantes de variantes estructurales en
cada muestra individual, por lo que combinamos llamadas entre muestras usando SURVIVOR (67) y llamadas combinadas que eran
del mismo tipo (p. ej., eliminación, inserción) y tenían puntos de interrupción con una diferencia de 1 kb entre sí. Dentro de estos
archivos VCF fusionados, usamos el campo "SUPP_VEC", creado por SURVIVOR, en la columna INFO. Este campo indica la presencia o
ausencia de una variante estructural dentro de cada muestra, pero no proporciona información sobre el genotipo. Para centrarnos
en las variantes que podrían explicar la diferencia en el estado de los colmillos entre los elefantes sin colmillos y con colmillos, solo
consideramos las variantes que tenían frecuencias apreciablemente más altas dentro de las muestras sin colmillos, es decir,
frecuencias que eran al menos un 50 % más altas en las muestras sin colmillos en comparación con las muestras con colmillos.
Elementos reguladores cis candidatos (cCRE)–Para evaluar si los SNP no exónicos ocurrieron dentro de los
elementos reguladores potenciales, buscamos cCRE en la base de datos ENCODE dentro del genoma humano
(compilación hg19). Para cada gen dentro de las dos ventanas candidatas mencionadas anteriormente (una en el
cromosoma X, otra en el cromosoma 1), reunimos todos los cCRE que estaban (1) dentro del cuerpo del gen, (2)
entre el primer y el último sitio de inicio de la transcripción, o (3) hasta 50 kb aguas arriba del sitio de inicio de la
transcripción. Luego tomamos las coordenadas de estos cCRE y usamos la herramienta de elevación del navegador
del genoma UCSC para traducir las coordenadas del genoma hg19 al genoma del elefante (loxAfr3.0). Si bien este
análisis brinda información adicional sobre el impacto potencial de las variantes no exónicas, no sabemos si estos
cCRE son biológicamente significativos en los elefantes africanos sin trabajo experimental adicional.
Texto complementario
Región candidata en el cromosoma X -Dentro de la región candidata scaffold_39: 19913744 - 20932476
(Figura 2d), hay un solo sitio con una frecuencia de alelos del 50% en muestras sin colmillos pero
completamente ausente en muestras con colmillos. Este sitio contiene una indel multialélica intergénica en la
que el alelo de referencia es T y los alelos alternativos son TG, TGG y TGGG. Cuatro muestras sin colmillos eran
heterocigotas para el alelo TGG, una muestra era heterocigota para el alelo TG, una muestra era heterocigota
para los alelos TG y TGG (es decir, sin alelo de referencia) y la última muestra sin colmillos era homocigota
para el alelo de referencia. En comparación, las muestras de dos individuos con dos colmillos eran
heterocigóticas para el alelo TGGG, que no se encontró en ninguna de las muestras sin colmillos. Más lejos
hay tres SNP aguas arriba del gen HCCS, y estos implican la no referencia bialélica.
alelos exclusivos de las muestras sin colmillos. Uno de estos SNP se encuentra dentro de un elemento regulador cis
candidato (cCRE, descrito anteriormente). En estos tres sitios, cinco muestras sin colmillos eran heterocigotas para el alelo
alternativo, una era homocigota y una tenía el alelo de referencia en los tres sitios.En términos de variantes estructurales
grandes, DELLY detectó una eliminación de 2,6 kb dentro del intervalo, pero estuvo presente en 4 de 7 muestras sin colmillos
y 1 de 6 muestras con colmillos. Fuera del intervalo candidato, tanto DELLY como LUMPY detectaron una eliminación de 1 kb
133 kb aguas arriba del intervalo, presente en 3 de 7 muestras sin colmillos y ausente en las muestras con colmillos.
Región candidata en el cromosoma 1 -Dentro de la región candidata scaffold_0:117320001-117350000 (Figura 3), hay un
conjunto de 246 variantes (SNP e indeles pequeños) que se fijaron todas en muestras sin colmillos pero presentes con una
frecuencia del 50 % en muestras con colmillos. Estas variantes son en su mayoría intergénicas entre MEP1a y PLA2G7 (norte
=126) o dentro de los intrones MEP1a (norte=86). Si bien ninguna de estas variantes se encuentra dentro del promotor de
MEP1a,norte=89 están dentro de cCRE (descritos a continuación), incluidos 7 pequeños indels. Además de estas variantes no
codificantes, existen 3 mutaciones no sinónimas en los exones de MEP1a. Dentro de toda esta ventana, cuatro de las seis
muestras con colmillos eran heterocigotas en todas estas posiciones variantes, mientras que una era homocigota para
todos los alelos alternativos y otra era homocigota para todos los alelos de referencia (totalizando una frecuencia alélica del
50 % en las muestras con colmillos en todas estas posiciones). variantes). Las muestras sin colmillos eran homocigóticas
para los alelos alternativos en todas las variantes de esta región.En términos de grandes variantes estructurales, no existía
nada dentro de la región candidata. Sin embargo, Manta detectó una inserción de 104 pb 411 kb aguas arriba de la región
que estaba presente en 5 de 7 muestras sin colmillos pero ausente en las muestras con colmillos. Además, DELLY detectó
una eliminación de 28 pb 364 kb aguas abajo de la región que estaba presente en 4 de 7 muestras sin colmillos pero
ausente en las muestras con colmillos.
Un mecanismo candidato para la ausencia de colmillos de base genética en los elefantes africanos -En los
mamíferos, las hembras portan el doble de cromosomas X que los machos, pero logran la misma dosis de genes a
través de la inactivación del cromosoma X, en la que una copia del cromosoma X se silencia al principio del desarrollo.
Las células somáticas pueden variar en qué copia silencian, lo que da como resultado un mosaicismo fenotípico
cuando alelos diferentes en cromosomas X alternativos dan lugar a fenotipos diferentes. En condiciones maleletales
dominantes ligadas a X, incluidas las asociadas con AMELX, la inactivación de X se ha identificado o se ha planteado la
hipótesis como un factor que determina los resultados fenotípicos de la variación genética.68–71). La variabilidad
fenotípica asociada con la inactivación del cromosoma X (XCI) es un fuerte candidato como un mecanismo que
impulsa la variabilidad en la herencia y presentación de la morfología del colmillo en Gorongosa. La región candidata
ligada al sexo es sinténica con la región Xp22.2 del cromosoma X humano, que está asociada con el síndrome de la
enfermedad MIDAS/XAI (71). La herencia letal masculina dominante ligada al cromosoma X de esta enfermedad se
asocia comúnmente con aneuploidía segmentaria en toda la región. Sin embargo, la variación en la expresión de
rasgos es muy variable dentro y entre familias. En ciertos casos, las madres asintomáticas pueden tener hijas que
presenten el rasgo (70). En tales casos, madre e hija comparten el genotipo causal, pero el individuo asintomático
muestra una inactivación preferencial del cromosoma afectado.70). Dada la homología entre la región
humanaXp.22.2 con la región candidata ligada al cromosoma X del genoma del elefante, proponemos que la
inactivación sesgada de X puede explicar instancias aberrantes en Gorongosa en las que una madre con dos colmillos
se asoció con una hija sin colmillos, como se observó en una sola caso en este estudio.
Encuesta literaria–Para evaluar la frecuencia de ausencia de colmillos en las poblaciones de elefantes africanos,
revisamos sistemáticamente la literatura (Tabla_S3.xlsx). El 20 de abril de 2020, buscamos en Google Scholar usando
la frase ["elefante africano" sin colmillos] e ISI Web of Science usando la frase [(sin colmillos* O "colmillo perdido*"
O "falta* de colmillo*" O "pérdida de colmillo*" ) Y elefante*]. Estas búsquedas arrojaron 308 y 18 resultados,
respectivamente, todos los cuales escaneamos; aquellos que parecían probables presentar
Se recuperó y revisó información cuantitativa sobre la carencia de colmillos en los elefantes africanos. Exploramos más la literatura citada en estas referencias para identificar fuentes potenciales
adicionales, y recuperamos datos de una base de datos en línea administrada por JP y PG. No se pudo obtener una pequeña cantidad de fuentes potenciales. En total, revisamos 87 libros, artículos y
disertaciones; Table_S3.xlsx incluye información de 24 fuentes, lo que produjo 58 registros en total. Los estudios difirieron en lo que informaron y cómo. Siempre que fue posible, extrajimos datos
sobre las frecuencias generales y específicas del sexo sin colmillos; los tamaños de muestra de los que se derivaron estas frecuencias; si se informaron individuos con un solo colmillo; la presión
cinegética en el lugar (evaluada cualitativamente a partir de la información de las fuentes); si los autores informaron indicios de que la falta de colmillos es hereditaria; la hipótesis de los autores sobre
el mecanismo de ausencia de colmillos; el formato en que se presentaron los datos; y cualquier nota pertinente sobre y/o citas de la fuente principal. Las celdas en blanco indican que la información
no se presentó y no se pudo inferir de la fuente. Los signos de interrogación indican incertidumbre (ya sea por calificaciones explícitas presentes en la fuente original o porque no estábamos seguros
de nuestra interpretación de la información en la fuente; en estos casos, explicamos más en las notas). Solo dos fuentes informaron frecuencias apreciables (> 2%) de falta de colmillos en elefantes
africanos machos ( y cualquier nota pertinente sobre y/o citas de la fuente principal. Las celdas en blanco indican que la información no se presentó y no se pudo inferir de la fuente. Los signos de
interrogación indican incertidumbre (ya sea por calificaciones explícitas presentes en la fuente original o porque no estábamos seguros de nuestra interpretación de la información en la fuente; en
estos casos, explicamos más en las notas). Solo dos fuentes informaron frecuencias apreciables (> 2%) de falta de colmillos en elefantes africanos machos ( y cualquier nota pertinente sobre y/o citas
de la fuente principal. Las celdas en blanco indican que la información no se presentó y no se pudo inferir de la fuente. Los signos de interrogación indican incertidumbre (ya sea por calificaciones
explícitas presentes en la fuente original o porque no estábamos seguros de nuestra interpretación de la información en la fuente; en estos casos, explicamos más en las notas). Solo dos fuentes
informaron frecuencias apreciables (> 2%) de falta de colmillos en elefantes africanos machos ( en estos casos, explicamos más adelante en las notas). Solo dos fuentes informaron frecuencias
apreciables (> 2%) de falta de colmillos en elefantes africanos machos ( en estos casos, explicamos más adelante en las notas). Solo dos fuentes informaron frecuencias apreciables (> 2%) de falta de
colmillos en elefantes africanos machos (72,73), y consideramos que estas dos estimaciones son inciertas según las declaraciones de las fuentes originales.
Fenotipo intermedio de un solo colmillo-Los individuos con un solo colmillo representan una categoría ambigua para
nuestros análisis porque las lesiones pueden resultar en la pérdida de un colmillo, mientras que la pérdida adquirida y
simétrica de ambos colmillos es extremadamente rara; por lo tanto, se puede suponer con seguridad que la falta de
colmillos bilateral es congénita (74, 75), pero los casos de falta de colmillos unilaterales son menos seguros.Varios estudios
incluyen informes anecdóticos de este tipo de lesiones, que pueden conducir a la infección y la pérdida del diente (76–79).
Hasta donde sabemos, los datos cuantitativos más extensos son de Steenkamp et al. 2007 (75), quien evaluó fotografías de
~1500 animales de 15 poblaciones en 7 países en busca de evidencia de fractura de ambos colmillos (norte=2695 colmillos
examinados) y ausencia de colmillos congénita (ver Table_S3.xlsx). La proporción de colmillos fracturados fue muy variable
entre las poblaciones (rango 0–44 %), con un promedio de 9,4 % entre las hembras y 7,4 % en general (la rotura bilateral fue
muy rara, lo que sugiere que hasta ~19 % de las hembras y ~15 % de todos los individuos). tenía un colmillo roto). Sin
embargo, la proporción de colmillos fracturados disminuyó exponencialmente con la lluvia; considerando solo las
poblaciones con precipitaciones ≥ 800 mm año-1(comparable a Gorongosa, 840 mm año-1), solo 6 de 301 colmillos estaban
fracturados (2%, lo que sugiere ~4% de los individuos). Por lo tanto, la frecuencia de hembras con un solo colmillo en
Gorongosa (~7,5 %) cae dentro del rango de frecuencias de fractura de colmillo en toda África, pero es más alta de lo
esperado dadas las lluvias en Gorongosa. De hecho, los registros en la base de datos de elefantes de Gorongosa indican que
la frecuencia de rotura de colmillos unilaterales entre sexos es del 3,1% (norte=260), predominantemente en animales
mayores, lo que sugiere que la lesión por sí sola no puede explicar el fenotipo de un solo colmillo en nuestro estudio.
La evidencia de varios estudios también indica que los colmillos únicos tienden a darse en familias, a
menudo aquellas con altas proporciones de hembras sin colmillos, lo que sugiere que existe un fenotipo
intermedio hereditario (72,73,80); pero mira (29). No obstante, debido a la incertidumbre sobre el origen congénito
versus adquirido de los colmillos únicos para cualquier observación dada y las comparativamente pocas
observaciones de hembras con un colmillo en los datos de la encuesta madre-hijo (norte=0 madres, 8 terneros),
nos abstenemos de intentar inferencias sobre la base y el patrón de herencia de la carencia de colmillos unilateral;
consideramos que esta es una vía importante para futuras investigaciones.
Figura S1. Genotipo LD vs distancia entre SNPs.Desequilibrio de ligamiento (DL) entre genotipos, medido comor2,
decayó a aproximadamente 0,5 (línea de puntos horizontal) para los SNP que estaban separados por 10 kb (línea de
puntos vertical). El panel superior muestra la descomposición der2vs distancia inter-SNP para el cromosoma X,
mientras que el panel inferior muestra este decaimiento para el cromosoma 1.
Figura S2. LD a lo largo de scaffold_0.El panel superior muestra el desequilibrio de ligamiento (LD), medido con vectores
SNP, utilizando ventanas de 50 SNP que tienen un tamaño medio de 7,8 kb. El panel inferior muestra patrones de LD a
mayor escala utilizando ventanas 1000-SNP que tienen un tamaño medio de 343 kb. La región candidata en scaffold_0 (línea
roja), que forma parte del cromosoma 1, está anidada dentro de una región más grande (~6Mb; líneas punteadas) que
muestra una LD elevada.
referencias y notas
13. RL Beyers, JA Hart, ARE Sinclair, F. Grossmann, B. Klinkenberg, S. Dino, Resource wars and
conflict ivory: The impact of civil conflict on elephants in the Democratic
República del Congo—el caso de la Reserva Okapi.MÁS UNO6, e27129 (2011).
doi:10.1371/journal.pone.0027129 Medline
14. V. Varun, T. Ewing, “La maldición del marfil: la militarización y profesionalización de la caza furtiva en
África” (Born Free USA, 2014);www.bornfreeusa.org/campaigns/wildlife-trade/ivoryscurse-the-
militarization-and-professionalization-of-poaching-in-africa-2014/ .
15. ME Stalmans, TJ Massad, MJS Peel, CE Tarnita, RM Pringle, Colapso inducido por la guerra y
recuperación asimétrica de poblaciones de grandes mamíferos en el Parque Nacional
Gorongosa, Mozambique.MÁS UNO14, e0212864 (2019).doi:10.1371/journal.pone.0212864
Medline
16. JOJ Hatton, M. Couto, “Biodiversity and war: A case study of Mozambique” (Programa
de Apoyo a la Biodiversidad, 2001).
17. AM Whitehouse, Sin colmillos en la población de elefantes del Parque Nacional de Elefantes
Addo.S.Afr. J. Zool.257, 249–254 (2002).
18. F. Tajima, Método estadístico para probar la hipótesis de mutación neutral por polimorfismo de
ADN.Genética123, 585–595 (1989).doi:10.1093/genética/123.3.585 Medline
19. JH Poole, “Los efectos de la caza furtiva en las estructuras de edad y los patrones sociales y reproductivos de
poblaciones seleccionadas de elefantes de África Oriental” enEl Comercio de Marfil y el Futuro del Elefante
Africano. Volumen 2: Informes técnicos(Grupo de Revisión del Comercio de Marfil, 1989).
20. H. Jachmann, PSM Berry, H. Imae, Sin colmillos en elefantes africanos: una tendencia futura.
Afr. J. Ecol.33, 230–235 (1995).doi:10.1111/j.1365-2028.1995.tb00800.x
21. G. Steenkamp, SM Ferreira, MN Bester, Sin colmillos y fracturas de colmillos en elefantes de
sabana africana (Loxodonta africana) en libertad.JS Af. Veterinario. Asoc.78, 75–80 (2007).
doi:10.4102/jsava.v78i2.294 Medline
23. PC Lee, V. Fishlock, CE Webber, CJ Moss, Las ventajas reproductivas de una larga vida: Longevidad y
senescencia en elefantes africanos hembra salvajes.Comportamiento Ecol. Sociobiol.70, 337–345
(2016).doi:10.1007/s00265-015-2051-5 Medline
24. N. Chhabra, M. Goswami, A. Chhabra, Bases genéticas de la agenesia dental: genética molecular que
modela la odontología clínica.Medicina. Patol Oral. Cirugía Oral. bucal19, e112–e119 (2014).
doi:10.4317/medoral.19158 Medline
25. M. Varela, MJ Trujillo-Tiebas, P. García-Camba, Gemelos idénticos que revelan hipodoncia discordante.
La justificación de las diferencias de los arcos dentales en gemelos monocigóticos.EUR. Arco.
pediatra Mella.12, 318–322 (2011).doi:10.1007/BF03262831 Medline
26. N. Alachiotis, P. Pavlidis, RAiSD detecta la selección positiva basada en múltiples firmas de un
barrido selectivo y vectores SNP.común Biol.1, 79 (2018).doi:10.1038/s42003-018-0085-8
Medline
27. KE Holsinger, BS Weir, Genética en poblaciones estructuradas geográficamente: definición,
estimación e interpretación de F (ST).Nat. Rev. Genet.10, 639–650 (2009). doi:10.1038/
nrg2611 Medline
28. M. Nei, JC Miller, Un método simple para estimar el número promedio de sustituciones de
nucleótidos dentro y entre poblaciones a partir de datos de restricción.Genética125, 873–879
(1990).doi:10.1093/genetics/125.4.873 Medline
34. S. Malik, S. Gupta, V. Wadhwan, GP Suhasini, Displasia de dentina tipo I: una entidad rara.J.
Oral Maxilofacial. Patol.19, 110 (2015).doi:10.4103/0973-029X.157220 Medline
35. P. Arnold, L. Koopmann, F. Peters, F. Birkenfeld, SV-L. Goff, T. Damm, C. Qin, C. Moali,
R. Lucius, C. Becker-Pauly, La deficiencia de las enzimas meprin α y meprin β que escinden
DSPP no produce malformación de la dentina en ratones.Res. de tejido celular.367, 351–358
(2017).doi:10.1007/s00441-016-2498-3 Medline
36. SR, J. Sarukhan, A. Whyte, R. Hassan, “Evaluación de ecosistemas del milenio: ecosistemas y
bienestar humano” (Island Press, 2006).
37. CJ Moss,Memorias de elefantes: trece años en la vida de una familia de elefantes(William
Morrow, 1988).
38. IJ Whyte, A. Hall-Martin, Características de crecimiento de colmillos de elefantes en el Parque Nacional
Kruger.Paquidermo55, 31–40 (2018).
39. E. Abe, Tusklessness entre los elefantes del Parque Nacional Queen Elizabeth, Uganda.
Paquidermo22, 46–47 (1996).
40. JS Brashares, B. Abrahms, KJ Fiorella, CD Golden, CE Hojnowski, RA Marsh, D.
J. McCauley, TA Nuñez, K. Seto, L. Withey, Decadencia de la vida silvestre y conflicto social.
Ciencias345, 376–378 (2014).doi:10.1126/science.1256734 Medline
41. JBS Haldane, La eliminación selectiva de los zorros plateados en el este de Canadá.J. Genet.44,
296–304 (1942).doi:10.1007/BF02982833
42. KH Andersen, K. Brander, La tasa esperada de evolución inducida por la pesca es lenta.proc. nacional
Academia ciencia EE.UU106, 11657–11660 (2009).doi:10.1073/pnas.0901690106 Medline
43. M. Heino, B. Díaz Pauli, U. Dieckmann, Evolución inducida por la pesca.año Rev. Ecol. Evol. sist.46
, 461–480 (2015).doi:10.1146/annurev-ecolsys-112414-054339
44. FW Allendorf, JJ Hard, Evolución inducida por humanos causada por selección no natural a través de la cosecha
de animales salvajes.proc. nacional Academia ciencia EE.UU106(suplementario 1), 9987–9994 (2009).
doi:10.1073/pnas.0901069106 Medline
45. JS Weir, Propiedades químicas y ocurrencia en la arena de Kalahari de lameduras de sal creadas
por elefantes.J. Zool.158, 293–310 (1969).doi:10.1111/j.1469-7998.1969.tb02148.x
53. CJ Moss, H. Croze, PC Lee, Eds., “Los elefantes de Amboseli: Introducción” enLos elefantes de
Amboseli: una perspectiva a largo plazo de un mamífero longevo(Universidad de Chicago
Press, 2011).
55. Sikes RS; Comité de Cuidado y Uso de Animales de la Sociedad Americana de Mamíferos, 2016
Directrices de la Sociedad Americana de Mamíferos para el uso de mamíferos salvajes en
investigación y educación.J. Mamífero.97, 663–688 (2016). doi:10.1093/jmammal/gyw078
Medline
56. S. Chen, Y. Zhou, Y. Chen, J. Gu, fastp: un preprocesador FASTQ todo en uno ultrarrápido.
Bioinformática34, i884–i890 (2018).doi:10.1093/bioinformatics/bty560 Medline
57. H. Li, R. Durbin, Alineación de lectura corta rápida y precisa con la transformada de Burrows-
Wheeler. Bioinformática25, 1754-1760 (2009).doi:10.1093/bioinformatics/btp324 Medline
58. P. Cingolani, A. Platts, L. Wang, M. Coon, T. Nguyen, L. Wang, SJ Land, X. Lu, DM Ruden, Un
programa para anotar y predecir los efectos de los polimorfismos de un solo nucleótido,
SnpEff .Volar6, 80–92 (2012).doi:10.4161/fly.19695 Medline
59. P. Danecek, A. Auton, G. Abecasis, CA Albers, E. Banks, MA DePristo, RE Handsaker, G. Lunter,
GT Marth, ST Sherry, G. McVean, R. Durbin, Grupo de análisis del proyecto 1000 Genomes ,
El formato de llamada variante y VCFtools.Bioinformática27, 2156–2158 (2011).doi:10.1093/
bioinformatics/btr330 Medline
60. B. Charlesworth, Medidas de divergencia entre poblaciones y el efecto de las fuerzas que
reducir la variabilidad.mol. Biol. Evol.15, 538–543 (1998).
doi:10.1093/oxfordjournals.molbev.a025953 Medline
61. TE Cruickshank, MW Hahn, Reanalysis sugiere que las islas genómicas de especiación se deben a una
diversidad reducida, no a un flujo genético reducido.mol. Ecol.23, 3133–3157 (2014). doi:10.1111/
mec.12796 Medline
63. J. Smith, MR Kronforst, ¿Las especies de mariposas Heliconius intercambian alelos miméticos?Biol.
Letón.9, 20130503 (2013).doi:10.1098/rsbl.2013.0503 Medline
64. BJ Arnold, B. Lahner, JM DaCosta, CM Weisman, JD Hollister, DE Salt, K. Bomblies, L. Yant,
Alelos prestados y convergencia en la adaptación serpentina.proc. nacional Academia
ciencia EE.UU113, 8320–8325 (2016).doi:10.1073/pnas.1600405113 Medline
65. A. Kuzniar, J. Maassen, S. Verhoeven, L. Santuari, C. Shneider, WP Kloosterman, J. de Ridder, sv-callers: un
flujo de trabajo paralelo altamente portátil para la detección de variantes estructurales en datos de
secuencias del genoma completo.PeerJ8, e8214 (2020).doi:10.7717/peerj.8214 Medline
66. DC Jeffares, C. Jolly, M. Hoti, D. Speed, L. Shaw, C. Rallis, F. Balloux, C. Dessimoz, J. Bähler, FJ
Sedlazeck, Las variaciones estructurales transitorias tienen fuertes efectos en los rasgos
cuantitativos y Aislamiento reproductivo en levaduras de fisión.Nat. común8, 14061 (2017).
doi:10.1038/ncomms14061 Medline
Ciencias(ISSN ) es una publicación de la Asociación Estadounidense para el Avance de la Ciencia. 1200 New York Avenue NW, Washington, DC 20005.
El títuloCienciases una marca registrada de AAAS.
Copyright © 2021 Los autores, algunos derechos reservados; Licenciatario exclusivo Asociación Americana para el Avance de la Ciencia. Sin reclamo de obras
originales del gobierno de EE. UU.