Articulo de referencia

Reconstrucción de la secuencia ancestral

La reconstrucción de secuencias ancestrales ( ASR ), también conocida como reconstrucción / resurrección de genes / secuencias ancestrales , es una técnica utilizada en el estud...

La reconstrucción de secuencias ancestrales ( ASR ), también conocida como reconstrucción / resurrección de genes / secuencias ancestrales , es una técnica utilizada en el estudio de la evolución molecular . El método utiliza secuencias relacionadas para reconstruir un gen "ancestral" a partir de un alineamiento de secuencias múltiples . [ 1 ]

El método puede utilizarse para «resucitar» proteínas ancestrales y fue propuesto en 1963 por Linus Pauling y Emile Zuckerkandl . [ 2 ] En el caso de las enzimas, este enfoque se ha denominado paleoenzimología . Algunos esfuerzos iniciales se realizaron en las décadas de 1980 y 1990, liderados por el laboratorio de Steven A. Benner , que demostraron el potencial de esta técnica. [ 3 ] Gracias a la mejora de los algoritmos y de mejores técnicas de secuenciación y síntesis, el método se desarrolló aún más a principios de la década de 2000 para permitir la resurrección de una mayor variedad de genes, y mucho más antiguos. [ 4 ] Durante la última década, la resurrección de proteínas ancestrales se ha desarrollado como una estrategia para revelar los mecanismos y la dinámica de la evolución de las proteínas. [ 5 ]

Principios

Una ilustración de un árbol filogenético y cómo influye en la conceptualización de cómo se lleva a cabo el reconocimiento automático del habla (ASR).
Algoritmo para reconstruir las secuencias ancestrales 1, 2 y 3 (ver figura anterior). La secuencia ancestral 1 se puede reconstruir a partir de B y C, siempre que se disponga de al menos un grupo externo, como D o E. Por ejemplo, las secuencias B y C difieren en la posición 4, pero dado que las secuencias D y E contienen una C en esa posición, es muy probable que la secuencia 1 también la tuviera. La secuencia 3 no se puede reconstruir completamente sin una secuencia externa adicional (la incertidumbre se indica con una "X").

La ASR se basa en la observación de que las especies estrechamente relacionadas tienen secuencias de ADN similares (véase la Figura 2). Por ejemplo, si dos especies difieren en un nucleótido, como A en humanos y G en chimpancés, podemos asumir con seguridad que un ancestro tenía A o G y que este nucleótido mutó en uno de los linajes ("con seguridad" porque es estadísticamente muy improbable que un nucleótido mute y luego vuelva a mutar). ¿Cómo podemos determinar si el ancestro tenía A o G? ¡Observamos uno o más grupos externos! Si tanto los gorilas como los orangutanes tienen A, podemos asumir con seguridad que este era el nucleótido ancestral y que la mutación A→G ocurrió en el linaje que dio origen a los chimpancés (véanse las Figuras 1 y 2).

Verificación experimental. La mayoría de los estudios ASR se realizan in vitro y han revelado propiedades de proteínas ancestrales que parecen ser rasgos evolutivamente deseables, como una mayor termoestabilidad, actividad catalítica y promiscuidad catalítica. Estos datos se han atribuido a artefactos de los algoritmos ASR, así como a ilustraciones indicativas del entorno de la Tierra antigua; a menudo, la investigación ASR debe complementarse con controles extensos (generalmente experimentos ASR alternativos) para mitigar el error algorítmico. No todas las proteínas ASR estudiadas exhiben esta llamada "superioridad ancestral". [ 6 ] El campo naciente de la " bioquímica evolutiva " se ha visto reforzado por el reciente aumento de estudios ASR que utilizan los ancestros como formas de investigar la aptitud de los organismos dentro de ciertos contextos celulares, probando efectivamente proteínas ancestrales in vivo . [ 7 ] Debido a las limitaciones inherentes en este tipo de estudios, principalmente la falta de genomas suficientemente antiguos para encajar estos ancestros, el pequeño repertorio de sistemas modelo de laboratorio bien categorizados y la incapacidad de imitar entornos celulares antiguos; Se han realizado muy pocos estudios de ASR in vivo . A pesar de los obstáculos mencionados anteriormente, las ideas preliminares sobre esta línea de investigación, a partir de un artículo de 2015, han revelado que la "superioridad ancestral" observada in vitro no se recapituló in vivo de una proteína determinada. [ 8 ] ASR representa uno de los pocos mecanismos para estudiar la bioquímica de la era precámbrica de la vida (>541 Ma ) y, por lo tanto, se usa a menudo en " paleogenética "; de hecho, Zuckerkandl y Pauling originalmente concibieron ASR como el punto de partida de un campo que denominaron "paleobioquímica".

Metodología

Primero, se alinean varios homólogos relacionados de la proteína de interés en un alineamiento múltiple de secuencias (MSA). Luego, se construye un árbol filogenético con las secuencias inferidas en los nodos de las ramas. Estas secuencias son las denominadas ancestros. Las secuencias ancestrales se calculan típicamente mediante máxima verosimilitud , aunque también se implementan métodos bayesianos . Dado que los ancestros se infieren a partir de una filogenia, la topología y la composición de esta desempeñan un papel fundamental en las secuencias ASR resultantes. Estas secuencias se comparan y, a menudo, se expresan y estudian varias (~10) por nodo filogenético. ASR no pretende recrear la secuencia real de la proteína/ADN ancestral, sino una secuencia que probablemente sea similar a la que se encontraba realmente en el nodo. Esto no se considera una deficiencia de ASR, ya que se ajusta al modelo de " red neutral " de evolución de proteínas, en el que en las uniones evolutivas (nodos) existía una población de secuencias de proteínas genotípicamente diferentes pero fenotípicamente similares en la población de organismos actual. Por lo tanto, es posible que ASR genere una de las secuencias de la red neutral de un nodo y, si bien puede que no represente el genotipo del último ancestro común de las secuencias actuales, es probable que sí represente el fenotipo. [ 7 ] Esto se ve respaldado por la observación actual de que muchas mutaciones en el sitio no catalítico/funcional de una proteína causan cambios menores en las propiedades biofísicas. Por consiguiente, ASR permite investigar las propiedades biofísicas de proteínas pasadas y es indicativo de genética antigua.

Los métodos de máxima verosimilitud (ML) funcionan generando una secuencia donde se predice que el residuo en cada posición es el más probable de ocupar dicha posición, utilizando una matriz de puntuación (similar a las utilizadas en BLAST o MSA) calculada a partir de secuencias existentes. Los métodos alternativos incluyen la máxima parsimonia (MP), que construye una secuencia basándose en un modelo de evolución de secuencias ; generalmente, la idea de que el número mínimo de cambios de secuencia de nucleótidos representa la ruta más eficiente para que la evolución la siga y, según la navaja de Occam, es la más probable. La MP se considera a menudo el método menos fiable para la reconstrucción, ya que, posiblemente, simplifica demasiado la evolución hasta un grado que no es aplicable a la escala de miles de millones de años.

Los llamados métodos bayesianos se utilizan a veces para complementar los métodos de ML, pero suelen producir secuencias más ambiguas. En ASR, el término "ambigüedad" se refiere a posiciones de residuos donde no se puede predecir una sustitución clara; a menudo, en estos casos, se producen varias secuencias ASR que abarcan la mayoría de las ambigüedades y se comparan entre sí. ML ASR suele necesitar experimentos complementarios para indicar que las secuencias derivadas son más que simples consensos de las secuencias de entrada. Esto es particularmente necesario en la observación de la "superioridad ancestral". [ 9 ] En la tendencia de aumento de la termoestabilidad, una explicación es que ML ASR crea una secuencia consenso de varios mecanismos paralelos diferentes que evolucionaron para conferir una termoestabilidad proteica menor a lo largo de la filogenia, lo que conduce a un efecto aditivo que resulta en una termoestabilidad ancestral "superior". [ 10 ]

Validación experimental . La expresión de secuencias de consenso y ASR paralelas mediante métodos no basados ​​en ML a menudo se requieren para desmantelar esta teoría por experimento. Otra preocupación que plantean los métodos de ML es que las matrices de puntuación se derivan de secuencias modernas y las frecuencias de aminoácidos particulares que se ven hoy pueden no ser las mismas que en la biología precámbrica, lo que resulta en una inferencia de secuencia sesgada. Varios estudios han intentado construir matrices de puntuación antiguas mediante diversas metodologías y han comparado las secuencias resultantes y las propiedades biofísicas de sus proteínas. Si bien estas secuencias modificadas dan como resultado secuencias ASR algo diferentes, las propiedades biofísicas observadas no parecieron variar fuera del error experimental. [ 11 ] Debido a la naturaleza "holística" de ASR y la intensa complejidad que surge cuando se consideran todas las posibles fuentes de error experimental, la comunidad experimental considera que la medida definitiva de la fiabilidad de ASR es la comparación de varias reconstrucciones ASR alternativas del mismo nodo y la identificación de propiedades biofísicas similares. Si bien este método no ofrece una medida estadística y matemática sólida de confiabilidad, se basa en la idea fundamental utilizada en ASR de que las sustituciones individuales de aminoácidos no causan cambios significativos en las propiedades biofísicas de una proteína, un principio que debe cumplirse para poder superar el efecto de la ambigüedad de la inferencia. [ 12 ]

Los candidatos utilizados para ASR a menudo se seleccionan en función de la propiedad de interés particular que se estudia, por ejemplo, la termoestabilidad . [ 6 ] Al seleccionar secuencias de cualquiera de los extremos del rango de una propiedad (por ejemplo, proteínas psicrófilas y proteínas termófilas) pero dentro de una familia de proteínas, ASR puede utilizarse para investigar los cambios de secuencia específicos que confirieron el efecto biofísico observado, como las interacciones estabilizadoras. Consideremos en el diagrama que, si la secuencia 'A' codifica una proteína que era óptimamente funcional a pH neutro y la 'D' en condiciones ácidas, los cambios de secuencia entre '5' y '2' pueden ilustrar la explicación biofísica precisa de esta diferencia. Dado que los experimentos de ASR pueden extraer ancestros que probablemente tienen miles de millones de años, a menudo hay decenas, si no cientos, de cambios de secuencia entre los propios ancestros y los ancestros y las secuencias existentes; por ello, estos estudios evolutivos de secuencia-función pueden requerir mucho trabajo y una dirección racional. [ 1 ] [ 13 ] [ 14 ]

Proteínas resucitadas

Existen numerosos ejemplos de proteínas ancestrales que han sido reconstruidas computacionalmente, expresadas en líneas celulares vivas y, en muchos casos, purificadas y estudiadas bioquímicamente.

Otros ejemplos son los pigmentos visuales ancestrales en vertebrados, [ 19 ] enzimas en levaduras que descomponen azúcares (800 Ma); [ 20 ] enzimas en bacterias que proporcionan resistencia a los antibióticos (2 – 3 Ga ); [ 21 ] las ribonucleasas involucradas en la digestión de rumiantes; las alcohol deshidrogenasas (Adhs) involucradas en la fermentación de levaduras (~85 Ma); [ 12 ] y RuBisCO en Solanáceas . [ 22 ]

La "edad" de una secuencia reconstruida se determina utilizando un modelo de reloj molecular , y a menudo se emplean varios. [ 9 ] [ 23 ] Esta técnica de datación suele calibrarse utilizando puntos de tiempo geológicos (como constituyentes oceánicos antiguos o BIF ) y, si bien estos relojes ofrecen el único método para inferir la edad de una proteína muy antigua , tienen amplios márgenes de error y son difíciles de defender contra datos contrarios. Por ello, la "edad" de ASR debería utilizarse solo como una característica indicativa y a menudo se supera por completo para una medición del número de sustituciones entre las secuencias ancestrales y modernas (el fundamento sobre el que se calcula el reloj). [ 6 ] Dicho esto, el uso de un reloj permite comparar los datos biofísicos observados de una proteína ASR con el entorno geológico o ecológico de la época. Por ejemplo, los estudios ASR sobre EF-Tus bacterianos (proteínas involucradas en la traducción , que probablemente rara vez están sujetas a HGT y que suelen exhibir Tms ~2C mayor que Tenv) indican una Tierra precámbrica más caliente que coincide muy estrechamente con los datos geológicos sobre las temperaturas oceánicas de la Tierra antigua basados ​​en los niveles isotópicos de oxígeno-18 . [ 11 ] Los estudios ASR de Adhs de levadura revelan que la aparición de Adhs subfuncionalizadas para el metabolismo del etanol (no solo la excreción de desechos) surgió en un momento similar al inicio de la fruta carnosa en el período Cámbrico y que antes de esta aparición, Adh servía para excretar etanol como subproducto del exceso de piruvato . [ 12 ] El uso de un reloj también podría indicar que el origen de la vida ocurrió antes de lo que indican los fósiles moleculares más antiguos (>4.1 Ga), pero dada la discutible fiabilidad de los relojes moleculares, tales observaciones deben tomarse con precaución. [ 23 ] [ 24 ]

Tiorredoxina

Un ejemplo es la reconstrucción de enzimas tiorredoxina de organismos de hasta 4 mil millones de años de antigüedad. [ 25 ] Si bien la actividad química de estas enzimas reconstruidas era notablemente similar a la de las enzimas modernas, sus propiedades físicas mostraron una estabilidad térmica y ácida significativamente elevada. Estos resultados se interpretaron como una indicación de que la vida antigua pudo haber evolucionado en océanos mucho más cálidos y ácidos que los actuales. [ 25 ]

Significado

Estos experimentos abordan varias preguntas importantes en biología evolutiva: ¿ la evolución procede en pequeños pasos o en grandes saltos?; ¿es reversible la evolución?; ¿cómo evoluciona la complejidad ? Se ha demostrado que pequeñas mutaciones en la secuencia de aminoácidos de los receptores hormonales determinan un cambio importante en sus preferencias por las hormonas. Estos cambios implican grandes avances en la evolución del sistema endocrino . Por lo tanto, cambios muy pequeños a nivel molecular pueden tener consecuencias enormes. El laboratorio de Thornton también ha podido demostrar que la evolución es irreversible estudiando el receptor de glucocorticoides . Este receptor fue modificado por siete mutaciones en un receptor de cortisol, pero revertir estas mutaciones no devolvió el receptor original. Esto indica que la epistasis juega un papel importante en la evolución de las proteínas, una observación que, en combinación con las observaciones de varios ejemplos de evolución paralela, apoya el modelo de red neutral mencionado anteriormente. [ 7 ] Otras mutaciones neutrales anteriores actuaron como un trinquete e hicieron que los cambios en el receptor fueran irreversibles. [ 26 ] Estos diferentes experimentos con receptores muestran que, durante su evolución, las proteínas se diferencian enormemente y esto explica cómo puede evolucionar la complejidad. Un análisis más detallado de los distintos receptores hormonales ancestrales y de las diversas hormonas revela que, a nivel de interacción entre residuos de aminoácidos individuales y grupos químicos hormonales, se producen cambios muy pequeños pero específicos. El conocimiento de estos cambios podría, por ejemplo, conducir a la síntesis de equivalentes hormonales capaces de imitar o inhibir la acción de una hormona, lo que abriría nuevas posibilidades terapéuticas.

Dado que la ASR ha revelado una tendencia hacia la termoestabilidad antigua y la promiscuidad enzimática, se presenta como una herramienta valiosa para los ingenieros de proteínas que a menudo buscan estas características (produciendo efectos a veces mayores que las herramientas actuales, dirigidas racionalmente). [ 10 ] La ASR también promete «resucitar» «organismos antiguos» fenotípicamente similares, lo que a su vez permitiría a los bioquímicos evolutivos investigar la historia de la vida. Los defensores de la ASR, como Benner, afirman que, a través de estos y otros experimentos, a finales del siglo actual se alcanzará un nivel de comprensión en biología análogo al que surgió en la química clásica en el siglo pasado. [ 12 ]

Referencias

  1. 1 2 Thornton JW (mayo de 2004). "Resucitando genes antiguos: análisis experimental de moléculas extintas". Nature Reviews. Genetics . 5 (5): 366– 375. doi : 10.1038/nrg1324 . PMID 15143319 . S2CID 205482979 .  
  2. Pauling L, Zuckerkandl E (1963). "Paleogenética química: estudios de restauración molecular de formas de vida extintas" . Acta Chemica Scandinavica . 17 : S9– S16.]
  3. Jermann TM, Opitz JG, Stackhouse J, Benner SA (marzo de 1995). "Reconstrucción de la historia evolutiva de la superfamilia de la ribonucleasa artiodáctila". Nature . 374 ( 6517): 57– 59. Bibcode : 1995Natur.374...57J . doi : 10.1038/374057a0 . PMID 7532788. S2CID 4315312 .  
  4. 1 2 Thornton JW, Need E, Crews D (septiembre de 2003). "Resucitando el receptor de esteroides ancestral: origen antiguo de la señalización de estrógenos". Science . 301 (5640): 1714– 1717. Bibcode : 2003Sci...301.1714T . doi : 10.1126/science.1086185 . PMID 14500980 . S2CID 37628350 .  
  5. Pearson H (marzo de 2012). "Proteínas prehistóricas: resucitando a los muertos" . Nature . 483 (7390). Londres: 390–3 . Bibcode : 2012Natur.483..390P . doi : 10.1038/483390a . PMID 22437590 . 
  6. 1 2 3 4 Lim SA, Hart KM, Harms MJ, Marqusee S (noviembre de 2016). "Tendencia evolutiva hacia la estabilidad cinética en la trayectoria de plegamiento de las RNasas H" . Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 113 (46): 13045– 13050. Bibcode : 2016PNAS..11313045L . doi : 10.1073/pnas.1611781113 . PMC 5135364. PMID 27799545 .  
  7. 1 2 3 Harms MJ, Thornton JW (agosto de 2013). "Bioquímica evolutiva: revelando las causas históricas y físicas de las propiedades de las proteínas" . Nature Reviews. Genetics . 14 (8): 559– 571. doi : 10.1038/nrg3540 . PMC 4418793. PMID 23864121 .  
  8. Hobbs JK, Prentice EJ, Groussin M, Arcus VL (octubre de 2015). "Las enzimas ancestrales reconstruidas imponen un costo de aptitud a las bacterias modernas a pesar de exhibir propiedades bioquímicas favorables". Journal of Molecular Evolution . 81 ( 3–4 ): 110–120 . Bibcode : 2015JMolE..81..110H . doi : 10.1007/s00239-015-9697-5 . hdl : 1721.1 /105120 . PMID 26349578. S2CID 18833850 .  
  9. 1 2 Wheeler LC, Lim SA, Marqusee S, Harms MJ (junio de 2016). " La termoestabilidad y especificidad de las proteínas antiguas" . Current Opinion in Structural Biology . 38 : 37–43 . doi : 10.1016/j.sbi.2016.05.015 . PMC 5010474. PMID 27288744 .  
  10. 1 2 Risso VA, Gavira JA, Sanchez-Ruiz JM (junio de 2014). "Proteínas precámbricas termoestables y promiscuas". Microbiología ambiental . 16 (6): 1485– 1489. Bibcode : 2014EnvMi..16.1485R . doi : 10.1111/1462-2920.12319 . hdl : 10481/87187 . PMID 25009840 . 
  11. 1 2 Gaucher EA, Govindarajan S, Ganesh OK (febrero de 2008). "Tendencia de paleotemperatura para la vida precámbrica inferida a partir de proteínas resucitadas". Nature . 451 (7179): 704– 707. Bibcode : 2008Natur.451..704G . doi : 10.1038/nature06510 . PMID 18256669 . S2CID 4311053 .  
  12. 1 2 3 4 Liberles DA, ed. (26 de julio de 2007). Reconstrucción de la secuencia ancestral . Oxford, Nueva York: Oxford University Press. ISBN 978-0-19-929918-8.
  13. Anderson DP, Whitney DS, Hanson-Smith V, Woznica A, Campodonico-Burnett W, Volkman BF, et al. (enero de 2016). "Evolución de una antigua función proteica involucrada en la multicelularidad organizada en animales" . eLife . 5 e10147. doi : 10.7554/eLife.10147 . PMC 4718807. PMID 26740169 .   
  14. Figura 1 de la referencia Harms MJ, Thornton JW (agosto de 2013). "Bioquímica evolutiva: revelando las causas históricas y físicas de las propiedades de las proteínas" . Nature Reviews. Genetics . 14 (8): 559– 571. doi : 10.1038/nrg3540 . PMC 4418793. PMID 23864121 .  
  15. Eick GN, Colucci JK, Harms MJ, Ortlund EA, Thornton JW (2012). "Evolución de la especificidad mínima y la promiscuidad en los receptores de hormonas esteroides" . PLOS Genetics . 8 (11) e1003072. doi : 10.1371/journal.pgen.1003072 . PMC 3499368. PMID 23166518 .  
  16. Harms MJ, Eick GN, Goswami D, Colucci JK, Griffin PR, Ortlund EA, Thornton JW (julio de 2013). "Mecanismos biofísicos para mutaciones de gran efecto en la evolución de los receptores de hormonas esteroides" . Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 110 (28): 11475– 80. Bibcode : 2013PNAS..11011475H . doi : 10.1073/pnas.1303930110 . PMC 3710831. PMID 23798447 .  
  17. Finnigan GC, Hanson-Smith V, Stevens TH, Thornton JW (enero de 2012). "Evolución de la complejidad aumentada en una máquina molecular" . Nature . 481 ( 7381): 360–364 . Bibcode : 2012Natur.481..360F . doi : 10.1038/nature10724 . PMC 3979732. PMID 22230956 .  
  18. Hart KM, Harms MJ, Schmidt BH, Elya C, Thornton JW, Marqusee S (noviembre de 2014). "Deriva del sistema termodinámico en la evolución de las proteínas" . PLOS Biology . 12 (11) e1001994. doi : 10.1371/journal.pbio.1001994 . PMC 4227636. PMID 25386647 .  
  19. Shi Y, Yokoyama S (julio de 2003). "Análisis molecular de la importancia evolutiva de la visión ultravioleta en vertebrados" . Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 100 (14): 8308– 8313. Bibcode : 2003PNAS..100.8308S . doi : 10.1073/pnas.1532535100 . PMC 166225. PMID 12824471 .  
  20. Voordeckers K, Brown CA, Vanneste K, van der Zande E, Voet A, Maere S, Verstrepen KJ (2012). "Reconstrucción de enzimas metabólicas ancestrales revela mecanismos moleculares subyacentes a la innovación evolutiva a través de la duplicación genética" . PLOS Biology . 10 (12) e1001446. doi : 10.1371/journal.pbio.1001446 . PMC 3519909. PMID 23239941 .  
  21. Risso VA, Gavira JA, Mejia-Carmona DF, Gaucher EA, Sanchez-Ruiz JM (febrero de 2013). "Hiperestabilidad y promiscuidad de sustrato en resurrecciones de laboratorio de β-lactamasas precámbricas". Journal of the American Chemical Society . 135 (8): 2899– 2902. Bibcode : 2013JAChS.135.2899R . doi : 10.1021/ja311630a . hdl : 11336/22624 . PMID 23394108. S2CID 207092445 .  
  22. Lin, Myat T.; Salihovic, Heidi; Clark, Frances K.; Hanson, Maureen R. (15 de abril de 2022). "Mejorando la eficiencia de Rubisco resucitando a sus ancestros en la familia Solanaceae" . Science Advances . 8 (15) eabm6871. Bibcode : 2022SciA....8M6871L . doi : 10.1126/sciadv.abm6871 . PMC 9012466. PMID 35427154 .  
  23. 1 2 Battistuzzi FU, Feijao A, Hedges SB (noviembre de 2004). "Una escala de tiempo genómica de la evolución de los procariotas: perspectivas sobre el origen de la metanogénesis, la fototrofia y la colonización de la tierra" . BMC Evolutionary Biology . 4 44. doi : 10.1186/1471-2148-4-44 . PMC 533871. PMID 15535883 .  
  24. Bell EA , Boehnke P, Harrison TM, Mao WL (noviembre de 2015). "Carbono potencialmente biogénico preservado en un circón de 4.100 millones de años" . Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 112 (47): 14518–14521 . Bibcode : 2015PNAS..11214518B . doi : 10.1073/pnas.1517557112 . PMC 4664351. PMID 26483481 .  
  25. 1 2 Pérez-Jiménez R, Inglés-Prieto A, Zhao ZM, Sánchez-Romero I, Alegre-Cebollada J, Kosuri P, et al. (mayo de 2011). " La paleoenzimología de molécula única explora la química de las enzimas resucitadas" . Nature Structural & Molecular Biology . 18 (5): 592– 596. doi : 10.1038/nsmb.2020 . PMC 3087858. PMID 21460845 .   
  26. Bridgham JT, Ortlund EA, Thornton JW (septiembre de 2009). "Un trinquete epistático restringe la dirección de la evolución del receptor de glucocorticoides" . Nature . 461 ( 7263): 515– 519. Bibcode : 2009Natur.461..515B . doi : 10.1038/nature08249 . PMC 6141187. PMID 19779450 .