11.07.2015 Views

Alvarez-Valín, F. 2000. Evolución molecular - Facultad de Ciencias

Alvarez-Valín, F. 2000. Evolución molecular - Facultad de Ciencias

Alvarez-Valín, F. 2000. Evolución molecular - Facultad de Ciencias

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Evolución <strong>molecular</strong>: neutralismo y seleccionismoEVOLUCIÓN MOLECULAR:NEUTRALISMO Y SELECCIONISMOFernando Álvarez-Valín 1LOS FACTORES QUE GOBIERNAN EL CAMBIO A NIVEL MOLECULAR, ASÍCOMO la variabilidad genética intra-específica, han sido tema <strong>de</strong> <strong>de</strong>bateintenso <strong>de</strong>s<strong>de</strong> finales <strong>de</strong> la década <strong>de</strong> 1960. Al menos dos visionesalternativas han surgido en relación al problema. Por un lado, lo queconocemos como Teoría Neutralista <strong>de</strong> la Evolución Molecular y por otrolado lo que llamaremos visiones Seleccionistas, las cuales engloban unaamplia gama <strong>de</strong> posibilida<strong>de</strong>s. En este capítulo discutiremos los tópicosmás sobresalientes <strong>de</strong> esta controversia.Teoría Neutralista <strong>de</strong> la Evolución MolecularEsta teoría sostiene que la inmensa mayoría <strong>de</strong>l cambio a nivel <strong>molecular</strong>es adaptativamente neutro, es <strong>de</strong>cir, que las nuevas variantes no son nimás ni menos adaptadas que las variantes preexistentes; por tanto sonselectivamente neutras. 2 Como <strong>de</strong>s<strong>de</strong> el punto <strong>de</strong> vista funcional lasvariantes son equivalentes, la selección natural no podría actuarfavoreciendo unas en relación a otras; en consecuencia, otra fuerza distintaa ella <strong>de</strong>bería ser la que gobierna el cambio. En concreto, el neutralismopropone que la mutación genética y la <strong>de</strong>riva genética al azar, son lasfuerzas fundamentales en el proceso <strong>de</strong> cambio a nivel <strong>molecular</strong>. Deacuerdo a esta visión la selección jugaría un papel muy importanteeliminando las variantes <strong>de</strong>letéreas (selección negativa), pero la fijación <strong>de</strong>variantes beneficiosas mediante selección natural seria un evento muypoco frecuente. Para po<strong>de</strong>r enten<strong>de</strong>r esta teoría resulta necesario repasarlos conceptos <strong>de</strong>: mutación, sustitución genética y <strong>de</strong>riva genética.Mutaciones y sustituciones1. Sección Biomatemática, Instituto <strong>de</strong> Biología, <strong>Facultad</strong> <strong>de</strong> <strong>Ciencias</strong>. Iguá 4225. Montevi<strong>de</strong>o 11400.Uruguay. E-mail: falvarez@fcien.edu.uy.2. Esta teoría fue propuesta in<strong>de</strong>pendientemente por Motoo Kimura (1968): Evolutionary rate at the<strong>molecular</strong> level, Nature 217: 624-626; y por King JL & Jukes TH (1969): Non-Darwinian evolution,Science 164: 788-798.243


Fernando Álvarez-ValínUna mutación es cualquier cambio en el material genético que surge enun gameto. Estos cambios pue<strong>de</strong>n ser <strong>de</strong> distinto tipo, como mutacionespuntuales (cambio <strong>de</strong> una base nitrogenada por otra), <strong>de</strong>leción o inserción<strong>de</strong> un segmento <strong>de</strong> ADN o inversiones. En general, las mutacionespuntuales son el tipo <strong>de</strong> cambio más común. Es importante <strong>de</strong>stacar quecuando una mutación surge, la misma está presente en un solo individuoen la población (aquél que se originó <strong>de</strong>l gameto mutante). A<strong>de</strong>más, para elgen mutante, el individuo es heterocigoto, pues es altamente improbableque el gameto que recibió <strong>de</strong>l otro parental también haya mutado enexactamente el mismo sitio. Por tanto po<strong>de</strong>mos <strong>de</strong>cir que la frecuenciaoriginal <strong>de</strong> cada mutación es ½ N, siendo N el tamaño <strong>de</strong> la población. Esmuy importante resaltar este hecho <strong>de</strong> que en general cada mutaciónpue<strong>de</strong> ser consi<strong>de</strong>rada como un evento único. Esto se basa en el siguienterazonamiento simple: para un gen codificante <strong>de</strong> una proteína <strong>de</strong> 300aminoácidos <strong>de</strong> largo y que por tanto consiste en 900 nucleótidos, existen4 900 (~10 541 ) formas alternativas <strong>de</strong> ese gen. Este es un número realmentegran<strong>de</strong>, incluso si se lo compara con, por ejemplo, el número <strong>de</strong> átomosque hay en el Universo que es alre<strong>de</strong>dor <strong>de</strong> 10 81 . Por otra parte, dado unalelo particular <strong>de</strong>l gen en cuestión, éste pue<strong>de</strong> transformarse medianteuna mutación que afecte a una sola base en 3x900 (=2700) alelos distintos.Sobre la base <strong>de</strong> este razonamiento, Crow & Kimura propusieron lo que seconoce como sistema <strong>de</strong> alelos infinitos, que básicamente sostiene que lamutación no es un evento recurrente, y en consecuencia por mutaciónningún alelo pue<strong>de</strong> alcanzar una frecuencia superior a ½N. 3¿Cuál es el <strong>de</strong>stino <strong>de</strong> una mutación que surge en la población? Enbuena medida <strong>de</strong>pen<strong>de</strong> <strong>de</strong> si la mutación es <strong>de</strong>letérea (<strong>de</strong>sventajosa),favorable o funcionalmente equivalente en relación al alelo “salvaje”. Si lamutación es <strong>de</strong>sfavorable, el individuo que porte la mutaciónseguramente tenga menor chance <strong>de</strong> sobrevivir o posea menor fertilidad.En dicho caso el <strong>de</strong>stino más probable <strong>de</strong> la mutación es que <strong>de</strong>saparezcarápidamente <strong>de</strong> la población. Esto es lo que llamamos selección naturalnegativa. Si la mutación es ventajosa, el individuo que porta la mutacióntendrá mayor chance o <strong>de</strong> sobrevivir o <strong>de</strong> <strong>de</strong>jar mayor <strong>de</strong>scen<strong>de</strong>ncia. El<strong>de</strong>stino más probable <strong>de</strong> esta mutación (si el coeficiente <strong>de</strong> selección essuficientemente alto) es que se impondrá en la población, o lo que esequivalente alcanzará una frecuencia igual a 1 (todos los individuos <strong>de</strong> lapoblación poseerán la nueva variante). El proceso mediante el cual unamutación se establece en la población es lo que llamamos fijación, y lamutación que se ha fijado se le llama sustitución. En otras palabras lafijación es el proceso mediante el cual una mutación se fija en una especieo en una población. La última posibilidad es que la mutación sea neutra,es <strong>de</strong>cir funcionalmente equivalente en relación con el alelo “salvaje”. El<strong>de</strong>stino para este tipo <strong>de</strong> mutaciones <strong>de</strong>pen<strong>de</strong> exclusivamente <strong>de</strong>l azar. Esprobable que dicha mutación <strong>de</strong>saparezca en la siguiente generación, porejemplo si el individuo que porta la mutación no tiene <strong>de</strong>scen<strong>de</strong>ncia, o sitransfiere en su gameto el alelo no mutante. Sin embargo es posible quela mutación “tenga suerte”, y que sea transmitida a la <strong>de</strong>scen<strong>de</strong>ncia y que3. Kimura M & Crow JF (1964): The number of alleles that can be maintained in a finite population.Genetics 49: 725-738.244


Fernando Álvarez-Valínexpresarse en los siguientes términos:A tiempo infinito: P(fA 1) + P(fA 2) + P(fA 3) + P(fA 4)+.......+ P(fA 2N-1) +P(fA 2N)=1don<strong>de</strong> P(fA i) es la probabilidad <strong>de</strong> fijación <strong>de</strong>l alelo A i.Dado quetenemos queP(fA 1 ) = P(fA 2) = P(fA 3) =............=P(fA 2N-1) = P(fA 2N),P(fA i) = ½N.Si queremos calcular la probabilidad <strong>de</strong> que se fije A 1 o A 2 o A 3, esto va aestar dado por la suma <strong>de</strong> las probabilida<strong>de</strong>s <strong>de</strong> cada uno por separado,puesto que son eventos excluyentes. Es <strong>de</strong>cir: P(fA 1|fA 2|fA 3) = P(fA 1) +P(fA 2) + P(fA 3).Ahora supongamos que consi<strong>de</strong>ramos a A 1, A 2 y A 3 como una mismacosa, llamémosle A 123; tenemos pues que la frecuencia <strong>de</strong> A 123 será igual ala frecuencia <strong>de</strong> A 1 más la frecuencia A 2 más la frecuencia <strong>de</strong> A 3, lo que eslo mismo que 3 / 2N; la probabilidad <strong>de</strong> fijación <strong>de</strong> A 123 será P(fA 1) + P(fA 2) +P(fA 3) = 3 / 2N. Vemos entonces que la probabilidad <strong>de</strong> fijación <strong>de</strong> un alelo encondiciones <strong>de</strong> <strong>de</strong>riva genética es igual a la frecuencia <strong>de</strong> ese alelo en lapoblación. Cabe <strong>de</strong>stacar que para el caso particular <strong>de</strong> una mutaciónrecién aparecida en la población, cuya frecuencia es ½N, también tiene unaprobabilidad <strong>de</strong> fijación <strong>de</strong> ½N.El Reloj MolecularUno <strong>de</strong> los <strong>de</strong>scubrimientos más importantes en lo que concierne a laevolución <strong>de</strong> las moléculas es el que realizaron Emile Zuckerkandl y LinusPauling 4 en 1964 analizaron el grado <strong>de</strong> divergencia aminoacídica enrelación al tiempo <strong>de</strong> divergencia. 5 Como se muestra en la Fig. 1, el grado<strong>de</strong> divergencia aminoacídica incrementa en forma lineal con el tiempo <strong>de</strong>divergencia entre las especies que se comparan. El tiempo <strong>de</strong> divergencia asu vez fue estimado basándose en el registro fósil. Distintas proteínasdivergen a diferente velocidad, pero la tasa <strong>de</strong> divergencia es constante a lolargo <strong>de</strong>l tiempo para cada una <strong>de</strong> ellas. Es interesante por ejemplo que elgrado <strong>de</strong> divergencia entre el hombre y la carpa se produzca a la mismavelocidad que entre hombre y ratón. Si consi<strong>de</strong>ramos que la carpa se hamantenido prácticamente incambiada <strong>de</strong>s<strong>de</strong> el punto <strong>de</strong> vista morfológico4. Emile Zuckerkandl (n.1922), evolucionista austríaco, luego estadouni<strong>de</strong>nse y francés, se doctoró enbioquímica en la Sorbonne en 1959, y <strong>de</strong>s<strong>de</strong> ese año fue colaborador <strong>de</strong>l químico estadouni<strong>de</strong>nseLinus Pauling. Tuvo un rol fundamental en el <strong>de</strong>sarrollo <strong>de</strong> las i<strong>de</strong>as sobre la evolución <strong>molecular</strong>; enparticular fue pionero en la <strong>de</strong>l “reloj <strong>molecular</strong>”, según la cual cada gen o región génica evoluciona auna tasa estadísticamente constante. En 1971 fue nombrado editor jefe <strong>de</strong>l Journal of MolecularEvolution. Pauling (1901-1994) abordó los problemas <strong>de</strong> la química atómica sobre la base <strong>de</strong> lamecánica cuántica, y precisó la naturaleza <strong>de</strong> las ligaduras químicas y la estructura <strong>de</strong> las moléculas.También colaboró en el <strong>de</strong>scubrimiento <strong>de</strong> una “enfermedad <strong>molecular</strong>” <strong>de</strong> la hemoglobina. Ganó suPremio Nobel en 1954. Sus campañas por el control <strong>de</strong> armas nucleares y contra las pruebas nucleareslo enemistaron con varios colegas, pero le merecieron también el Premio Nobel <strong>de</strong> la Paz en 1962. (N.<strong>de</strong> E.)5. Zuckerkandl E & Pauling L (1965): Evolutionary divergence and convergence in proteins, pp. 97-166 <strong>de</strong> Bryson V & Hogel J (eds.): Evolving genes and proteins, Aca<strong>de</strong>mic Press, New York.246


Evolución <strong>molecular</strong>: neutralismo y seleccionismoen los últimos 400 millones <strong>de</strong> años mientras que los mamíferos hansufrido una gran diferenciación morfológica en los últimos 80 millones <strong>de</strong>años, uno esperaría que el linaje que conduce a la carpa evolucione máslentamente que el <strong>de</strong> los mamíferos. Sin embargo, las comparaciones <strong>de</strong> lassecuencias muestran que el grado <strong>de</strong> divergencia es proporcional al tiempo.¿Cómo explica la teoría neutralista el reloj <strong>molecular</strong>?La tasa <strong>de</strong> sustituciones K <strong>de</strong> alelos selectivamente equivalentes estaría<strong>de</strong>terminada por los siguientes parámetros: el número <strong>de</strong> nuevos mutantesque en cada generación son introducidos en la población, y la probabilidad<strong>de</strong> fijación (es <strong>de</strong>cir <strong>de</strong> transformarse en una sustitución) <strong>de</strong> esos nuevosmutantes; <strong>de</strong> esta manera tenemos que K = n° nuevos mutantes xprobabilidad <strong>de</strong> fijación. El número <strong>de</strong> nuevos mutantes está <strong>de</strong>terminadopor la tasa <strong>de</strong> mutación así como por el tamaño <strong>de</strong> la población. Es <strong>de</strong>ciren cada generación aparecen u2N nuevos mutantes (siendo u la tasa <strong>de</strong>mutación). Como ya se ha visto anteriormente, cuando discutimos elmo<strong>de</strong>lo <strong>de</strong> alelos infinitos, cada mutación pue<strong>de</strong> consi<strong>de</strong>rarse como unevento único, por lo que cada nuevo alelo tiene una frecuencia original <strong>de</strong>½N. En consecuencia su probabilidad <strong>de</strong> fijación en condiciones <strong>de</strong> <strong>de</strong>rivagenética es también ½N. De esta forma tenemos: K=u2N x ½N=u.En otras palabras: la tasa <strong>de</strong> sustituciones K es igual a la tasa <strong>de</strong>mutación u. Vemos entonces que K <strong>de</strong>pen<strong>de</strong> exclusivamente <strong>de</strong> u y escompletamente in<strong>de</strong>pendiente <strong>de</strong>l tamaño <strong>de</strong> la población. Esta afirmaciónpue<strong>de</strong> parecer paradójica, pues es <strong>de</strong> esperar que en una población pequeñala <strong>de</strong>riva genética tenga mayor inci<strong>de</strong>ncia y por lo tanto la tasa <strong>de</strong>sustituciones <strong>de</strong>bería ser mayor. Debemos tener en cuenta sin embargo,que si bien en poblaciones chicas la <strong>de</strong>riva genética tiene mayor inci<strong>de</strong>ncia,y <strong>de</strong> hecho la probabilidad <strong>de</strong> fijación es mayor (dado que ½N es mayor amenor N), en estas poblaciones el número <strong>de</strong> nuevos mutantes que sonintroducidos en cada generación es menor, y por tanto ambos efectos secancelan mutuamente.247


Fernando Álvarez-Valín1.31.21.1Fibrinopéptidos(6 millones <strong>de</strong> años)Globinas(6 millones <strong>de</strong> años)Número <strong>de</strong> sustituciones por sitio1.00.90.80.70.60.50.40.3Citocromo C(20 millones <strong>de</strong> años)0.20.10.00 200 400 600 800 1000 1200 1400Millones <strong>de</strong> años<strong>de</strong> divergenciaHistona IV(600 millones <strong>de</strong> años)Radiación adaptativa<strong>de</strong> los mamíferosSeparación entremamíferos y reptilesSeparación entrecordados e insectosSeparación entreplantas y animalesFigura 1.- Relojes <strong>molecular</strong>es en cuatro proteínas. En esta figura se ha graficado el grado<strong>de</strong> divergencia aminoacídica (número <strong>de</strong> cambios por sitio) versus el tiempo <strong>de</strong> divergenciaentre las especies que se comparan. El tiempo <strong>de</strong> divergencia se estima a partir <strong>de</strong>l registrofósil. Se incluye a<strong>de</strong>más la línea <strong>de</strong> regresión para cada una <strong>de</strong> las proteínas analizadas. Enalgunos casos también se incluye una barra que representa los límites <strong>de</strong> confianza ( ) <strong>de</strong>la estimación <strong>de</strong> divergencia. Esta barra se incluye con el objetivo <strong>de</strong> dar una i<strong>de</strong>a <strong>de</strong>l errorexperimental. Los cuadrados ( ) representan a los fibrinopéptidos, los triángulosorientados hacia abajo ( ) a las globinas, los triángulos orientados hacia arriba ( ) alcitocromo C y los círculos ( ) a la histona IV. 66. Modificado <strong>de</strong> Dickerson (1972): The structure and history of an ancient protein, ScientificAmerican, abril.248Separación entrepeces teleósteos y crondícteos


Evolución <strong>molecular</strong>: neutralismo y seleccionismoRestan por explicar sin embargo dos aspectos. En primer lugar: a qué se<strong>de</strong>be que la tasa sea constante por unidad <strong>de</strong> tiempo en lugar <strong>de</strong> constantepor generación. Téngase en cuenta que la tasa <strong>de</strong> mutación mi<strong>de</strong> el número<strong>de</strong> mutaciones por generación. Organismos con tiempos <strong>de</strong> generación máscortos tienen una tasa <strong>de</strong> mutaciones por año superior a aquellosorganismos con tiempo <strong>de</strong> generación más prolongado. Dicho <strong>de</strong> otraforma: si la tasa <strong>de</strong> mutación por generación es la misma para dosorganismos que presenten tiempos <strong>de</strong> generación diferentes, esperaríamosque aquel que posee duración generacional más corta evolucione másrápidamente, puesto que por unidad <strong>de</strong> tiempo (digamos un millón <strong>de</strong>años) tendremos un mayor número <strong>de</strong> generaciones <strong>de</strong> éstos que <strong>de</strong> losorganismos que presentan duración <strong>de</strong> generación mayor. Una posibleexplicación a esta inconsistencia <strong>de</strong> la teoría neutralista la veremos cuandodiscutamos la teoría <strong>de</strong> los alelos levemente <strong>de</strong>letéreos. En segundolugar, resta por explicar a qué se <strong>de</strong>be que distintas proteínas evolucionena diferentes velocida<strong>de</strong>s. El parámetro u, que habíamos <strong>de</strong>finido como tasa<strong>de</strong> mutación <strong>de</strong>l gen por generación, pue<strong>de</strong> en realidad <strong>de</strong>scomponerse endos partes: v la tasa <strong>de</strong> mutación por gen, y f la proporción <strong>de</strong> sitios quepue<strong>de</strong>n aceptar mutaciones sin que se altere la función <strong>de</strong> la proteína,siendo u=vf. De esta manera u representa la proporción <strong>de</strong> nuevasmutaciones que son selectivamente neutras, o lo que es lo mismo, u es latasa <strong>de</strong> mutaciones neutras. El parámetro f, si se espera que cambiesustancialmente <strong>de</strong> una proteína a otra. En base a la ecuación dada arribapodríamos reescribir K=vf, con lo que esperamos que la tasa <strong>de</strong>sustituciones, K, sea constante para cada gen, mientras que se espera queK varíe <strong>de</strong> una proteína a otra en virtud <strong>de</strong> que f sí varía <strong>de</strong> un gen a otro.Debemos tener en cuenta que si f es la proporción <strong>de</strong> aminoácidos quepue<strong>de</strong>n variar, 1-f representa la proporción <strong>de</strong> aminoácidos don<strong>de</strong> no seaceptan cambios, ya sea porque los cambios en dichos aminoácidos alteranla estructura o la función <strong>de</strong> la proteína. En otras palabras 1-f es laproporción <strong>de</strong> aminoácidos funcionalmente importantes que no aceptancambios y por lo tanto <strong>de</strong>berían encontrarse conservados cuandocomparamos las secuencias <strong>de</strong> esta proteína entre diferentes especies.Esto nos lleva al siguiente punto: cómo clasificamos los distintos tipos<strong>de</strong> mutaciones, y con qué frecuencia esperamos que aparezcan cada una <strong>de</strong>ellas.Tipos <strong>de</strong> mutaciones y sus frecuencias <strong>de</strong> apariciónComo ya ha sido mencionado anteriormente, las mutaciones pue<strong>de</strong>n serneutras (adaptativamente equivalentes al alelo salvaje), <strong>de</strong>sventajosas, oventajosas. Los dos últimos tipos <strong>de</strong> mutaciones son afectados por laselección natural.Ahora po<strong>de</strong>mos plantearnos la siguiente pregunta: dado un gencualquiera ¿cuál es la frecuencia esperada para cada uno <strong>de</strong> estos trestipos <strong>de</strong> mutaciones?. Como veremos, éste es el punto neurálgico en lapolémica neutralismo-seleccionismo.249


Fernando Álvarez-ValínDes<strong>de</strong> el punto <strong>de</strong> vista <strong>de</strong> la teoría neutralista realizaríamos el siguienterazonamiento: los genes codifican para proteínas cuyas funciones han sidoestablecidas muy tempranamente en la evolución. Por ejemplo, los genes quecodifican para enzimas que catalizan la glucólisis en células humanastambién están presentes en la bacteria Escherichia coli. Vemos pues que losgenes codifican para proteínas cuyas funciones ya fueron establecidas hacevarios cientos (o incluso miles) <strong>de</strong> millones <strong>de</strong> años. En consecuencia, esaltamente improbable que una mutación nueva pueda introducir algunamejora en estas proteínas cuya funcionalidad ha sido probada por un período<strong>de</strong> tiempo tan prolongado. Mucho más probable, en cambio, es que lasnuevas mutaciones disminuyan o incluso eliminen la funcionalidad <strong>de</strong> laproteína; es <strong>de</strong>cir, esperaríamos que la gran mayoría sean <strong>de</strong>letéreas. En estesentido, cabe resaltar que los estudios experimentales muestran queefectivamente la gran mayoría <strong>de</strong> las nuevas mutaciones son <strong>de</strong>letéreas.Algunas mutaciones podrían no alterar la funcionalidad <strong>de</strong> la proteína, porejemplo los cambios entre aminoácidos similares <strong>de</strong>s<strong>de</strong> el punto <strong>de</strong> vistafísico-químico (por ejemplo ácido aspártico–ácido glutámico), o aquellasmutaciones que se ubiquen en regiones no esenciales <strong>de</strong> la proteína. Estasmutaciones son, por <strong>de</strong>finición, mutaciones neutras. El panorama que quedaplanteado es entonces el siguiente: la gran mayoría <strong>de</strong> la nuevas mutacionesserían <strong>de</strong>letéreas; éstas, claro, son eliminadas por la selección naturalnegativa y por lo tanto no llegan a fijarse, y por tanto no contribuyen a laevolución. Una proporción bastante menor, sería adaptativamenteequivalente al alelo salvaje (neutras); éstas pue<strong>de</strong>n llegar a fijarse mediante<strong>de</strong>riva genética, y por lo tanto pue<strong>de</strong>n contribuir a la evolución. Por último lasmutaciones ventajosas, que también contribuyen a la evolución pues pue<strong>de</strong>nllegar a fijarse al ser favorecidas (selección positiva), surgirían con frecuenciasextremadamente bajas.Dado que únicamente las mutaciones ventajosas y las neutras pue<strong>de</strong>ncontribuir a la evolución, el eje <strong>de</strong> la polémica neutralismo-seleccionismotiene que ver con la proporción relativa entre mutaciones (y sustituciones)que efectivamente pue<strong>de</strong>n contribuir a la evolución, es <strong>de</strong>cir la proporciónrelativa entre variantes neutras y favorables. La visión seleccionista planteaque la frecuencia <strong>de</strong> aparición <strong>de</strong> mutaciones favorables no es tan bajacomo predice la teoría neutralista. De hecho, se han presentado evi<strong>de</strong>nciasindicando que en varios genes la frecuencia <strong>de</strong> sustituciones que pue<strong>de</strong>natribuirse a la selección positiva es inusitadamente alto. Más a<strong>de</strong>lanteveremos algunos ejemplos <strong>de</strong> selección positiva.Evi<strong>de</strong>ncias <strong>de</strong> la existencia <strong>de</strong> las sustituciones selectivamenteneutrasUna metodología bastante usada para aislar genes consiste en transformarcepas mutantes nulas (una mutación nula es aquella cuyo efecto fenotípicoes la pérdida <strong>de</strong> función) para un gen en particular usando una librería <strong>de</strong>genes que contenga el alelo <strong>de</strong> tipo salvaje. Aquellos clones que hayanrecuperado la función perdida representan o bien mutaciones que250


Evolución <strong>molecular</strong>: neutralismo y seleccionismorestauraron la función, o (lo más común) que han sido transformados conel alelo salvaje. Debido al hecho que nuestro gen problema se encuentrainserto en el vector <strong>de</strong> clonación resulta fácil aislarlo. Esta metodología hasido también usada para aislar genes humanos. Por ejemplo muchos <strong>de</strong> losgenes que participan en las vías <strong>de</strong> reparación <strong>de</strong>l ADN en humanos hansido aislados por su habilidad <strong>de</strong> restaurar la capacidad <strong>de</strong> reparación envarias cepas <strong>de</strong>ficientes <strong>de</strong> la levadura Saccharomyces cerevisiae. Veamosun ejemplo concreto. Algunas cepas <strong>de</strong> Saccharomyces son <strong>de</strong>ficientes en eltransporte hacia el aparato <strong>de</strong> Golgi. 7 Esto se <strong>de</strong>be a una mutación en elgen que codifica la proteína Sec23. 8 El gen humano codificante <strong>de</strong> Sec23 escapaz <strong>de</strong> complementar esta disfunción. Si alineamos la secuenciaaminoacídica <strong>de</strong> la proteína Sec23 entre hombre y levadura, po<strong>de</strong>mos verque existen varias diferencias entre ambas proteínas como lo muestra laFig. 2 (pág. siguiente). Sin embargo el gen humano funciona losuficientemente bien en levaduras como para restaurar la función normal.Esto nos lleva a pensar que las sustituciones que observamos entre ambasespecies no afectan <strong>de</strong> manera fundamental la funcionalidad <strong>de</strong>l gen, <strong>de</strong> locontrario el gen humano sería incapaz <strong>de</strong>HumanoMTTYLEFIQQNEERDGVRFSWNVWPSSRLEATRMVVPVAALFTPLKERPDLPPIQYEPVLLevadura----MDF-ETNEDINGVRFTWNVFPSTRSDANSNVVPVGCLYTPLKEYDELNVAPYNPVV::* : **: :****:***:**:* :*. ****..*:***** :* *:**:HumanoCSRTTCRAVLNPLCQVDYRAKLWACNFCYQRNQFPPSYAGISELNQPAELLPQFSSIEYVLevaduraCSGPHCKSILNPYCVIDPRNSSWSCPICNSRNHLPPQYTNLSQENMPLEL--QSTTIEYI** . *:::*** * :* * . *:* :* .**::**.*:.:*: * * ** * ::***:HumanoVLRGPQMPLIFLYVVDTCMEDEDLQALKESMQMSLSLLPPTALVGLITFGRMVQVHELGCLevaduraTNKPVTVPPIFFFVVDLTSETENLDSLKESIITSLSLLPPNALIGLITYGNVVQLHDLSS. : :* **::*** * *:*::****: *******.**:****:*.:**:*:*..HumanoEGISKSYVFRGTKDLSAKQLQEML-GLSKV-PVTQATRGPQVQ---QPPPSNRFLQPVQKLevaduraETIDRCNVFRGDREYQLEALTEMLTGQKPTGPGGAASHLPNAMNKVTPFSLNRFFLPLEQ* *.:. **** :: . : * *** * . . * *:: *:. * . ***: *:::7. Organoi<strong>de</strong> citoplasmático <strong>de</strong>scubierto en 1909 por el médico, neurólogo y embriólogo italianoCamillo Golgi (1844-1926) gracias a una técnica <strong>de</strong> teñido que él mismo <strong>de</strong>sarrolló. A fin <strong>de</strong> eseaño se le otorgó el Premio Nobel <strong>de</strong> Medicina y Fisiología, compartido con el histólogo españolSantiago Ramón y Cajal (1852-1934). (N. <strong>de</strong> E.)8. Paccaud JP, Reith W, Carpentier JL, Ravazzola M, Amherdt M, Schekman R & Orci L (1996):Cloning and functional characterization of mammalian homologues of the COPII component Sec23,Mol Biol Cell 7: 1535-1546.251


Fernando Álvarez-ValínHumanoIDMNLTDLLGELQRDPWPVPQGKRPLRSSGVALSIAVGLLE-CTFPNTGARIMMFIGGPALevaduraVEFKLNQLLENLSPDQWSVPAGHRPLRATGSALNIASLLLQGC-YKNIPARIILFASGPG:::*.:** :*. * *.** *:****::* **.** **: * : * ***::* .**.HumanoTQGPGMVVGDELKTPIRSWHDIDKDNAKYVKKGTKHFEALANRAATTGHVIDIYA-CALDLevaduraTVAPGLIVNSELKDPLRSHHDIDSDHAQHYKKACKFYNQIAQRVAANGHTVDIFAGC-YD* .**::*..*** *:** ****.*:*:: **. *.:: :*:*.*:.**.:**:* * *HumanoQTGLLEMKCCPNLTGGYMVMGDSFNTSLFKQTFQRVFTKDMHGQFKMGFGGTLEIKTSRELevaduraQIGMSEMKQLTDSTGGVLLLTDAFSTAIFKQSYLRLFAKDEEGYLKMAFNGNMAVKTSKD* *: *** .: *** ::: *:*.*::***:: *:*:** .* :**.*.*.: :***::HumanoIKISGAIGPCVSLN-SKGPCVSENEIGTGGTCQWKICGLSPTTTLAIYFEVVN-QHN-APLevaduraLKVQGLIGHASAVKKTDANNISESEIGIGATSTWKMASLSPYHSYAIFFEIANTAANSNP:*:.* ** . ::: :.. :**.*** *.*. **:..*** : **:**:.* * *Humano---IP-QGGR---GAIQFVTQYQHSSGQRRIRVTTIARNWADAQTQIQNIAASFDQEAAALevaduraMMSAPGSADRPHLAYTQFITTYQHSSGTNRIRVTTVA-N-QLLPFGTPAIAASFDQEAAA* ...* . **:* ****** .******:* * ***********HumanoILMARLAIYRAETEEGPDVLRWLDRQLIRLCQKFGEYHKDDPSSFRFSETFSLYPQFMFHLevaduraVLMARIAVHKAETDDGADVIRWLDRTLIKLCQKYADYNKDDPQSFRLAPNFSLYPQFTYY:****:*:::***::*.**:***** **:****:.:*:****.***:: .******* ::HumanoLRRSSFLQVFNNSPDESSYYRHHFMRQDLTQSLIMIQPILYAYSFSGPPEPVLLDSSSILLevaduraLRRSQFLSVFNNSPDETAFYRHIFTREDTTNSLIMIQPTLTSFSMEDDPQPVLLDSISVK****.**.********:::*** * *:* *:******* * ::*:.. *:****** *:HumanoADRILLMDTFFQILIYHGETIAQWRKSGYQDMPEYENFRHLLQAPVDDAQEILHSRFPMPLevaduraPNTILLLDTFFFILIYHGEQIAQWRKAGYQDDPQYADFKALLEEPKLEAAELLVDRFPLP.: ***:**** ******* ******:**** *:* :*: **: * :* *:* .***:*252


Evolución <strong>molecular</strong>: neutralismo y seleccionismoHumanoRYIDTEHGGSQARFLLSKVNPSQTHNNMYAWGQESGAPILTDDVSLQVFMDHLKKLAVSSLevaduraRFIDTEAGGSQARFLLSKLNPSDNYQDM-ARG--GSTIVLTDDVSLQNFMTHLQQVAVSG*:**** ***********:***:.:::* * * ..: :******** ** **:::***.Figura 2.- Alineamiento <strong>de</strong> la proteína Sec23 entre humano y levadura (Saccharomycescerevisiae). Los asteriscos (*) indican que el aminoácido se ha mantenido incambiado,los dos puntos (:) que se trata <strong>de</strong> una sustitución conservadora (entre aminoácidosbioquímicamente similares), mientras que los puntos indican que se trata <strong>de</strong> unasustitución mediana.complementar la mutación pérdida <strong>de</strong> función. Es <strong>de</strong>cir dichassustituciones son selectivamente neutras.Otra línea <strong>de</strong> evi<strong>de</strong>ncia en favor <strong>de</strong> que las sustituciones selectivamenteneutras son predominantes en muchos genes está dada por lo siguienteobservación. Cuando alineamos y comparamos las secuencias proteicas o<strong>de</strong> ADN entre genes homólogos provenientes <strong>de</strong> distintas especies po<strong>de</strong>mosobservar que existen zonas más conservadas (incambiadas) y zonas másdivergentes. Las zonas conservadas coinci<strong>de</strong>n con los aminoácidosfuncionalmente importantes como lo <strong>de</strong>muestran varios estudios <strong>de</strong>mutaciones. Ciertamente se conoce la secuencia nucleotídica <strong>de</strong> lasversiones mutantes <strong>de</strong> varios alelos que producen alteraciones genéticastanto en el hombre como en otros organismos. Dichas mutaciones se ubicanen su inmensa mayoría en las regiones <strong>de</strong>l gen que se encuentranevolutivamente conservadas, indicando entonces que dichos aminoácidosson funcionalmente importantes pues mutaciones en los mismos dan lugara alelos <strong>de</strong>letéreos. Por otra parte muy pocas (o ninguna en algunos genes)se ubican en regiones evolutivamente divergentes. 9 Esto no quiere <strong>de</strong>cir queel gen no mute en estas zonas, sino que las mismas no producenalteraciones fenotípicas <strong>de</strong>tectables. Resulta claro que estos resultados soncompatibles con la teoría neutralista pues los mismos indican quesolamente las zonas funcionalmente menos importantes evolucionan (don<strong>de</strong>esperamos que las nuevas mutaciones sean neutras o casi neutras),mientras que las zonas funcionalmente importantes tien<strong>de</strong>n a serevolutivamente conservadas por acción <strong>de</strong> la selección natural negativa. Siel gen evolucionara bajo efecto <strong>de</strong> la selección (positiva) sería <strong>de</strong> esperar quelas zonas funcionalmente importantes tendieran a evolucionar más rápido.Selección positivaUn tema central en la polémica neutralismo-seleccionismo tiene que vercon la proporción <strong>de</strong> sustituciones que resultan por efecto <strong>de</strong> la selecciónpositiva. Claro está que este tipo <strong>de</strong> selección solo pue<strong>de</strong> llevar a la fijacióna aquellas mutaciones que son beneficiosas. Como hemos discutidoanteriormente, la frecuencia <strong>de</strong> aparición <strong>de</strong> mutaciones ventajosas se9. Krazewak M & Cooper DN (1996): Mutational processes in pathology and evolution, en M.Jackson, Strachan T & Dover G (eds.): Human genome evolution, BIOS Scientific Publishers,Oxford.253


Fernando Álvarez-Valínespera que sea baja en tanto la función <strong>de</strong>l gen permanezca incambiada.Es improbable introducir mejoras en proteínas que ya funcionan bien. Sinembargo po<strong>de</strong>mos consi<strong>de</strong>rar situaciones en las cuales la función <strong>de</strong>l gen, yla proteína por éste codificada, haya cambiado leve o profundamente en sufunción. También es probable que la selección positiva juegue un papelsignificativo en aquellos genes cuyos productos sean sensibles a loscambios ambientales y/o <strong>de</strong> nicho ecológico. Uno <strong>de</strong> los primeros casosdon<strong>de</strong> se pudo presentar evi<strong>de</strong>ncia rotunda <strong>de</strong> la existencia <strong>de</strong> selecciónpositiva, lo constituye el ejemplo <strong>de</strong>l gen que codifica la lisozima enrumiantes y en los langures, estos últimos pertenecientes al or<strong>de</strong>n <strong>de</strong> losprimates. 10 En estos dos grupos <strong>de</strong> mamíferos el enzima es secretada en elestómago anterior don<strong>de</strong> cumple la función <strong>de</strong> digerir las pare<strong>de</strong>sbacterianas lo cual permite utilizar la celulosa <strong>de</strong>gradada por las bacterias.En otros grupos <strong>de</strong> mamíferos este enzima no se secreta en el estomago.Las secuencias proteicas <strong>de</strong> estos enzimas contienen algunos aminoácidosque están presentes solamente en estos dos grupos <strong>de</strong> mamíferos. Otrosgrupos <strong>de</strong> primates como los homínidos y los gran<strong>de</strong>s monos (apes) carecen<strong>de</strong> dichos aminoácidos. Sin duda no po<strong>de</strong>mos atribuir la similitud que seobserva para esta proteína entre langures y rumiantes a ancestría común(esto es son similares porque el ancestro común entre ambos ya poseíadichos aminoácidos), puesto que esta explicación implicaría asumir que loslangures son más próximos a los rumiantes que a los restantes primates.La explicación más razonable es que en ambos linajes filogenéticos estosaminoácidos se adquirieron in<strong>de</strong>pendientemente, por lo que la lizosimarepresentaría un ejemplo <strong>de</strong> convergencia a nivel <strong>molecular</strong>. Podríamosargumentar que esta convergencia se <strong>de</strong>be a sustituciones neutras que seadquirieron al azar y en paralelo en ambos grupos <strong>de</strong> mamíferos, sinembargo esta alternativa es altamente improbable si tenemos en cuentaque las mismas involucran por lo menos siete sitios aminoacídicos, cadauno <strong>de</strong> los cuales tiene una probabilidad <strong>de</strong> (1/19) <strong>de</strong> ser convergente porazar. La probabilidad conjunto para las siete sustituciones convergentes esmenor <strong>de</strong> 10 -6 . Mucho más razonable, en cambio, es la explicación quepropone que dichos cambios aminoacídicos paralelos son el resultado <strong>de</strong>una misma respuesta adaptativa, es <strong>de</strong>cir darle capacidad a la lizosima <strong>de</strong>funcionar correctamente a bajo pH, como el que se encuentra en elestómago anterior <strong>de</strong> estos animales. La lizosima humana, por ejemplo,funciona en forma muy ineficiente en condiciones <strong>de</strong> aci<strong>de</strong>z. A<strong>de</strong>más se hapodido <strong>de</strong>terminar que algunas <strong>de</strong> estas sustituciones <strong>de</strong> aminoácidosconfieren mejor performance en condiciones ácidas.Lo visto en el párrafo anterior representa un ejemplo <strong>de</strong> respuestaadaptativa a un cambio en el nicho ecológico. Esta situaciónprobablemente sólo afecta a una minoría <strong>de</strong> genes. Sin embargo elsurgimiento <strong>de</strong> nuevas funciones no solamente involucra a cambios en elmedio ambiente, el incremento <strong>de</strong> la complejidad celular y tisular pue<strong>de</strong>ser consi<strong>de</strong>rado como un cambio ambiental a escala <strong>molecular</strong> o celular.10. Stewart CB & Wilson AC (1987): Sequence convergence and functional adaptation of stomachlyzozymes from foregut fermenters, Cold Spring Harbor Symp. Quant. Biol. 52: 891-899.254


Evolución <strong>molecular</strong>: neutralismo y seleccionismo¿Cuál es la base genética que genera la complejidad celular? La repuestaa esta pregunta no es simple, pero po<strong>de</strong>mos afirmar que la aparición ydiversificación <strong>de</strong> las familias multigénicas está estrechamente ligado alincremento <strong>de</strong> la complejidad. Una familia multigénica es un grupo <strong>de</strong>genes que codifican para el mismo tipo <strong>de</strong> proteína (pue<strong>de</strong> que no seangenes codificantes <strong>de</strong> proteínas, por ejemplo los genes ribosomalesconstituyen una familia multigénica). Los distintos miembros <strong>de</strong> unafamilia pue<strong>de</strong>n ser muy similares entre sí o incluso idénticos, pero enmuchos casos los distintos miembros <strong>de</strong> la familia génica codifican paraproteínas relacionadas pero funcionalmente divergentes. Los distintosmiembros <strong>de</strong> una familia multigénica se forman por duplicación <strong>de</strong> genes(crossing-over <strong>de</strong>sigual, transposición, etc.) y posterior diversificación. Elmecanismo evolutivo por el cual se genera la diversidad en una familia noestá bien entendido, pero es muy probable que la selección positivajuegue un rol <strong>de</strong>terminante. Un mo<strong>de</strong>lo ampliamente citado es elsiguiente: luego <strong>de</strong> la duplicación tenemos dos copias <strong>de</strong> un gen, lo cualpermite que una <strong>de</strong> ellas acumule mutaciones libremente (dado que lasmutaciones en el gen extra serían neutras) pues hay un gen redundante.En la gran mayoría <strong>de</strong> los casos el gen redundante pier<strong>de</strong> funcionalidad(es <strong>de</strong>cir, se transforma en lo que conocemos como seudogén) y sigueacumulando cambios hasta que <strong>de</strong>genera completamente. Sin embargo,ocasionalmente este conjunto <strong>de</strong> mutaciones pue<strong>de</strong> conducir a unafunción nueva, la cual pue<strong>de</strong> ser ampliamente favorable, por ejemplo lacreación <strong>de</strong> un nuevo sitio activo o una región <strong>de</strong> interacción con otrasproteínas. Cabe aclarar que bajo este mo<strong>de</strong>lo cada una <strong>de</strong> las mutacionesindividuales sería neutra (o <strong>de</strong>letérea si el gen no estuviera duplicado), loque sería beneficioso es el conjunto <strong>de</strong> mutaciones que crean la nuevafunción. Este mo<strong>de</strong>lo ha sido llamado evolución <strong>de</strong> proteínasfuncionalmente nuevas durante el periodo <strong>de</strong> no-funcionalidad. 11 Sinembargo las evi<strong>de</strong>ncias indican que este mo<strong>de</strong>lo es incorrecto.Ciertamente varios estudios muestran que los genes duplicados no mutanlibremente, <strong>de</strong> hecho la mayoría <strong>de</strong> las mutaciones son <strong>de</strong>letéreas. 12 Elmo<strong>de</strong>lo alternativo (seleccionista) propone que no sólo el conjunto, sinocada uno <strong>de</strong> los cambios individuales que llevan a la nueva función esadaptativamente favorable. Esto pue<strong>de</strong> ser testado <strong>de</strong> la siguiente forma:si la evolución ocurrió mediante sustituciones favorables, éstas <strong>de</strong>beránafectar las regiones funcionalmente importantes <strong>de</strong>l gen, lo cual esexactamente lo contrario a lo que se observaría bajo evolución neutra quepredice que las regiones importantes <strong>de</strong>ben estar muy conservadas. Ensegundo lugar, la velocidad <strong>de</strong> evolución en estas regiones <strong>de</strong>be sersuperior a lo que predice la teoría neutral puesto que la selección positivaes muy eficiente. Ambas predicciones pue<strong>de</strong>n ser testadas. En relación ala primera <strong>de</strong> estas predicciones po<strong>de</strong>mos realizar la siguientecomparación. Si una región es funcionalmente importante, entonces lamisma <strong>de</strong>be estar conservada cuando comparamos entre distintasespecies miembros <strong>de</strong> la familia multigénica que cumplen la mismafunción (por ejemplo cuando comparamos dos globinas ? <strong>de</strong> diferentesespecies), pero la misma <strong>de</strong>bería ser muy divergente cuando comparamosdistintas variantes funcionales <strong>de</strong> la misma familia génica (por ejemplo la11. Ohno S (1973): Ancient linkage groups and frozen acci<strong>de</strong>nts, Nature 244: 259-262.12. Ver Hughes AL (1994): The evolution of functionally novel proteins after gene duplication, Proc. R.Soc. London 256: 119-124.255


Fernando Álvarez-Valíncomparación entre ? y ? globinas). La segunda predicción, involucra lavelocidad <strong>de</strong> evolución relativa a la velocidad <strong>de</strong> evolución neutra, sepue<strong>de</strong> testar <strong>de</strong> la siguiente forma: po<strong>de</strong>mos asumir que los cambiossinónimos (entre codones que codifican el mismo aminoácido) sonselectivamente neutros pues los mismos no producen alteraciones <strong>de</strong>ningún tipo en la secuencia <strong>de</strong> aminoácidos. Por lo tanto la tasa <strong>de</strong>evolución sinónima (K s) <strong>de</strong>bería ser una buena estimación <strong>de</strong> la tasa <strong>de</strong>evolución neutra. Po<strong>de</strong>mos comparar entonces la tasa <strong>de</strong> cambioaminoacídico (K a) en una región que sospechamos se encuentra bajoefecto <strong>de</strong> selección positiva con la tasa <strong>de</strong> cambio sinónimos para lamisma región <strong>de</strong>l gen. Sí el segmento <strong>de</strong>l gen en cuestión estáevolucionando a gran velocidad <strong>de</strong>bido a la acumulación <strong>de</strong> sustitucionesbeneficiosas, entonces K a <strong>de</strong>bería ser mayor que K s. Varios estudios hanconfirmado esta predicción. Numerosas familias multigénicas presentanevi<strong>de</strong>ncias claras <strong>de</strong> selección positiva en regiones <strong>de</strong> sus genes quecodifican segmentos funcionalmente importantes <strong>de</strong> la proteína (ver Tabla1 en pág. siguiente).Por último es importante mencionar que también se ha <strong>de</strong>tectadoselección positiva en aquellos genes y familias multigénicas en los cuales lavariabilidad tiene un valor <strong>de</strong> por sí. En estos casos la selección favorece lageneración <strong>de</strong> la diversidad a nivel aminoacídico, tanto entre los alelos <strong>de</strong>una población, como entre los miembros <strong>de</strong> la familia multigénica. Losejemplos mejor entendidos son aquellos que están relacionados conproteínas <strong>de</strong>l sistema inmunitario, ya sea inmunoglobulinas, receptores <strong>de</strong>linfocitos T, así como en los genes <strong>de</strong>l sistema <strong>de</strong> histocompatibilidad. Enestos genes las llamadas regiones variables, que participan en lainteracción con antígenos (lo que permite el reconocimiento <strong>de</strong> éstos), estánsujetas a selección que favorece la generación <strong>de</strong> diversidad. 13 Resulta claroque la diversidad en estas regiones tiene un valor <strong>de</strong> por sí, pues es lamisma diversidad la que permite el reconocimiento <strong>de</strong> una amplia gama <strong>de</strong>antígenos. En el otro lado <strong>de</strong> la moneda tenemos que muchos parásitos yvirus presentan selección positiva para incrementar la diversidad enaquellos genes codificantes <strong>de</strong> proteínas antigénicas. 14 Esta estrategiapermite a los parásitos evadir o al menos retardar una respuestainmunitaria efectiva <strong>de</strong>l huésped.Tabla 1.- Algunos ejemplos <strong>de</strong> selección positiva.Genes pertenecientes a familiasmultigénicasInhibidor <strong>de</strong> la serín-proteasaRegión variable <strong>de</strong> las inmunoglobulinasReceptores olfatoriosTipo <strong>de</strong> evi<strong>de</strong>nciaK n >K sK n >K sK n >K s , P R >P C13. Tanaka T & Nei (1989): Positive darwinian selection observed at the variable region genes ofinmunoglobulins. Mol. Biol. Evol., 6: 447-459. Hughes , AL & Yeager M (1998): Naturalselection at major histocompatibility complex loci of vertebrates. Annu. Rev. Genetics, 32: 415-435.14. Hughes, AL (1991): Circumsporozoite protein genes of malaria parasites (Plasmodium spp.):evi<strong>de</strong>nce for positive selection on inmunogenic regions. Genetics, 127: 345-353.256


Evolución <strong>molecular</strong>: neutralismo y seleccionismoHistaminaca<strong>de</strong>na ? <strong>de</strong> la integrinaCisteín proteinasaRibonucleasas <strong>de</strong> primates (genes ECP y EDN)Genes <strong>de</strong>l complejo mayor <strong>de</strong> histocompatibilidadK n >K sP R >P CP R >P CK n >K sK n >K sGenes en los que ha cambiado la funcióny genes codificantes <strong>de</strong> proteínas antigénicasAlcohol <strong>de</strong>shidrogenasa en DrosophilaInterleucina 2Enzimas proteolíticas (Calicreina, ? 1 -antripsina, Serpina)Proteínas <strong>de</strong>l Circunesporozoito <strong>de</strong> PlasmodiumGen nef <strong>de</strong>l virus VIHs N n / s N s > p N n / p N sK n >K sK n >K sK n >K sK n >K sNota: K n >K s se refiere al hecho <strong>de</strong> que al menos en algunas regiones <strong>de</strong>l gen, la tasa <strong>de</strong>sustituciones aminoacídicas (K n ) es superior a la tasa <strong>de</strong> cambio sinónimo, mientras queP R >P C , se refiere al hecho <strong>de</strong> que el número <strong>de</strong> cambios aminoacídicos radicales (es <strong>de</strong>cirentre aminoácidos bioquímicamente disímiles ) es superior al <strong>de</strong>l número <strong>de</strong> cambiosconservadores. Por último, s N n / s N s > p N n / p N s significa que la relación entre el número <strong>de</strong>cambios aminoacídicos sobre el número <strong>de</strong> cambios sinónimos es mayor en los sitios quepresentan estas diferencias inter-específicamente (sustituciones) que en los sitios quepresentan la diferencia intra-específicamente (polimorfismos).Cambios sinónimos: ¿paradigma <strong>de</strong> la evolución neutra?Una <strong>de</strong> las primeras predicciones <strong>de</strong> la teoría neutralista era que loscambios sinónimos (entre codones que codifican el mismo aminoácido)<strong>de</strong>berían estar completamente exentos <strong>de</strong> selección natural. Estapresunción estaba basada en el simple hecho <strong>de</strong> que dichos cambios noalteran la naturaleza codificante <strong>de</strong> los genes al no implicar modificaciones<strong>de</strong> ningún tipo en la estructura primaria <strong>de</strong> las proteínas por elloscodificadas. Debido a esta naturaleza supuestamente inocua <strong>de</strong> lasmutaciones sinónimas, era <strong>de</strong> esperar que las diferencias <strong>de</strong> índolesinónima se acumularan a una tasa muy alta en el proceso <strong>de</strong>diferenciación entre genes y especies. Esta afirmación <strong>de</strong> la teoríaneutralista representa un caso particular <strong>de</strong> uno <strong>de</strong> sus postuladosbásicos: la tasa <strong>de</strong> cambio a nivel <strong>molecular</strong> es inversamente proporcionalal grado <strong>de</strong> restricción funcional.Sin embargo, a medida que se acumularon datos <strong>de</strong> secuenciasnucleotídicas, resultó evi<strong>de</strong>nte que los distintos codones sinónimosaparecen en los genes con frecuencias que claramente se apartan <strong>de</strong> unadistribución al azar. Este fenómeno conocido como “uso <strong>de</strong> codonessinónimos” se observa tanto en organismos procariotas como eucariotas. 1515. Grantham R, Gautier C, Gouy M, Mercier R & Pave R (1980): Codon catalog usage and the genomehypothesis, Nucleic Acids Res. 8: 49-62.257


Fernando Álvarez-ValínPor otro lado, Toshimichi Ikemura (<strong>de</strong>l Instituto Nacional <strong>de</strong> Genética,Mishima, Japón) presentó evi<strong>de</strong>ncia que indica que en Escherichia coli ySaccharomyces cerevisiae el uso <strong>de</strong> codones está, en gran medida,<strong>de</strong>terminado por la abundancia relativa <strong>de</strong> los ARN <strong>de</strong> transferencia (ARNt)que reconocen a los correspondientes codones. 16 Esto es, los codones másfrecuentes (codones “óptimos”) son reconocidos por los ARNt másabundantes. Esta correlación entre abundancia <strong>de</strong> codones y abundancia<strong>de</strong> ARNt es especialmente evi<strong>de</strong>nte en los genes que codifican proteínas queestán presentes en altas concentraciones. Esta observación llevóinmediatamente a postular que en estos microorganismos el uso <strong>de</strong>codones estaría <strong>de</strong>terminado por la selección para incrementar la eficienciatraduccional. Este incremento en la eficiencia <strong>de</strong> la síntesis proteica sería<strong>de</strong>bido a que el número medio <strong>de</strong> interacciones entre el ARNt y el sitio-A <strong>de</strong>lribosoma se reduce consi<strong>de</strong>rablemente en un sistema que posea sesgo en eluso <strong>de</strong> codones y sesgo en las poblaciones <strong>de</strong> ARNt en relación a unsistema en el cual todos los codones y ARNt son equiprobables. Comoresultado <strong>de</strong> este <strong>de</strong>cremento en el número <strong>de</strong> interacciones, se reduce eltiempo medio <strong>de</strong> espera (<strong>de</strong>l aminoacil-ARNt correcto) por aminoácido, loque a su vez resulta en una reducción neta <strong>de</strong>l número <strong>de</strong> ribosomasnecesarios para producir una <strong>de</strong>terminada cantidad <strong>de</strong> proteína por unidad<strong>de</strong> tiempo.El uso sesgado <strong>de</strong> codones acompañado con sesgos en las poblaciones <strong>de</strong>ARNt también ha sido implicado en el incremento <strong>de</strong> la fi<strong>de</strong>lidadtraduccional. Por un lado, se ha <strong>de</strong>mostrado experimentalmente en E. colique la sustitución <strong>de</strong> un codón mayor (es <strong>de</strong>cir reconocido por un ARNtabundante) por un codón menor, produce un incremento <strong>de</strong> casi 10 veces enla tasa <strong>de</strong> errores traduccionales en el aminoácido don<strong>de</strong> se realizó lasustitución. 17 Por otro lado Hiroshi Akashi, analizando genes <strong>de</strong> Drosophilamelanoganster y D. simulans mostró que los aminoácidos funcionalmenteimportantes (en general pertenecientes a motivos conservados o dominiosproteicos cuya función se sabe que es importante) tien<strong>de</strong>n a estar codificadospor codones óptimos en una frecuencia significativamente más alta que losaminoácidos no conservados. 18La posible existencia <strong>de</strong> selección en las posiciones sinónimas basada enla evi<strong>de</strong>ncia aportada por el uso <strong>de</strong> codones no azaroso y su vinculacióncon la disponibilidad <strong>de</strong> ARNts isoaceptores, llevó a postular que laevolución sinónima <strong>de</strong>bería estar sujeta a presión selectiva.Posteriormente Paul M. Sharp & Li Wen-Hsiung mostraron que enenterobacterias (E. coli, Salmonella typhimurium) la tasa <strong>de</strong> sustitucionessinónimas es inversamente proporcional al grado <strong>de</strong> sesgo en el uso <strong>de</strong>codones, esto es, genes con mayor sesgo en sus preferencias <strong>de</strong> codones (y16. Ikemura T (1981): Correlation between the abundance of Escherichia coli transfer RNAs and theoccurrence of the respective codons in proteins genes, J. Mol. Biol. 146: 1-21; y (1982): Correlationbetween the abundance of yeast transfer RNAs and the ocurrence of the respective codons in proteingenes, J. Mol. Biol. 158: 573-587.17. Precup J & Parker J (1987): Missense misreading of asparagine codons as a function of codoni<strong>de</strong>ntity and context, J. Biol. Chem. 262: 11351-11356.18. Akashi H (1994): Synonymous codon usage in Drosophila melanogaster, natural selection andtranslational accuracy. Genetics, 136: 927-935.258


Evolución <strong>molecular</strong>: neutralismo y seleccionismopor tanto expresados a altos niveles y con mayor <strong>de</strong>pen<strong>de</strong>ncia <strong>de</strong> lapoblación <strong>de</strong> ARNts), evolucionan (a nivel sinónimo) significativamente máslentamente que aquellos genes con menor sesgo (y por lo tanto expresadosmás débilmente y con menor <strong>de</strong>pen<strong>de</strong>ncia <strong>de</strong> los ARNts). 19 Resultadossimilares han sido obtenidos en otros organismos como Drosophila,Caenorhabditis y Mycobacterium. 20Por otra parte, varios autores han reportado in<strong>de</strong>pendientemente que lastasas <strong>de</strong> sustituciones sinónimas y no-sinónimas no son in<strong>de</strong>pendientes.Por el contrario, existe una correlación relativamente alta (r entre 0.5 y 0.6)y muy significativa entre ellas. 21 Esta correlación ha sido atribuida adistintas causas. Por un lado Ken Wolfe & Paul Sharp sostienen que lamisma podría ser explicada como el resultado <strong>de</strong> mutaciones en dobletes.Esto es, dos bases consecutivas que mutan simultáneamente comoresultado <strong>de</strong> un único evento. Sin duda, si estas mutaciones fueranfrecuentes podrían explicar en parte la correlación en cuestión, puesto queal cambiar dos bases consecutivas, en aproximadamente ½ <strong>de</strong> casos seproduciría una mutación sinónima y una no-sinónima. Sin embargoDominique Mouchiroud, C. Gautier & Giorgio Bernardi han <strong>de</strong>mostrado quesi se eliminan <strong>de</strong> los alineamientos las sustituciones consecutivas en lasposiciones <strong>de</strong>l codón 2-3 y 3-1 (es <strong>de</strong>cir sustituciones sinónimas y nosinónimasconsecutivas que podrían haber sido originadas a partir <strong>de</strong> unaúnica mutación en doblete) la correlación baja pero se mantiene altamentesignificativa. Como explicación alternativa, los mismos autores propusieronque la correlación en cuestión <strong>de</strong>bería ser la consecuencia <strong>de</strong> restriccionesfuncionales comunes a los cambios sinónimos y no-sinónimos. Un posiblevínculo entre ambos tipos <strong>de</strong> sustituciones (es <strong>de</strong>cir, restriccionesfuncionales comunes) pue<strong>de</strong> estar dado por la fi<strong>de</strong>lidad traduccional puestoque los aminoácidos más importantes <strong>de</strong>s<strong>de</strong> el punto <strong>de</strong> vista funcionaltien<strong>de</strong>n a presentar mayor conservación evolutiva como ya ha sido visto ya<strong>de</strong>más tien<strong>de</strong>n a estar codificados por codones mayores para disminuir latasa <strong>de</strong> errores traduccionales en estos codones. Dada esta situación es <strong>de</strong>esperar que estos codones mayores tiendan a mantenerse incambiados(evolutivamente conservados), dado que sustituciones en los mismosimplicaría pasar <strong>de</strong> un codón mayor a uno menor con el consiguienteaumento en la tasa <strong>de</strong> errores traduccionales. Las correlaciones entre lastasas <strong>de</strong> sustituciones sinónimas y no-sinónimas también se observan anivel intragénico. Estudios en genes <strong>de</strong> mamíferos y gramíneas muestran19. Sharp PM & Li W-H (1987): The rate of synonymous substitutions in enetrobacterial genes isinversely related to codon usage bias, Mol. Biol. Evol. 4: 222-230.20. Respectivamente: Sharp PM & Li W-H (1988): On the rate of DNA sequence evolution inDrosophila, J. Mol. Evol. 28: 398-402; Stenico M, Lloyd AT & Sharp PM (1994): Codon usage inCaenorhabditis elegans: <strong>de</strong>lineation of translational selection and mutational biases, Nucleic AcidRes. 22: 2437-2446; <strong>de</strong> Miranda A, Álvarez-Valín F, Jabbari K, Degrave WM & Bernardi G(2000): Gene expression, amino acid conservation and hydrophobicity are the main factors shapingcodon preferences in Mycobacterium tuberculosis and Mycobacterium leprae, J. Mol Evol. 50: 45-55.21. Mouchiroud D, Gautier C & Bernardi G (1995): Frequencies of synonymous substitutions inmammals are gene-specific and correlated with frequencies of non-synonymous substitutions, J.Mol. Evol. 40: 107-113; Ohta T & Ina Y (1995): Variation in synonymous substitutions ratesamong mammalian genes and correlations between synonymous and nonsynosymous divergences,J. Mol. Evol. 41: 717-720; Wolfe KH & Sharp PM (1993): Mammalian gene evolution: nucleoti<strong>de</strong>sequence divergence between mouse and rat, J. Mol. Evol. 37: 441-456.259


Fernando Álvarez-Valínque el número <strong>de</strong> genes que presentan coeficientes <strong>de</strong> correlaciónestadísticamente significativos es mayor <strong>de</strong> lo que se esperaría por azar. 22Cuando hablamos <strong>de</strong> correlación intragénica entre ambas tipos <strong>de</strong>sustituciones queremos <strong>de</strong>cir que aquellas regiones <strong>de</strong>l gen que son másconservadas <strong>de</strong>s<strong>de</strong> el punto <strong>de</strong> vista aminoacídico también son másconservadas a nivel sinónimo, mientras que las zonas con mayordivergencia en términos <strong>de</strong> sustituciones <strong>de</strong> aminoácidos también presentan0,80Distancia nucleotídica, CAI0,700,600,500,400,300,200,100,0015 65 115 165 215 265 315 365 415 465 515 565Posición en el genmayor divergencia sinónima (Fig. 3).Figura 3.- Perfiles <strong>de</strong> divergencia no-sinónima (línea gruesa), sinónima (línea <strong>de</strong>lgada) ysesgo en el uso <strong>de</strong> codones (línea punteada) en el gene codificante <strong>de</strong> la metaloproteinasa<strong>de</strong> membrana (GP63) <strong>de</strong> Leishmania. Estos perfiles fueron obtenidos midiendo lafrecuencia <strong>de</strong> cambios <strong>de</strong> tipo no-sinónimo y sinónimo (corregida para sustitucionesmúltiples y paralelas) <strong>de</strong>ntro <strong>de</strong> ventanas móviles <strong>de</strong> 30 codones <strong>de</strong> largo.Resultados más recientes aportan evi<strong>de</strong>ncias claras en favor <strong>de</strong> lahipótesis <strong>de</strong> las restricciones comunes. Por un lado M.L. Chiusano ycolaboradores han reportado que la estructura secundaria <strong>de</strong> las proteínas(<strong>de</strong>finidas en términos <strong>de</strong> alfa hélice, hoja beta y estructuras aperiódicas)presentan tasas <strong>de</strong> cambio sinónimos y no sinónimos diferentes. 23 Por otrolado Álvarez-Valín y colaboradores han mostrado que en el gen codificante<strong>de</strong> la proteína GP63 <strong>de</strong> Leishmania las sustituciones sinónimas y nosinónimas presentan una distribución claramente no aleatoria en laestructura tridimensional <strong>de</strong> la proteína. 2422. Álvarez-Valín F, Jabbari K & Bernardi G (1998): Synonymous and nonsynonymous substitutionsin mammalian genes: intragenic correlations, J. Mol. Evol. 46: 37-44; Álvarez-Valín F, Jabbari K,Carels N & Bernardi G (1999): Synonymous and nonsynonymous substitutions in genes fromGraminae: intragenic correlations, J. Mol. Evol. 49: 330-342.23. Chiusano ML, D’Onofrio G, Álvarez-Valín F, Jabbari K, Colonna G & Bernardi G (1999):Correlations of nucleoti<strong>de</strong> substitution rates and base composition of mammalian codingsequences with protein structure, Gene 238: 23-31.24. Álvarez-Valín F, Tort JF & Bernardi G (2000): Non-random spatial distribution of synonymoussubstitutions in the GP63 gene from Leishmania, Genetics, 155: 1683-1692.260


Evolución <strong>molecular</strong>: neutralismo y seleccionismoTodo este cúmulo <strong>de</strong> resultados indica con claridad que las posicionessinónimas <strong>de</strong> los codones, lejos <strong>de</strong> estar exentas <strong>de</strong> presiones selectivas –como se pensó originalmente– están sujetas a una amplia gama <strong>de</strong>restricciones funcionales. Muchas <strong>de</strong> estas restricciones son las mismas quetambién afectan a los cambios aminoacídicos.Teoría <strong>de</strong> los alelos casi neutros levemente <strong>de</strong>letéreosHasta ahora nos hemos manejado con el siguiente esquema simplista<strong>de</strong>l <strong>de</strong>stino <strong>de</strong> los distintos tipos <strong>de</strong> mutaciones: las mutaciones<strong>de</strong>sventajosas son eliminadas <strong>de</strong> la población por efecto <strong>de</strong> la selecciónnatural negativa, las variantes ventajosas en cambio son favorecidas por laselección (selección positiva) y tienen gran chance <strong>de</strong> fijarse en lapoblación. Por último las mutaciones neutras no son afectadas por laselección natural, su <strong>de</strong>stino <strong>de</strong>pen<strong>de</strong> exclusivamente <strong>de</strong> la <strong>de</strong>riva genética.Es importante tener en cuenta sin embargo, que la <strong>de</strong>riva genética tambiéninfluye en el <strong>de</strong>stino <strong>de</strong> las mutaciones no-neutras, sean éstas favorables o<strong>de</strong>sfavorables. En poblaciones pequeñas, por efecto <strong>de</strong> la <strong>de</strong>riva genética,una mutación favorable pue<strong>de</strong> llegar a per<strong>de</strong>rse <strong>de</strong> la población, o una<strong>de</strong>letérea pue<strong>de</strong> llegar a fijarse. Po<strong>de</strong>mos afirmar que una mutación secomporta como neutra si el producto |S*N|


Fernando Álvarez-Valín<strong>de</strong> que la tasa <strong>de</strong> aparición <strong>de</strong> mutaciones efectivamente neutras (es <strong>de</strong>cir quese comportan como tales aunque en un sentido estricto no lo sean) <strong>de</strong>pen<strong>de</strong><strong>de</strong>l tamaño poblacional, siendo esta tasa mayor en poblaciones pequeñas queen poblaciones <strong>de</strong> gran tamaño.Esta teoría tiene implicancias en varios aspectos, pero la fundamental esque se relaciona con el reloj <strong>molecular</strong>. Como acabamos <strong>de</strong> ver, enpoblaciones pequeñas el porcentaje <strong>de</strong> alelos que se comportará comoefectivamente neutros es superior al porcentaje que esperamos parapoblaciones gran<strong>de</strong>s. Visto <strong>de</strong>s<strong>de</strong> otro punto <strong>de</strong> vista, la tasa <strong>de</strong> aparición<strong>de</strong> mutaciones efectivamente neutras es inversamente proporcional altamaño <strong>de</strong> poblaciones, consecuentemente la tasa <strong>de</strong> aparición <strong>de</strong>mutaciones que son “vistas” como <strong>de</strong>letéreas por la selección natural esdirectamente proporcional al tamaño poblacional. Como ya mostramosanteriormente, la tasa <strong>de</strong> sustituciones K <strong>de</strong>pen<strong>de</strong> exclusivamente <strong>de</strong> latasa <strong>de</strong> mutaciones neutras (K=u=vf). Bajo la hipótesis que estamos viendo,u varía en forma inversamente proporcional al tamaño poblacional, <strong>de</strong>manera que esperaríamos que la tasa <strong>de</strong> sustituciones K también seainversamente proporcional al tamaño <strong>de</strong> las poblaciones. Por otro lado,tenemos que existe también una relación inversa entre el tiempo <strong>de</strong>generación y el tamaño poblacional. Organismos que poseen tiempos <strong>de</strong>generación cortos tienen tamaños poblacionales mayores (consi<strong>de</strong>remospor ejemplo la comparación entre elefantes y ratones). De esta formapo<strong>de</strong>mos enten<strong>de</strong>r cómo la tasa <strong>de</strong> sustituciones es proporcional al tiempoy no al número <strong>de</strong> generaciones, ya que en organismos con poblacionespequeñas la tasa <strong>de</strong> mutaciones efectivamente neutras sería mayor (y portanto K sería mayor por generación) pero también es mayor la duración <strong>de</strong>cada generación. En poblaciones pequeñas en cambio la tasa <strong>de</strong>sustituciones K es menor por generación pero tenemos un mayor número<strong>de</strong> generaciones por unidad <strong>de</strong> tiempo. Claramente esta relación daríalugar a un reloj <strong>molecular</strong> que es proporcional al tiempo y no al número <strong>de</strong>generaciones.ConclusionesLa discusión <strong>de</strong> algunos <strong>de</strong> los tópicos más importantes en la polémicaseleccionismo-neutralismo nos permite ver que existen evi<strong>de</strong>ncias en uno yotro sentido. Si bien este no es un tema cerrado, estamos en condiciones <strong>de</strong>afirmar que muchos genes, y particularmente aquellos que codifican parafunciones que se han mantenido incambiadas a lo largo <strong>de</strong>l tiempo,evolucionan <strong>de</strong> acuerdo a las predicciones <strong>de</strong> la teoría neutral. Otro grupo<strong>de</strong> genes en cambio parecen evolucionar bajo la influencia <strong>de</strong> selecciónpositiva. Uno podría preguntarse si estos últimos representan una minoríasiendo lo predominante la evolución neutra. Por el momento no estamos encondiciones <strong>de</strong> respon<strong>de</strong>r a esta pregunta, pero es preciso tener en cuentaque una proporción muy importante <strong>de</strong> los genes <strong>de</strong> vertebrados formafamilias multigénicas en las cuales se observa diferenciación en diversassub-funciones. Es posible que en muchos <strong>de</strong> éstos la selección positivajuegue un rol prepon<strong>de</strong>rante.262


Evolución <strong>molecular</strong>: neutralismo y seleccionismoOtra fuente consultadaKrazewak M & Cooper DN (1996): Mutational processes in pathology an<strong>de</strong>volution. en Human Genome Evolution, ed. Jackson M, Strachan T yDover G. BIOS Scientific Publishers, Oxford.263

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!