Genes Hox
Genes Hox
Genes Hox
Facultat de Ciències
ÍNDICE
ÍNDICE DE TABLAS iii
ÍNDICE DE FIGURAS v
ABREVIATURAS vii
RESUMEN ix
ABSTRACT xi
1. INTRODUCCIÓN 1
1.1. LA UTILIDAD DE LOS ORGANISMOS MODELO 4
1.2. DROSOPHILA COMO MODELO 8
1.3. DROSOPHILA COMO MODELO EN BIOLOGÍA DEL DESARROLLO 12
1.3.1. Ciclo Vital 12
1.3.1.1. Embriogénesis 13
1.3.1.2. Metamorfosis 14
1.3.2. Genética del Desarrollo Embrionario 15
1.4. DROSOPHILA COMO MODELO EN GENÓMICA 17
1.4.1. El Número de Genes 18
1.4.2. Las Secuencias Reguladoras 19
1.4.3. Genómica Comparada 21
1.5. EVOLUCIÓN CROMOSÓMICA 22
1.5.1. Origen de las Reordenaciones Cromosómicas 24
1.5.2. Importancia Evolutiva de los Segmentos Conservados 24
1.6. EL COMPLEJO DE GENES HOX 26
1.6.1. Los Genes Homeóticos 26
1.6.2. Agrupación y Colinealidad 28
1.6.3. Origen del Complejo Hox 29
1.6.4. Los Genes Hox y la Evolución Morfológica 31
1.6.5. Estructura del Complejo en los Metazoos 35
1.7. LOS COMPLEJO DE GENES HOX EN DROSOPHILA 37
1.7.1. Estructura del Complejo 37
1.7.2. Activación y Mantenimiento de los Patrones de Expresión Hox 39
1.7.3. Función de los Genes Hox 40
1.7.4. Cambios Estructurales del HOM-C en el Género Drosophila 41
1.8. OBJETIVOS 44
1.8. OBJECTIVES 45
i
INDICE
2. RESULTADOS 47
2.1. ARTÍCULO 1 49
2.2. ARTÍCULO 2 65
3. DISCUSIÓN 115
3.1. EVOLUCIÓN DE LAS REGIONES CIS-REGULADORAS 117
3.2. EL HOM-C EN DROSOPHILA 119
3.2.1. Evolución de los Genes Hox: Estructura, Regulación y Función 119
3.2.1.1. Hox3: zen, bcd y zen2 119
3.2.1.2. lab 122
3.2.1.3. pb 124
3.2.1.4. abdA 127
3.2.2. Evolución de las Estructura del HOM-C en Drosophila 128
3.2.2.1. Estructura del HOM-C 128
3.2.2.2. Reconstrucción de la Evolución de la Estructura del
HOM-C 130
3.2.3. Mecanismo de las Reordenaciones 132
3.2.4. Causas y Consecuencias de las Roturas del HOM-C 133
3.2.5. ANT-C versus BX-C 136
3.3. CONSERVACIÓN Y DESINTEGRACIÓN DEL HOM-C EN LOS METAZOOS 138
3.3.1. Drosophila versus Vertebrados 138
3.3.2. Estructura del HOM-C en los Metazoos 139
3.3.3. ¿Qué Tienen en Común los Organismos con el Complejo
Fragmentado? 142
4. CONCLUSIONES 145
4. CONCLUSIONS 147
5. REFERENCIAS 149
ii
ÍNDICE DE TABLAS
ÍNDICE DE TABLAS
INTRODUCCIÓN
Tabla 1. Homologías cromosómicas entre algunas especies del género
Drosophila. 23
RESULTADOS
Artículo 1 – Tabla 1. Secuencia de los cebadores utilizados para la
amplificación por PCR de diferentes regiones del gen lab en varias
especies de Drosophila. 53
Artículo 1 – Tabla 2. Estructura del gen lab en tres especies de Drosophila
y Tribolium castaneum. 55
Artículo 1 – Tabla 3. Estimas por máxima verosimilitud del número de
substituciones sinónimas (dS) y no sinónimas (dN) por posición en
regiones codificantes del gen lab en comparaciones a pares entre D.
buzzatii, D. virilis y D. melanogaster. 58
Artículo 1 – Tabla 4. Test de la razón de máxima verosimilitud de la
variación en la evolución molecular y razón ω del gen lab en el género
Drosophila. 59
Artículo 1 – Tabla 5. Logaritmos de los valores de verosimilitud y estimas
de los parámetros bajo modelos de sitios fijados. 59
Artículo 2 – Tabla 1. Características de las regiones conservadas no
codificantes (CNS) detectadas con mVISTA en comparaciones de las
regiones de los genes Hox entre D. melanogaster, D. pseudoobscura y
D. buzzatii. 70
Artículo 2 – Tabla 2. Características de las regiones conservadas no
codificantes (CNS) detectadas con mVISTA en comparaciones de las
regiones de los genes no Hox entre D. melanogaster, D. pseudoobscura
y D. buzzatii. 71
Artículo 2 – Tabla S1. Análisis estadístico de las características de los
CNS entre las tres comparaciones a pares entre especies. 77
Artículo 2 – Tabla S2. Análisis estadístico de las características de los
CNS entre parejas de especies para las diferentes regiones. 77
iii
ÍNDICE DE TABLAS
iv
ÍNDICE DE FIGURAS
ÍNDICE DE FIGURAS
INTRODUCCIÓN
Figura 1. Filogenia de los metazoos. 7
Figura 2. Filogenia de los principales grupos de artrópodos e insectos. 7
Figura 3. Árbol filogenético de los principales grupos de Drosophila. 11
Figura 4. Ciclo vital de Drosophila. 12
Figura 5. Fases del desarrollo embrionario de Drosophila. 14
Figura 6. Discos y tejidos imaginales de Drosophila. 15
Figura 7. Cascada de activación de los genes de segmentación. 17
Figura 8. Transformaciones homeóticas en Drosophila. 26
Figura 9. Conservación de la organización genómica y patrones de
expresión de los genes Hox. 30
Figura 10. Distribución de los genes Hox en los Bilateria. 32
Figura 11. Jerarquía de cambios evolutivos en los genes Hox según su
presunto efecto morfológico. 34
Figura 12. Estructura del HOM-C en D. melanogaster. 38
Figura 13. Genes que intervienen en la regulación del patrón de expresión
de los genes Hox. 40
Figura 14. Reconstrucción filogenética de la organización de los genes
Hox en el género Drosophila. 42
RESULTADOS
Artículo 1 – Figura 1. Reconstrucción filogenética de la organización de
los genes Hox en el género Drosophila. 52
Artículo 1 – Figura 2. Localización cromosómica de los genes lab y abdA
en el grupo repleta. 54
Artículo 1 – Figura 3. Mapa de restricción y organización molecular de la
región lab en D. buzzatii. 56
Artículo 1 – Figura 4. Alineamiento de la 5’ UTR y secuencias 5’ de lab
en D. buzzatii, D. virilis y D. melanogaster. 57
Artículo 1 – Figura 5. Gráfico VISTA mostrando la similaridad en las
secuencias nucleotídicas del gen lab entre D. melanogaster y D.
buzzatii. 58
v
ÍNDICE DE FIGURAS
DISCUSIÓN
Figura 15. Modelo de evolución del gen Hox3. 122
Figura 16. Expresión del gen pb en embriones de cinco especies de
Drosophila. 126
Figura 17. Modelo de evolución del HOM-C en Drosophila. 131
Figura 18. Estructura del HOM-C en los Metazoos. 140
vi
ABREVIATURAS
ABREVIATURAS
vii
ABREVIATURAS
viii
RESUMEN
RESUMEN
Los genes homeóticos (Hox) codifican factores de transcripción involucrados en
la especificación de la identidad segmental a lo largo del eje anteroposterior en el
embrión de los metazoos. Estos genes se presentan habitualmente agrupados y
organizados en el mismo orden en el que se expresan a lo largo del eje anteroposterior
del cuerpo. La conservación de esta organización genómica a lo largo de la filogenia ha
sugerido la existencia de constricciones funcionales actuando sobre ella, pero la
desestructuración del complejo en Caenorhabditis elegans y las tres roturas descritas en
el género Drosophila ponen en cuestión que esta organización sea realmente necesaria
para su correcto funcionamiento en algunos linajes. En este trabajo se han estudiado las
consecuencias genómicas y funcionales de las dos reorganizaciones del complejo de
genes Hox presentes en Drosophila buzzatii, especie perteneciente al grupo repleta. En
la primera parte del trabajo se han clonado y secuenciado las regiones codificantes del
gen labial (lab) en D. buzzatii y D. virilis. Se han comparado las secuencias de estas dos
especies y la de D. melanogaster para comprobar si el cambio de posición del gen lab
ocurrido en el linaje de D. buzzatii había producido algún cambio en la estructura del
gen o en la evolución de su secuencia. Los resultados muestran una tasa de cambio
heterogénea a lo largo del gen pero homogénea a lo largo de la filogenia. La tasa de
cambio nucleotídico de lab se ha mantenido constante a pesar del cambio de posición.
En la segunda parte del trabajo se han secuenciado dos regiones del genoma de D.
buzzatii, una contiene los genes lab y abdominal A (abdA) y la otra incluye el gen
proboscipedia (pb), y se ha comparado su organización génica con D. melanogaster y
D. pseudoobscura para localizar con precisión los puntos de rotura. También se han
comparado las secuencias no codificantes de estas regiones, se ha observado una
elevada presencia de bloques conservados en los intrones y regiones adyacentes de los
genes Hox. La comparación de los bloques conservados con las regiones reguladoras
conocidas de los genes Hox (lab, pb y abdA) en D. melanogaster muestra que la
posición y orientación de las regiones reguladoras está conservada entre las tres
especies, con excepciones menores. Finalmente se analizó el patrón de expresión de los
tres genes Hox en embriones y discos imaginales de D. buzzatii, D. repleta, D. virilis, y
D. melanogaster, cuatro especies con diferentes organizaciones de los genes Hox. Las
dos roturas estudiadas se produjeron mediante inversiones paracéntricas, con los puntos
de rotura respetando las regiones reguladoras de ambos genes. De manera que las
regiones reguladoras y patrones de expresión de los genes Hox adyacentes se han
ix
RESUMEN
x
ABSTRACT
ABSTRACT
Homeotic (Hox) genes code for transcription factors involved in the
specification of segmental identity in the anteroposterior axis of the early metazoan
embryo. These genes are usually clustered and arranged in the same order as they are
expressed along the anteroposterior body axis. The conservation of this Hox gene
organization along the phylogeny has suggested the existence of functional constraints.
However, the partial disassembly of the Caenorhabditis elegans complex and the three
splits observed in the Drosophila genus question whether this organization is an
absolute necessity for proper function in some lineages. In this work, I analysed the
genomic and functional consequences of the two splits present in Drosophila buzzatii, a
member of the repleta species group. In the first part, the coding regions of the labial
(lab) gene were cloned in D. buzzatii and D. virilis. The sequences of these two species
were compared with that of D. melanogaster to test whether the change in position of
lab in de D. buzzatii lineage produced any change in gene structure or sequence
evolution. The results show that the substitution rate is heterogeneous along the gene
but homogeneous along the phylogeny. The nucleotide substitution rate of lab has been
constant in spite of the positional change. In the second part, two regions of the D.
buzzatii genome have been sequenced, one including the lab y abdominal A (abdA)
genes and the other containing the proboscipedia (pb) gene, and compared with the
genic organization of D. melanogaster and D. pseudoobscura to precisely locate the
breakpoints. The noncoding sequences of these regions have also been compared, and a
high presence of conserved blocks has been observed in the introns and surrounding
regions of Hox genes. The comparison of these conserved blocks with the known
regulatory regions of the Hox genes (lab, pb and abdA) in D. melanogaster shows that
the position and order of the regulatory regions is conserved between the three species,
with minor exceptions. Finally the expression pattern of the three Hox genes has been
analysed in embryos and imaginal discs of D. buzzatii, D. repleta, D. virilis, and D.
melanogaster, four species with different Hox gene arrangements. The two splits took
place through two paracentric inversions, with their breakpoints between the regulatory
regions of adjacent genes. So that the regulatory regions and expression patterns of
these Hox genes have been conserved in spite of the reorganizations. In Drosophila the
Hox gene complex seems to be composed by independent modules (including the gene
and its regulatory regions), whose association is not required for proper function. The
organization of these genes is modular and their clustering seems de result of
xi
ABSTRACT
xii
1. INTRODUCCIÓN
INTRODUCIÓN
3
INTRODUCCIÓN
4
INTRODUCIÓN
5
INTRODUCCIÓN
Para poder aplicar los conocimientos que se van adquiriendo en una especie
modelo a otros organismos o hacer comparaciones entre especies modelo es importante
tener una perspectiva evolutiva. Los estudios evolutivos requieren una buena base
filogenética para poder interpretar los cambios detectados. La aplicación de métodos
moleculares en la elaboración de filogenias en las últimas dos décadas a cambiado
notablemente la filogenia de los metazoos. Este hecho ha influido en la interpretación de
algunas comparaciones. Especialmente notable es el cambio de posición sufrido por los
nemátodos que han pasado de ser considerados un grupo primitivo, rama basal a
protóstomos y deuteróstomos, a ser un grupo hermano de los artrópodos (AGUINALDO et
al. 1997).
Para referencia a lo largo del presente trabajo, se incluyen aquí dos árboles
filogenéticos. El primero es una filogenia de los grandes grupos de metazoos (Figura 1).
El segundo sitúa al género Drosophila dentro de los artrópodos, indicando también la
posición de otros géneros mencionados en el presente trabajo (Figura 2).
6
INTRODUCIÓN
Figura 1. Filogenia de los Metazoos. Modificado de AGUINALDO et al. (1997), DE ROSA et al.
(1999) MARTINEZ et al. (1999) y DANCHIN y PONTAROTTI (2004).
Figure 1. Phylogeny of the Metazoan. Modified from AGUINALDO et al. (1997), DE ROSA et al.
(1999) MARTINEZ et al. (1999) and DANCHIN and PONTAROTTI (2004).
7
INTRODUCCIÓN
8
INTRODUCIÓN
9
INTRODUCCIÓN
como el sistema Flp-FRP (XU y RUBIN 1993; GOLIC y GOLIC 1996), que permiten la
generación de inversiones o la delección de regiones determinadas. En los últimos años
han proliferado los usos de las técnicas de transgénesis gracias a ingenios y variaciones
sobre estas técnicas, y se están desarrollando también nuevos vectores para evitar las
limitaciones de los sistemas actuales como son el tamaño de las secuencias que se
pueden insertar o los efectos de posición (VENKEN y BELLEN 2005). Estas herramientas
de manipulación genética se han beneficiado del mantenimiento sencillo y económico
de Drosophila, que permite mantener de forma rutinaria gran cantidad de cepas
“salvajes” o mutantes, tanto espontáneas como inducidas.
10
INTRODUCIÓN
11
INTRODUCCIÓN
indican la existencia de dos linajes principales dentro del género Drosophila, los cuales
divergieron hace entre 40 y 62 millones de años (Figura 3) (POWELL y DESALLE 1995;
RUSSO et al. 1995; TAMURA et al. 2004). Uno de los linajes dió lugar al subgénero
Sophophora el cual contiene unas 300 especies, entre ellas se encuentra D.
melanogaster. El segundo linaje dió lugar a los subgéneros Drosophila y Idiomyia
(Drosophilas Hawaianas), que incluyen unas 700 y 375 especies descritas
respectivamente (POWELL 1997) (http://taxodros.unizh.ch).
12
INTRODUCIÓN
13
INTRODUCCIÓN
Figura 5. Fases del desarrollo embrionario de Drosophila. Los embriones se muestran con la
región anterior en la izquierda y la dorsal arriba. El número sobre cada embrión indica el estadío
que representa. Tomado de HARTENSTEIN (1993).
Figure 5. Phases of embryonic development in Drosophila. The embryos are shown with
anterior to the left and posterior up. The number over each embryo corresponds to the stage.
Taken from HARTENSTEIN (1993).
14
INTRODUCIÓN
15
INTRODUCCIÓN
16
INTRODUCIÓN
polaridad segmental son de este tipo. Esta diferencia se debe principalmente a que
cuando estos genes se expresan las células ya se han formado, y por tanto es necesario
activar, entre otros, los mecanismos de comunicación celular para enviar y recibir
señales entre células (NASIADKA et al. 2002) lo que permitirá el desarrollo coordinado
de las diferentes partes.
Los datos genómicos sin ordenar son simplemente una retahíla de As, Ts, Gs y
Cs. En este estado, los datos de secuencia aportan muy poco conocimiento biológico.
Para entender el legado genético de un organismo debemos interpretar su secuencia
genómica, traduciendo la información que contiene en forma molecular en anotaciones
17
INTRODUCCIÓN
leíbles por el hombre (LEWIS et al. 2002). El proceso de anotación es la tarea de anclar
en la secuencia de ácidos nucleicos notaciones explicativas de sus características o
función: transcrito, exón, promotor, elemento transponible, región reguladora, etc. La
anotación funcional de las secuencias genómicas de los eucariotas es uno de los
principales retos de la biología moderna.
El pequeño tamaño del genoma de Drosophila con 176 Megabases (Mb)
(HOSKINS et al. 2002), la ha convertido en un modelo fundamental en Genómica,
disciplina centrada en el estudio de la secuencia, estructura y función del genoma.
Aprender cómo identificar, exponer, interrogar e interpretar las características del
genoma en un organismo bien caracterizado como Drosophila, es crucial para entender
los genomas de organismos más complejos, incluyendo humanos. Añadir anotaciones a
la secuencia de un genoma de forma rigurosa y consistente es un prerrequisito para el
uso eficiente de esa secuencia en la investigación biológica (MISRA et al. 2002).
La secuencia anotada del genoma de D. melanogaster es resultado de la primera
aplicación del método “shotgun” para la secuenciación de un genoma animal (ADAMS et
al. 2000), y supone un modelo para la anotación a gran escala de las secuencias
genómicas. En los últimos años se han desarrollado multitud de programas informáticos
encaminados a identificar los genes y las regiones cis-reguladoras que controlan su
expresión, con el objetivo final de automatizar el proceso de anotación de los genomas.
18
INTRODUCIÓN
ligeramente, 13.379 genes, a pesar de presentar una secuencia más completa e incluir la
heterocromatina, ausente en el ensamblaje inicial (MISRA et al. 2002). Aunque entre el
Release 1 y el Release 3 el número de genes anotados ha variado poco, sí ha variado la
estructura de muchos de ellos. Con una predicción menos restrictiva usando Fgenesh,
análisis de la expresión por “microarrays” y validación por RT-PCR e hibridación in
situ, HILD et al (2003) identifican al menos 2,000 nuevos genes en el genoma de D.
melanogaster no incluidos en la anotación del Release 3. Aunque se ha liberado un
Release 4 de la eucromatina, donde se han cerrado 21 gaps de secuencia, la anotación
no ha variado respecto al Release 3 (http://www.fruitfly.org/annot/release4.html).
19
INTRODUCCIÓN
20
INTRODUCIÓN
gran limitación para la aplicación generalizada de este método. Como los CRM están
sujetos a constricciones funcionales para mantener la expresión del gen diana su
secuencia evoluciona más despacio que las secuencias no funcionales. Por lo tanto, la
conservación de secuencias no codificantes entre especies filogenéticamente distantes
(huella filogenética) puede utilizarse como guía para identificar secuencias reguladoras
(TAGLE et al. 1988). Varios estudios recientes avalan el uso del análisis comparativo de
secuencias y la caracterización de las secuencias no codificantes conservadas (CNS –
Conserved Noncoding Sequences) como una aproximación útil para la detección de
posibles CRM en Drosophila (BERGMAN y KREITMAN 2001; BERGMAN et al. 2002),
mamíferos (DERMITZAKIS et al. 2002; BOFFELLI et al. 2003; COOPER y SIDOW 2003;
NOBREGA et al. 2003; BOFFELLI et al. 2004; FRAZER et al. 2004b; GAFFNEY y
KEIGHTLEY 2004; MARTIN et al. 2004), vertebrados (SANTINI et al. 2003) y plantas
(GUO y MOOSE 2003; BUCHANAN et al. 2004). En Drosophila, la interpretación de los
CNS como probables regiones reguladoras está respaldada por la elevada tasa de
substitución nucleotídica neutra (MORIYAMA y GOJOBORI 1992) y la elevada tasa
intrínseca de pérdida de DNA (PETROV et al. 1996; SINGH y PETROV 2004). No se
espera que las secuencia no codificantes estén conservadas entre especies poco
relacionadas a no ser que estén sometidas a constricciones funcionales (O'NEIL y
BELOTE 1992).
21
INTRODUCCIÓN
En las especies del género Drosophila el cariotipo es muy variable por lo que refiere al
número haploide de cromosomas que está habitualmente entre 3 y 6. Gracias a la
información de mapas citológicos y de ligamiento, MULLER (1940) formuló la hipótesis
de que el cariotipo ancestral del género estaba formado por seis elementos
cromosómicos que se han conservado como cromosomas o brazos cromosómicos en
todo el género. Estableció homologías entre las diferentes especies de Drosophila
elaborando una nomenclatura que permite relacionar los cariotipos. Así los elementos
22
INTRODUCIÓN
23
INTRODUCCIÓN
24
INTRODUCIÓN
SKARMETA et al. 1996) o knirps y knirps related (LUNDE et al. 1998) están conservados
en D. repleta (GONZÁLEZ et al. 2002) lo que sugiere la acción de la selección natural al
menos en algunos casos.
Estudios recientes parecen indicar que la organización lineal de los genes en los
cromosomas no es aleatoria (KOSAK y GROUDINE 2004). Varios estudios mediante
microarrays o “análisis serial de la expresión génica” (SAGE – Serial Analysis of Gene
Expression) muestran agrupación genómica de genes coexpresados en Drosophila
(BOUTANAEV et al. 2002; SPELLMAN y RUBIN 2002), nemátodos, mamíferos (LERCHER
et al. 2002) o Arabidopsis (WILLIAMS y BOWLES 2004). Se ha detectado la agrupación
de genes “housekeeping” con nivel de transcripción elevado (LERCHER et al. 2002), de
genes con patrones de expresión similares (SPELLMAN y RUBIN 2002) o de genes de
tejidos específicos (BOUTANAEV et al. 2002). ¿Podría la coexpresión ser un factor en el
mantenimiento de la sintenia? Si la expresión coordinada de estos genes es ventajosa,
deberíamos encontrar en las regiones coexpresadas más sintenia de la esperada
(SPELLMAN y RUBIN 2002).
Las presiones selectivas que promueven la agrupación no están relacionadas
simplemente con la tasa de transcripción, sino que favorecen la unión de los genes que
son activos en un tejido determinado o en todos los tejidos (LERCHER et al. 2002). La
corregulación transcripcional en eucariotas superiores opera en extensos dominios
cromatínicos que comprenden múltiples genes (BOUTANAEV et al. 2002). La estructura
de la cromatina, y por tanto, la accesibilidad de diferentes regiones genómicas a la
maquinaria transcripcional, puede variar según el tipo celular. En esta situación, podría
ser ventajoso agrupar los genes “housekeeping” en una región que mantenga una
conformación accesible en cualquier tipo celular o los genes específicos de testículo en
una región que se active sólo en esta fase. Aunque todavía falta verificación
experimental, la regulación transcripcional coordinada al nivel de la cromatina podría
ser un sistema económico de remodelar las estructuras cromosómica (BOUTANAEV et al.
2002; LERCHER et al. 2002; UEDA et al. 2002). Por otro lado, la agrupación de genes
coexpresados podría también ser el resultado de la expresión residual de genes próximos
a los verdaderamente activos debido a la descompactación de la cromatina de la región.
25
INTRODUCCIÓN
26
INTRODUCIÓN
Desde principios de los 40, Lewis analizó las bases genéticas de las
transformaciones homeóticas, en un intento de contrastar la hipótesis de Bridges, basada
en evidencias citológicas, de que el genoma contiene duplicaciones génicas naturales, a
menudo en tandem. Empezó trabajando con diferentes mutantes ya descritos que
causaban transformaciones similares a bx y que, en muchos casos, complementaban a
esta mutación. Estas mutaciones tenían la particularidad de estar agrupadas en una
región del cromosoma, que fue bautizada como Complejo Bithorax (BX-C). Los genes
responsables de este fenómeno resultaron ser miembros de una familia génica que
controla la identidad de los segmentos a lo largo del eje anteroposterior del cuerpo.
Lewis trabajó en estos problemas durante décadas y en 1978 publicó un resumen de sus
resultados. En esta revisión formuló teorías sobre cómo interactúan los genes
homeóticos, cómo el orden de estos genes se corresponde al orden de los segmentos a lo
largo del eje corporal y sobre cómo se expresan los genes individuales.
Paralelamente, Kaufman y colaboradores, estudiando mutantes homeóticos que
afectaban a las regiones anteriores del organismo, describieron un segundo complejo, el
Complejo Antennapedia (ANT-C) e hicieron predicciones similares a las del BX-C
(LEWIS et al. 1980b; KAUFMAN et al. 1990).
Reglas de la regulación en cis de los genes del BX-C (LEWIS 1978; LEWIS 1998):
1.- Colinealidad (COL): El orden de los loci en el cromosoma es igual al orden en que
esos loci se expresan a lo largo del eje anteroposterior del cuerpo.
2.- Inactivación en cis (CIN): Las mutaciones de pérdida de función asociados con una
región cis-reguladora determinada tienden a inactivar la región siguiente del complejo
en sentido distal.
3.- Sobreexpresión en cis (COE): Las mutaciones de pérdida de función asociados con
una región cis-reguladora determinada tienden a ir acompañados por una sobreexpresión
de la función asociada con la región cis-reguladora proximal a ella.
27
INTRODUCCIÓN
28
INTRODUCIÓN
29
INTRODUCCIÓN
30
INTRODUCIÓN
Pero ¿por qué estos genes han permanecido agrupados y en el mismo orden
lineal durante más de 500 millones de años de evolución? Se considera que la
conservación de estas agrupaciones génicas es reflejo de la necesidad funcional para la
correcta activación y función de estos genes. Entre estas necesidades funcionales
encontraríamos la existencia de regiones reguladoras compartidas entre genes, los
efectos de posición o, más recientemente, la colinealidad temporal.
31
INTRODUCCIÓN
fue el resulado de la aparición de un nuevo gen que inhibía la formación del ala, sino la
adaptación de un gen existente que a través de pequeños cambios había ido modificando
los procesos morfogenéticos que llevan a la formación del ala. Poco a poco, el
descubrimiento de que estos genes son mucho más antiguos y la presencia de dotaciones
complejas en todos los grupos ha ampliado el abanico de hipótesis acerca de su
repercusión sobre la evolución morfológica (DENELL et al. 1996).
Figura 10. Distribución de los genes Hox en los bilateria. El árbol muestra la filogénia del DNA
ribosómico 18S obtenida por AGUINALDO et al. (1997). B, ancestro común a los bilateria
(Urbilateria); D, deuteróstomos; E, ecdisozoos; L, lofotrocozoos; P, protóstomos. Las líneas
horizontales indican datos de cartografía. Las barras blancas verticales delimitan genes o grupos
de genes ortólogos. Las relaciones de ortología inciertas están marcadas con un interrogante, los
rectángulos punteados indican pequeños fragmentos de PCR, los rectángulos sólidos indican
homeobox completos o casi completos, y los rectángulos a rayas indican genes Hox de
artrópodos que muestran evolución rápida. Tomado de DE ROSA et al. (1999).
Figure 10. Distribution of Hox genes in bilaterians. Phylogenetic tree obtained with ribosomic
DNA 18S by AGUINALDO et al. (1997). B, common bilaterian ancestor (Urbilateria); D,
deuterostomes; E, ecdisozoans; L, lofotrocozoans; P, protostomes. Horizontales black lines
indicate mapping data. Vertical white bars delineate orthologous genes or groups of genes.
Uncertain orthology relationships are indicated by question marks, dashed boxes indicate short
PCR fragment, solid boxes indicate complete or almost complete homeoboxes, and striped
colours indicate fast evolving arthropod Hox sequences. From de DE ROSA et al. (1999).
32
INTRODUCIÓN
Aunque la dotación del ancestro de los bilateria era ya compleja (ver apartado
1.6.3) el número y tipo de genes Hox centrales y posteriores difiere entre Ecdisozoos,
Lofotrocozoos y Deuteróstomos (DE ROSA et al. 1999; CARROLL 2000), lo cual sugiere
la relación del número y tipos de genes con las diferencias entre estos grandes grupos.
Incluso la aparición de los genes Hox centrales podría haber coincidido con el origen de
los animales triploblásticos (GARCIA-FERNÀNDEZ 2004).
Dentro de los Artrópodos, la diferencia entre crustáceos, miriápodos e insectos
no es el número o el tipo de genes Hox, sino sus dominios de expresión. El límite
anterior de expresión de los genes Ubx or abdominal-A (abdA) está correlacionada con
transiciones en la morfología de los apéndices entre esto grupos (GRENIER et al. 1997).
En insectos, cambios en la morfología del aparato bucal se han relacionado con cambios
en el patrón de expresión de proboscipedia (pb) (ROGERS et al. 2002). En cordados,
cambios en el área de expresión de genes Hox específicos se han utilizado para explicar
la transición de las aletas de los peces a las extremidades de los anfibios (SHUBIN et al.
1997), o los cambios en el número relativo de vértebras cervicales, torácicas y lumbares
en anfibios, pájaros y mamíferos (BURKE et al. 1995).
33
INTRODUCCIÓN
Figura 11. Jerarquía de cambios evolutivos en los genes Hox según su presunto efecto
morfológico. Sobre los recuadros, se muestran los cambios en los planes corporales asociados a
los cambios en los genes Hox. Bajo los recuadros, probables mecanismos moleculares
responsables de estos cambios. Tomado de GELLON y MCGINNIS (1998).
Figure 11. Hierarchy of evolutionary changes in Hox genes arranged by their presumed
morphological impact. Above the boxes, modifications of body plans associated with changes in
Hox genes. Below the boxes, likely molecular mechanims responsible of such changes. From
GELLON y MCGINNIS (1998).
En otros casos han sido cambios en las proteínas codificadas por los genes Hox
lo que ha modificado su acción sobre los genes diana y, en consecuencia su efecto
fenotípico. Así, la diversificación morfológica de los segmentos torácicos en insectos se
ha producido gracias a la ganancia de nuevos genes diana (WEATHERBEE et al. 1999)
como resultado de la adquisición de un nuevo dominio en una proteína Hox (GALANT y
CARROLL 2002).
Cambios sutiles tienen también efecto fenotípico. El momento y lugar de
expresión de un gen Hox dentro de un dominio pueden ser cruciales para su efecto. Por
un lado la expresión de un mismo gen en una misma célula tendrá efectos diferentes
según el momento en que se produzca, por otro lado pequeños cambios cuantitativos
podrán también producir pequeños cambios fenotípicos. Dentro del mismo género,
34
INTRODUCIÓN
Aunque hace tiempo que sabemos que los genes homeóticos codifican factores
de transcripción capaces de cambiar el destino de un segmento entre vías de desarrollo
alternativas, todavía sabemos muy poco sobre cómo controlan la morfología detallada
de estos segmentos (AKAM 1998). Se creía que los genes Hox determinan los patrones
de expresión de los genes realizadores. Ahora se considera que lo que hacen es
modificar la acción de las vias de señalización, modificando esos patrones de expresión
según la región del organismo. A pesar de la cantidad de trabajos realizados sobre estos
genes todavía sabemos muy poco sobre sus genes diana. Los productos de los genes
Hox parecen ser generalistas versátiles. Actúan en numerosos tipos celulares y tejidos,
donde modulan, unas veces de forma dramática pero habitualmente sutil, una amplia
variedad de procesos de desarrollo (AKAM 1998). Pequeños cambios tanto en el
homeobox como en la secuencia diana pueden, en teoría, modificar la afinidad y
producir pequeños cambios de reconocimiento. Las nuevas herramientas genómicas se
están aplicando para la identificación de los genes diana de los genes Hox. Por un lado
el uso de microarrays para detectar genes con expresión diferencial en cepas “salvajes”
versus mutantes (LEEMANS et al. 2001). Por otro lado, el estudio del funcionamiento y
la evolución de las regiones cis-reguladoras donde se unen los genes Hox, para poder
predecir los genes diana a partir de las secuencias genómicas (EBNER et al. 2005; HERSH
y CARROLL 2005).
35
INTRODUCCIÓN
cuanto a su organización genómica (DE ROSA et al. 1999). Aunque en los últimos años
han aumentado los datos disponibles a este respecto, me centraré aquí en los datos
disponibles cuando se inició este trabajo, en diciembre de 1999 (Figura 10). La
organización cromosómica de los genes Hox ha sido caracterizada solamente en
insectos, cordados, un equinodermo y un nemátodo. No se ha analizado la organización
de estos genes en ningún Lofotrocozoo, a pesar de la gran variedad de phyla incluidos
en este grupo.
Dentro de los equinodermos se han estudiado los genes Hox en dos erizos de
mar. Heliocidaris erythrogramma posee al menos ocho genes Hox que parecen estar
36
INTRODUCIÓN
37
INTRODUCCIÓN
segmento abdominal, mientras que abdA y AbdB dan identidad a los restantes
segmentos abdominales (LEWIS 1978; SÁNCHEZ-HERRERO et al. 1985; DUNCAN 1987;
MARTIN et al. 1995).
De los genes no homeóticos contenidos en el ANT-C, tres son genes con
homeobox que desempeñan funciones importantes en el establecimiento del plan
corporal del embrión: el gen materno bcd codifica un morfógeno importante para
establecer la polaridad anteroposterior (BERLETH et al. 1988), ftz es un gen pair rule
(WEINER et al. 1984) y zen un gen zigótico involucrado en la diferenciación
dorsoventral (RUSHLOW et al. 1987). El gen zen2 contiene un homeobox similar al de
zen, pero no tiene función conocida (RUSHLOW et al. 1987). Genes como ama (SEEGER
et al. 1988) o los Ccp (PULTZ et al. 1988) no están implicados en el control del
desarrollo.
38
INTRODUCIÓN
dominios específicos, como regiones “barrera” que aislan las diferentes regiones
reguladoras.
39
INTRODUCCIÓN
Figura 13. Genes que intervienen en la regulación de los genes Hox: activación por la cascada
de segmentación y mantenimiento del patrón por polycomb y trithorax.
Figure 13. Genes involved in Hox gene regulation: activation by the segmentation cascade and
maintenance of the expression patterns by polycomb and trithorax.
40
INTRODUCIÓN
(LEWIS 1998). Las fenotipos homeóticos de pérdida de función se producen sólo cuando
hay otro gen Hox expresandose en la región, ya sea por que la expresión normal es
solapante o porque la ausencia de un gen desreprime la expresión de otro. Cuando esto
no ocurre, como en el caso de lab y Dfd, no se produce una transformación homeótica
sino una deficiencia estructural, como es la ausencia de involución cefálica en los
mutantes para lab (MCGINNIS y KRUMLAUF 1992).
41
INTRODUCCIÓN
42
INTRODUCIÓN
43
INTRODUCCIÓN
1.8. Objetivos
Los genes Hox codifican factores de transcripción cuya función es especificar la
identidad de los segmentos en el embrión temprano. Estos genes están habitualmente
agrupados y organizados en el mismo orden en el que se expresan a lo largo del eje
anteroposterior del cuerpo de los metazoos. Aunque la conservación de la organización
genómica ha sugerido la existencia de constricciones funcionales actuando sobre ella, la
explicación mecanística de la colinealidad ha sido esquiva y bien puede ser que no
exista una causa única y universal. La desestructuración del complejo en C. elegans y la
presencia de tres roturas en dos linajes de Drosophila ponen en cuestión que esta
organización genómica sea realmente necesaria para su correcto funcionamiento.
- Efectos funcionales de las roturas del HOM-C sobre los genes Hox adyacentes:
4. Determinar si la separación de lab del resto de genes del HOM-C (cambio de
posición) ha tenido algún efecto en la estructura del gen o en la evolución de su
secuencia.
5. Determinar si la rotura ha separado algún elemento regulador de los genes Hox
adyacentes.
6. Determinar si la rotura ha tenido algún efecto en la expresión de los genes Hox
adyacentes.
44
INTRODUCIÓN
1.8. Objetives
Hox genes code for transcription factors whose function is to specify segment
identity in the early embryo. Those genes are usually clustered and arranged in the
same order as they are expressed along the antero-posterior body axis of the metazoan.
Although the conservation of the genomic organization has suggested the existence of
functional constrictions acting over it, the mechanistic explanation of colinearity has
been elusive and it can well be that a unique and universal cause does not exist. The
disassembly of the complex in C. elegans and the presence of three splits in two lineages
of Drosophila call into question if this genomic organization is really necessary for
proper functioning.
The objective of this thesis was to investigate the causes and consequences of
natural HOM-C rearrangements in Drosophila. The concrete objectives are the
following:
45
2. RESULTADOS
RESULTADOS
2.1. Artículo 1
La primera parte de este trabajo describe el descubrimiento de una nueva rotura
del HOM-C en las especies del grupo repleta de Drosophila. Este resultado fue obtenido
previamente al inicio de esta tesis doctoral. Esta rotura ha recolocado el gen más
anterior del complejo, lab, en una posición cromosómica distante pero indistinguible de
la posición citológica de los genes Hox más posteriores, abdA y AbdB.
NEGRE B., J.M. RANZ, F. CASALS, M. CÁCERES y A. RUIZ (2003) A New Split of the
Hox Gene Complex in Drosophila: Relocation and Evolution of the Gene labial.
Molecular Biology and Evolution 20(12):2042-2054.
49
RESULTADOS
50
A New Split of the Hox Gene Complex in Drosophila: Relocation and Evolution
of the Gene labial
Bárbara Negre, José Marı́a Ranz,1 Ferran Casals, Mario Cáceres,2 and Alfredo Ruiz
Departament de Genètica i de Microbiologia, Universitat Autònoma de Barcelona, Bellaterra (Barcelona) Spain
Hox genes encode transcription factors involved in the specification of segment identity in the early metazoan embryo.
These genes are usually clustered and arranged in the same order as they are expressed along the anteroposterior body
axis. This conserved genomic organization has suggested the existence of functional constraints acting on the genome
organization. Partial disassembly of the Hox gene complex (HOM-C) in Caenorhabditis elegans and in two different
Drosophila lineages, however, calls into question whether this cluster organization is absolutely required for proper
function. Here we report a new split of the HOM-C discovered in the species of the Drosophila repleta group, which
relocated the most anterior gene of the complex, lab, to a distant chromosomal site near the two most posterior Hox
genes, abd-A and Abd-B. To investigate the evolutionary consequences of natural rearrangements of the Hox gene
complex, the gene lab has been cloned and sequenced in D. buzzatii, a member of the D. repleta group with the split, and
in D. virilis, a member of a different species group without the split. The results show that the structure of lab in D.
buzzatii is intact and place the breakpoint at least 8 kb from its transcription start site. The nucleotide sequence evolution
of lab in the genus Drosophila has been investigated by means of maximum likelihood methods. No significant variation
has been observed among lineages in the rate of nucleotide substitution or in the nonsynonymous/synonymous
substitution ratio. Seemingly, the relocation of lab has not induced a change in evolution rate or degree of functional
constraint. Nevertheless, further work is needed to ascertain whether the lab-pb split has had any effects on gene
expression.
Introduction
Homeotic (Hox) genes encode transcription factors order): labial (lab), proboscipedia (pb), Hox3, Deformed
that specify the future cellular identity along the ante- (Dfd), Sex comb reduced (Scr), fushi tarazu (ftz), Anten-
roposterior axis in the early metazoan embryo by activating napedia (Antp), Ultrabithorax (Ubx), abdominal-A (abd-
or repressing their target genes (Gehring and Hiromi 1986). A), and Abdominal-B (Abd-B). This was probably also the
In the genome, Hox genes are found clustered and arranged organization in the ancestor of the Drosophila genus,
in the same order as they are expressed along the although in insects Hox3 and ftz have lost their Homeotic
anteroposterior body axis (Lewis 1978; McGinnis and function (Hughes and Kaufman 2002). Two rearrangements
Krumlauf 1992). This colinearity appeared before the of this ancestral Hox complex are known to have occurred
separation between vertebrates and invertebrates, and it is during the evolution of the Drosophila genus (fig. 1). In the
widely conserved across the bilaterian phylogeny, which lineage leading to D. melanogaster a split between Antp and
has suggested the existence of functional constraints acting Ubx separated the genes in two complexes: the Anten-
on their genome organization (Duboule and Morata 1994). napedia complex (ANT-C) (Kaufman, Seeger, and Olsen
Partial disassembly of Hox complexes in Caenorhabditis 1990) and the Bithorax complex (BX-C) (Lewis 1978;
elegans (Ruvkun and Hobert 1998) and Drosophila Duncan 1987). These two complexes are located on the
(vonAllmen et al. 1996; Adams et al. 2000; Ranz, Casals right arm of D. melanogaster chromosome 3, separated by
and Ruiz 2001), together with reported experimental breaks approximately 9.6 Mb of euchromatic sequences. The
of the Drosophila complexes (Hazelrigg and Kaufman ANT-C spans ; 400 kb (Adams et al. 2000) and comprises
1982; Struhl 1984; Tiong, Whittle and Gribbin 1987), calls five Hox genes (lab, pb, Dfd, Scr, and Antp), which are
into question whether this organization is absolutely expressed in those segments that will become the head and
required for proper function. the anterior thorax. The BX-C spans ;350 kb (Martin et al.
A single Homeotic gene cluster (HOM-C) is postulated 1995; Adams et al. 2000) and contains three Hox genes
to have existed in the ancestor of the arthropods (Hughes (Ubx, abd-A, and Abd-B), which are expressed in the
and Kaufman 2002), although most sequences are partial posterior thorax and abdomen. A different arrangement of
and, usually, there are no data about the chromosomal the Hox genes was found in D. virilis, a species of the
localization of the genes (Finnerty and Martindale 1998; Drosophila subgenus (vonAllmen et al. 1996). In this
deRosa et al. 1999; Cook et al. 2001). The 10 Hox genes species Ubx mapped by in situ hybridization with Antp
included in this ancestral complex are (in anteroposterior instead of with abd-A, indicating that in the lineage leading
to this species a disruption of the Hox complex occurred
1
Present address: Department of Genetics, University of Cam- between the genes Ubx and abd-A (fig. 1).
bridge, Cambridge, United Kingdom. The diversity in Hox gene arrangement observed
2
Present address: Department of Human Genetics, Emory Univer- among Drosophila species along with the reported
sity School of Medicine, Atlanta, Georgia. flexibility of its genome (Stone 1962; Hartl and Lozovskaya
Key words: labial, Hox complex, Drosophila buzzatii, chromo- 1994; Powell 1997; Ranz, Casals, and Ruiz 2001) prompted
somal rearrangement, genome evolution.
us to investigate the organization of Hox genes in D. buzzatii
E-mail: [email protected]. and D. repleta, two representative species of the D. repleta
Mol. Biol. Evol. 20(12):2042–2054. 2003 species group of the Drosophila subgenus. The in situ
DOI: 10.1093/molbev/msg238
Molecular Biology and Evolution, Vol. 20, No. 12, hybridization to the salivary gland chromosomes of the
Society for Molecular Biology and Evolution 2003; all rights reserved. Antp, ftz, pb, Ubx, abd-A, and Abd-B genes (Ranz, Segarra,
2042
Evolution of the Drosophila Hox gene labial 2043
FIG. 1.—Phylogenetic reconstruction of Hox gene organization in the genus Drosophila. All known information of Hox gene localization has been
plotted in a phylogenetic tree to establish when each reorganization took place. The evolutionary relationships among species and divergence times are
from different sources (Powell 1997; Ranz, Segarra, and Ruiz 1997; Kwiatowski and Ayala 1999). Hox genes are shown as triangles when the
orientation (59 fi 39) is known; they are otherwise shown as rectangles. Hox genes expressed in the anterior and posterior parts of the embryo
appear in red and blue gradients, respectively. The genes ftz, zen1, and Pxd are also included (Wharton 1942; leCalvez 1953; vonAllmen et al.
1996; Powell 1997; Adams et al. 2000); the gene zen2 found in the Hox complex of D. melanogaster (Adams et al. 2000) is not represented.
Unless a connecting line appears, the order of the genes has been assumed as in D. melanogaster but has not yet been demonstrated. Arrows
indicate known splits in the Hox complex. Vertical lines indicate chromosomal discontinuities in the Hox complex. When known, map
position of genes is shown. Information for D. canalinea and D. immigrans: this work; D. melanogaster: Adams et al. 2000; D. pseudoobscura:
Terol, Perez-Alonso, and Frutos 1991, Randazzo et al. 1994; and D. virilis: Maier, Preiss, and Powell 1990, Maier, Sperlich, and Powell 1993,
vonAllmen et al. 1996, and this work. For all the D. repleta group species, a representative organization obtained from D. buzzatii (Ranz,
Segarra and Ruiz 1997; this work), D. hydei (Maier, Preiss, and Powell 1990; Maier, Sperlich, and Powell 1993; this work), D. mercatorum
(this work), and D. repleta (Ranz, Segarra, and Ruiz 1997; Ranz, Casals, and Ruiz 2001; this work) is shown.
2044 Negre et al.
Table 1
Sequence of Oligonucleotide Primers Used for PCR Amplification of Different Regions of the Gene labial in Different
Drosophila Species
Gene Region Species Forward (59–39) Reverse (59–39) Product Size
Exons 2–3 (homeobox) D. buzzatii TCCAGCACCGAAACTACCTG TTCACGCGCTTCTTCTGCTTCAT 1,122 bp
D. hydei 954 bp
D. mercatorum 725 bp
D. repleta 727 bp
D. virilis 1,172 bp
Exon 1 D. buzzatii ATGATGGACGTAAGCAGCATG CATTCCTCTTGAGTTGCATCCA 1,285 bp
D. virilis 1,261 bp
cDNA D. buzzatii CATCCGTGTCAGCTGCTAAC TTCACGCGCTTCTTCTGCTTCAT 570 bp
D. virilis CCTCGTCATCAACAACATCTG TTCACGCGCTTCTTCTGCTTCAT 548 bp
Intron 1 fragment D. buzzatii CCTGGAGCCATATAATAGGT GCTCGATTAATTGATGGATG ;3.5 kb
59 UTR D. virilis TGCTCGGCCAGCAATCAG AGGCAGTGGCTGTGGCATGA 818 bp
and Ruiz 1997; Ranz, Casals, and Ruiz 2001; Ranz et al. the amplified products. The clone of D. virilis was used to
2003) indicated that both species lack the Antp-Ubx split map lab in D. immigrans and D. canalinea.
present in D. melanogaster and share the split between Ubx Preparations of slides, probe labeling, and hybridiza-
and abd-A with D. virilis (fig. 1). However, the location of tion to salivary gland chromosomes were performed
the lab gene was not yet known in those species. Here, we as previously described (Ranz, Segarra, and Ruiz 1997).
have mapped the lab gene in different Drosophila species Homologous and heterologous hybridizations were per-
and show that a new split has separated this gene from the formed at 378C and 258C, respectively. Micrographs were
other anterior Hox genes and relocated it near abd-A and taken by phase contrast with a Nikon Optiphot-1
Abd-B, which direct the posterior development of the fly. To microscope and a Nikon H-III photomicrographic system
determine if the lab-pb split has altered in any way the at 6003 magnification using EKTAR-25 Kodak film and
structure or evolution of lab, we have cloned and sequenced a blue filter. The chromosome map of D. repleta (Wharton
this gene in D. buzzatii (which presents the split) and also 1942) was used to localize the in situ hybridization signals
partially in D. virilis (which lacks the split). The results in all D. repleta group species and D. canalinea; the
show that the molecular structure of lab is conserved cytological maps of D. immigrans (leCalvez 1953) or D.
between D. buzzatii, D. melanogaster, and D. virilis, and virilis (Gubenko and Evgen’ev 1984) were used in the
that the breakpoint must be located at least 8 kb upstream of remaining cases.
this gene. In addition, they indicate that lab is a relatively
fast evolving gene, although there was no evidence for Southern Hybridization and Library Screening
positive or relaxed selection either before or after the split.
Genomic DNA extraction was performed as described
in Piñol et al. (1988). DNA was digested with different
Materials and Methods restriction enzymes and transferred from the agarose gels to
Drosophila Stocks nylon membranes by capillarity transfer under neutral
The following species and stocks were used (code conditions (Sambrook, Fritsch, and Maniatis 1989). Probes
numbers of the Tucson Drosophila Species Stock Center are were labeled with digoxygenin-11-dUTP by random
given in parenthesis): D. buzzatii (st-1 and j-19; see primer with the DIG DNA Labeling and Detection Kit
Cáceres, Puig, and Ruiz 2001), D. canalinea (15050– (Roche). Hybridization was performed at 428C for 16 h
1201.0), D. hydei (Carboneras, Spain), D. immigrans using standard hybridization buffer with formamide (50%).
(15111–1731.0), D. mercatorum (Comarapa, Bolivia), D. Filters were washed twice for 15 min in 23 SSC 0.1%
repleta (15084–1611.2), and D. virilis (Tokyo, Japan). All sodium dodecyl sulfate (SDS) at room temperature and
stocks are fixed for the standard chromosomal arrangements twice for 20 min in 0.13 SSC 0.1% SDS at 688C.
except line j-19 of D. buzzatii, which is homozygous for A lambda genomic library derived from D. buzzatii
inversion 2j. line j-19 (Cáceres, Puig, and Ruiz 2001) was amplified
following Sambrook, Fritsch, and Maniatis (1989) and used
In Situ Hybridization to Salivary Gland Chromosomes for cloning of lab. Screening was done by plaque
hybridization under the same conditions as for the Southern
Clones used as probes to map by in situ hybridization hybridization. Two probes were used to clone lab in D.
the Antp, Ubx, abd-A, and Abd-B genes have been reported buzzatii, one corresponding to exons 2–3 and the other to
elsewhere (Ranz, Segarra, and Ruiz 1997; Ranz, Casals, and exon 1, separated in D. melanogaster by a large intron of
Ruiz 2001). To localize the lab gene in D. buzzatii, D. hydei, nearly 14 kb. The probe containing exons 2 and 3 was the
D. mercatorum, D. repleta, and D. virilis, homologous same 1,122 bp fragment used in the in situ hybridization
probes were generated in each species by polymerase chain (see above). The exon 1 probe consisted in a 1,285-bp
reaction (PCR) amplification of a 725–1,172-bp segment fragment amplified by PCR using specific primers (table 1).
encompassing exons 2–3 of the gene and including the Sequencing of the PCR product confirmed its homology
homeobox (table 1), followed by cloning and sequencing of with exon 1 of D. melanogaster lab.
Evolution of the Drosophila Hox gene labial 2045
Table 2
Structure of the labial Gene in Three Drosophila Species and Tribolium castaneum
Gene Region D. buzzatiia D. virilisa D. melanogasterb T. castaneumc
d
Promoter Inr þ DPE Inr þ DPE Inr þ DPE TATA
59 UTR 757 bp 738 bp 239 bp 181 bp
Exon 1 1,300 bp 1,276 bp 1,216 bp 630 bp
Intron 1 ;19 kb ND 14 kb 11.5 kb
Exon 2 380 bp 380 bp 416 bp 431 bp
Intron 2 687 bp 737 bp 245 bp —
Exon 3 288 bp .50 bp 258 bp —
39 UTR 619–831 bp ND 647 bp 144 bp
NOTE.—ND: not determined.
a
Data from this work.
b
Accession number NG_000062.
c
Data from Nie et al. (2001).
d
Inr: transcription initiator; DPE: downstream promoter element.
in situ hybridization to the chromosomes of those four including exons 2 and 3, and kj-19/68 and kj-19/89,
species, as well as those of D. canalinea and D. immigrans including exon 1. Because those phages did not contain the
(fig. 1). D. buzzatii, D. hydei, D. mercatorum, and D. whole gene, another phage was obtained (kj-19/77), and the
repleta are representative species of four subgroups of the remaining gap (3.5 kb) was closed by PCR amplification
D. repleta group (Wasserman 1992), and in all of them lab (table 1 and fig. 3). Overall, 12,069 bp of Db lab were
was localized at band G4g-5a of chromosome 2. Likewise, sequenced, corresponding to the entire coding sequence
in D. canalinea, which belongs to a closely related species (CDS), intron 2, the 59 UTR, the 39 UTR, and 6,241 bp from
group (Wasserman 1992), Antp and lab mapped to different the upstream regulatory region (table 2 and fig. 3). Intron 1,
chromosomal sites. However, in D. virilis, a phylogeneti- with an estimated size of ;19 kb, was only partially
cally more distant species, lab mapped to 24E, the same sequenced.
constricted region of chromosome 2 where the genes Antp, In the analysis of the genomic sequence of Db lab
zen, and Ubx have been previously localized (vonAllmen (12,069 bp) no other conserved coding regions were found
et al. 1996). Finally, five Hox genes (abd-A, Abd-B, Antp, apart from those corresponding to the lab gene. However,
lab, and Ubx) were mapped by in situ hybridization to the two transposable element insertions were identified: a new
chromosomes of D. immigrans, an even more distantly copy of the ISBu2 element (Cáceres, Puig, and Ruiz 2001)
related species. In this species, lab, Antp, and Ubx mapped of 738 bp and an ISBu3 element of 993 bp, inserted 5,500
to the same chromosomal site, whereas abd-A and Abd-B bp and 1,663 bp, respectively, upstream from the Db lab
mapped to another location of the same chromosome transcription start site (fig. 3). Both elements are closely
(fig. 1). Therefore, we conclude that the virilis-immigrans related to the ISBu1 element of D. buzzatii (Cáceres, Puig,
Hox gene organization (including the split between Ubx and Ruiz 2001), and all three belong to the mini-me class
and abd-A) appears to be ancestral within the Droso- described by Wilder and Hollocher (2001).
phila subgenus. The newly discovered lab-pb split occur- In D. virilis three genomic fragments of the lab gene
red before the radiation of the D. repleta group, 15–22 (Dv lab) where amplified by PCR and sequenced (table 2):
MYA, but after the divergence between the repleta and a 1,172-bp fragment corresponding to exons 2 and 3 and
virilis species groups, 20–30 MYA (Spicer 1988; Russo, intron 2, a 1,268-bp fragments corresponding to exon 1, and
Takezaki, and Nei 1995). To our knowledge, the split a 1,818-bp fragment containing the 59 UTR. In addition,
lab-pb is unique among the eukaryotes and provides the op- a 548-bp cDNA fragment comprising the area from the end
portunity to investigate the evolutionary consequences of of exon 1 to the beginning of exon 3 was also cloned and
natural Hox gene complex rearrangements. sequenced. Overall, 3,184 bp were sequenced, with only
intron 1, part of exon 3, and the 39 UTR remaining.
The molecular structure of Db lab and Dv lab is similar
Molecular Structure of labial in D. buzzatii
to that of Dm lab, comprising three exons and two introns
and D. virilis
(table 2 and fig. 3). Therefore, the structure of the lab
To determine if the lab-pb split has affected in any way transcription unit has not been altered by the lab-pb split in
the molecular structure of lab, this gene was cloned and D. buzzatii and the exact breakpoint must be located at least
sequenced in D. buzzatii and partially in D. virilis. The gene 8 kb upstream from the lab transcription start point. In
lab of D. melanogaster (Dm lab) has three exons and two contrast, Db lab is 32% larger than Dm lab, mostly because
introns, intron 1 being rather long (table 2 and fig. 3). of changes in noncoding regions. The total length of the
Cloning of lab in D. buzzatii (Db lab) was carried out by CDS is about 4% larger (1,968 bp versus 1,890 bp),
screening a lambda library of D. buzzatii genomic DNA whereas intron 1 is ;5 kb longer and intron 2 and the 59
with two probes corresponding to the two distantly located UTR are more than twice as large in D. buzzatii (table 2). In
exon 1 and exon 2–3 regions (see Materials and Methods). general, the D. virilis gene is more like that of D. buzzatii
Three positive phages coming from these two regions were than that of D. melanogaster (table 2), as expected from
subcloned and partially sequenced (fig. 3): kj-19/21, their phylogenetic relationships (fig. 1). Interestingly, the
Evolution of the Drosophila Hox gene labial 2047
FIG. 3.—Restriction map and molecular organization of the lab region in D. buzzatii. The lines underneath the map represent lambda clones studied
during gene cloning and fragments obtained by PCR amplification. Thick segments of lambda clones represent sequenced regions. The structure of lab
is shown above the map, for both D. buzzatii and D. melanogaster. The arrow indicates the start point and direction of transcription, black boxes
represent coding regions, and empty boxes show untranslated regions (UTRs). Gray boxes represent transposable elements found in D. buzzatii; none
has been described in this region in D. melanogaster. Restriction sites: B, BamHI; E, EcoRI; H, HindIII; P, PstI; S, SalI; X, XbaI.
lab structure differs widely between Drosophila and The open reading frame and translation start site for
Tribolium castaneum (table 2). In the latter species, the Db lab were predicted by identification of consensus
CDS is 1,061 bp long with two exons only, and intron 1, the sequences and homology with Dm lab. To find the promoter
59 UTR, and the 39 UTR are also significantly shorter (Nie and transcription start sites, 1,308 bp of genomic DNA
et al. 2001). sequence from the 59 region of Db lab were aligned with
630 bp from the 59 region of Dm lab (X13104 X12834)
Expression of labial and Identification of Regulatory including the 59 UTR and 400 additional nucleotides
Sequences upstream (fig. 4). The most conserved region between the
two sequences corresponds to the transcription start site
To verify the predicted exon structure of Db lab, total
described for Dm lab (Mlodzik, Fjose, and Gehring 1988;
RNA from 0–12-h-old embryos was amplified by RT-PCR.
Diederich et al. 1989). In Db lab, we find in this conserved
Primers were designed on exon 1 and exon 3 (table 1 and
fig. 3), allowing us to examine both exon junctions at the region the sequence TCAGTC that fits well the proposed
same time. The amplified cDNA fragment had the expected consensus TCAþ1G/TTT/C of the initiator (Inr) of Drosophila
size of 570 bp, and its sequence confirmed precisely the genes, where Aþ1 is the transcription start site (Arkhipova
predicted exon-intron junctions. In certain D. melanogaster 1995). No other sequence similar to the Inr consensus was
lines, an alternative splicing of lab intron 1 that produces found in the 59 upstream region of Db lab. An additional
two mRNAs with 18 nucleotides of difference has been alignment including the 738-bp upstream sequence of Dv
described (Mlodzik, Fjose, and Gehring 1988). This lab and using the ClustalW algorithm (Thompson, Higgins,
alternative splicing is caused by a small duplication at the and Gibson 1994) with a gap extension penalty of 0.05
end of intron 1 that provides a new functional splice indicated also a high identity around both the transcription
acceptor site and gives rise to the longer form of the mRNA (Inr) and the translation (Met) start sites between the three
(FlyBase 2002). The analysis of the genomic and cDNA Drosophila species (fig. 4). There is a high homology
sequences of Db lab indicates that in D. buzzatii there is no between the Inr and positionþ32 in the three species (fig. 4),
alternative splicing and only the short form of mRNA is with the sequence CACG located between positions þ29
present. In D. virilis, the sequenced cDNA fragment also and þ32 being fully conserved. This agrees with the
corresponds to the short mRNA. To identify the end of the characterization of this sequence as the Downstream
39 UTR, polyadenylation signals were searched along the Promoter Element (DPE) in Dm lab (Kutach and Kadonaga
sequence 39 from the stop codon. Four possible polyA 2000). Because no TATA box was found upstream of the
signals were identified between 623 and 825 bp from the transcription start site, we conclude that Db lab and Dv lab
stop codon of Db lab. No conserved sequences were found present a TATA-less promoter with Downstream Promoter
between Dm lab and Db lab on the 39 UTR. Element (DPE), as does Dm lab and all other homeotic
2048 Negre et al.
FIG. 4.—Alignment of 59 UTR and upstream sequences of lab from D. buzzatii (Db), D. virilis (Dv) and D. melanogaster (Dm). Alignment was
performed with ClustalW and revised by hand. The primer used to amplify this region in D. virilis is indicated in lowercase. The Initiator (Inr),
downstream promoter element (DPE), and first codon of the protein (Met) are shown. Shaded regions show conservation of three nucleotides in
a window of four in at least two sequences.
genes described in Drosophila (Nie et al. 2001). Analysis of 82.8%, and a second segment 108 bp long in the 59 region
with the UTRscan software (Pesole and Liuni 1999) of the with 84.3% identity. In addition, downstream of this second
59 UTRs of the three Drosophila species identified internal conserved segment there is a ;750-bp region with 50%–
ribosome entry site (IRES) elements in all of them. These 75% identity. A Blast search with those conserved
IRES elements promote cap-independent translation in sequences from D. buzzatii shows similarity only with the
different genes (Hellen and Sarnow 2001). In the Hox genes same sequence in D. melanogaster.
Ubx and Antp, they regulate translation during development We have compared the position of those noncoding
(Ye et al. 1997). conserved sequences with the regulatory elements de-
Given the time elapsed since the separation of the scribed for Dm lab (Chouinard and Kaufman 1991). The
Sophophora and Drosophila subgenera (40–60 Myr), no conserved regions upstream of lab are inside the 2.4-kb
significant nucleotide identity should be expected in fragment that promotes the initiation of procephalon
noncoding regions in the absence of selective pressure to expression. The small fragment with 50%–65 % identity
conserve the function (O’Neil and Belote 1992). Thus, located 4,800 bp upstream of Db lab corresponds to the 1.2-
conserved structures may be interpreted as functionally kb Dm lab fragment that drives expression in the anterior
important for gene expression. The entire sequence for midgut via a decapentaplegic (DPP)-regulated enhancer
Db lab (12,069 bp; fig. 3) was aligned with Dm lab (lab550 – Marty et al. 2001). Finally, the 180-bp conserved
(NG_000062) with the AVID algorithm and visualized with segment at the end of intron 1 corresponds to a fragment not
the VISTA software (Mayor et al. 2000; Dubchak et al. tested by Chouinard and Kaufman (1991), which might
2000). The resulting graph (fig. 5) highlights those contain the regulatory element required for embryonic
segments with a nucleotide identity .50% between Db peripheral nervous system (PNS) expression.
lab and Dm lab. Conserved regions both in coding and
noncoding sequences are shown relative to their position in
Nucleotide Sequence Comparisons
Db lab. There is a high variability in conservation along the
coding sequence. Several peaks of conserved sequence The coding sequences of Db lab and Dv lab were
appear along exon 1, which could be a consequence of the aligned and compared to that of Dm lab (NM_057265).
folding pattern of the protein. Two peaks of conserved Table 3 shows the number of synonymous and non-
sequence at the end of exon 2 and the beginning of exon 3 synonymous substitutions per site (dS and dN) for pairwise
correspond to the homeobox. In the noncoding regions, two comparisons between the three Drosophila species. These
prominent peaks of conserved sequence show up, one numbers have been estimated using maximum likelihood
segment 180-bp long at the end of intron 1 with an identity methods (Yang and Bielawski 2000; Yang and Swanson
Evolution of the Drosophila Hox gene labial 2049
FIG. 5.—VISTA plot showing similarity between the nucleotide sequences of the lab gene of D. melanogaster and D. buzzatii. Conserved regions
are shown relative to their position in Db lab (horizontal axis), and their percent identities (50%–100%) are indicated on the vertical axes. The location
of coding exons (black rectangles) and 59 and 39 UTRs (empty rectangles) are shown above the profile, where the thin arrow indicates the direction of
transcription and gray boxes represent transposable elements. Peaks with an identity of at least 70% over 100 bp are shaded.
2002) for the entire coding sequence (CDS) as well as for complex. Likewise, the x ratio is not different in the
three separate (non-overlapping) regions: exon 1, exon 2 þ phylogenetic branch leading to D. buzzatii than in the rest of
3 (excluding the homeobox) and the homeobox. The ratio Drosophila, and there is no evidence that positive selection
between the number of nonsynonymous and synonymous or relaxation of selection has occurred after the relocation of
substitutions (x ¼ dN/dS) provides information on the type lab by the lab-pb split.
of selection acting on amino acid sequences. The low Finally, we tested for differences among the three
average x values for the entire gene (0.0696–0.0910) separate coding regions, fitting a series of fixed-sites models
suggest a relatively high degree of functional constraint, but as described by Yang and Swanson (2002). A detailed list
there is wide variation among regions (see below). of the models and their results is given in table 5. We were
To investigate if the molecular evolution of lab had particularly interested in the among-region variation in x
been affected by the lab-pb split in the D. repleta species ratio. To test the hypothesis of no variation, model D (which
group, three hypotheses related to the molecular evolution assumes different rs, j, and x) was compared with model B
of lab were tested by comparing the likelihood of different (which assumes different rs). The log-likelihood difference,
models. First, we tested the constancy of the evolutionary 2(lD lB)¼71.65, is highly significant (df¼1915¼4; P ,
rates by comparing the base model (which assumes 0.001) indicating differences either in the transition/trans-
a molecular clock and a single x ratio for all lineages) version ratio j or the ratio x. A modified D model
with a model assuming no clock. The results were (designated D9) with a fixed j value was compared to model
nonsignificant for the entire CDS and the three separate B, and the resulting log-likelihood difference, 2(lD9 lB) ¼
coding regions (table 4). Second, we tested for variation in 64.36, showed that the variation in x ratio was indeed
the x ratio among lineages by comparing the base ‘‘one significant (df ¼ 1615 ¼ 1; P , 0.001). Similar tests can be
ratio’’ model with a ‘‘free-ratio’’ model, involving as many conducted comparing model E (which assumes different rs,
x ratios as branches in the tree (four). The results again were j and x, and ps) with model C (which assumes different rs
nonsignificant both for the entire CDS and for the three and ps). The difference in log-likelihood, 2(lE lC)¼85.06,
separate coding regions (table 4). Therefore, the evolution- is highly significant (df ¼ 3733 ¼ 4; P , 0.001), indicating
ary rate of lab has not accelerated or decelerated again significant differences either in the transition/trans-
significantly after its separation from the anterior Hox gene version ratio j or the ratio x. If a modified E model
Table 3
Maximum Likelihood Estimates of the Number of Synonymous (dS) and Nonsynonymous (dN) Substitutions per Site in
the Coding Regions of the Gene labial for Pairwise Comparisons Between D. buzzatii, D. virilis, and D. melanogaster
Comparison Regiona Codons t j S dS N dN x
D. buzzatii-D. virilis CDS 545 0.7982 1.9050 424.2 0.8140 1,210.8 0.0741 0.0910
Exon 1 413 0.9009 1.7786 340.3 0.9053 898.7 0.0712 0.0786
Exon 2 þ 3 78 0.7204 2.7602 56.9 0.4878 177.1 0.1605 0.3289
homeobox 54 0.7578 0.4796 20.7 1.9608 141.3 0.0020 0.0010
D. buzzatii-D. melanogaster CDS 595 2.4155 1.5918 419.3 2.6722 1,365.7 0.2318 0.0868
Exon 1 395 3.0631 1.4996 286.7 3.4913 898.3 0.2325 0.0666
Exon 2 þ 3 140 2.0560 1.9818 111.7 1.4648 308.3 0.4028 0.2750
homeobox 60 0.8506 1.1505 23.3 2.1726 156.7 0.0022 0.0010
D. virilis-D. melanogaster CDS 525 2.5782 1.4040 345.4 3.1407 1,229.6 0.2186 0.0696
Exon 1 395 3.6802 1.2985 274.2 4.4807 910.8 0.2470 0.0551
Exon 2 þ 3 76 1.7824 1.8819 49.9 1.5643 178.1 0.3223 0.2061
homeobox 54 1.7454 4.8507 34.7 2.7077 127.3 0.0027 0.0010
NOTE.—t: branch length; j: transition/transversion rate ratio; S: number of synonymous sites; N: number of nonsynonymous sites; x: dN/dS ratio.
a
Exon 2 þ 3 does not include homeobox.
2050 Negre et al.
Table 5
Log-Likelihood Values and Parameter Estimates Under Fixed-Sites Models
Model p l r2 r3 j x
A (homogeneous model) 13 5147.69 1 1 1.6552 0.0971
B (different rs) 15 5125.17 1.569 0.387 1.6433 0.0862
C (different rs and ps) 33 5087.45 1.398 0.277 1.6547 0.0813
D (different rs, j, and x) 19 5089.34 0.942 0.983 j1 ¼ 1.4036; j2 ¼ 0.0853; x1 ¼ 0.2857; x2 ¼ 0.7401;
j3 ¼ 2.5631 x3 ¼ 0.0001
D9 (j fixed to 1.643) 16 5092.99 0.888 0.966 j ¼ 1.643 x1 ¼ 0.0883; x2 ¼ 0.2718;
x3 ¼ 0.0001
E (different rs, j and (x), 37 5044.93 0.726 0.613 j1 ¼ 1.4864; j2 ¼ 0.0809; x1 ¼ 0.3255; x2 ¼ 1.1217;
and ps) j3 ¼ 2.3396 x3 ¼ 0.0001
E9 (j fixed to 1.655) 34 5046.94 0.703 0.670 j ¼ 1.655 x1 ¼ 0.0824; x2 ¼ 0.3107;
x3 ¼ 0.0001
F (separate analysis) 39 5042.39 j1 ¼ 1.4978; j2 ¼ 2.3044; x1 ¼ 0.0766; x2 ¼ 0.3374;
j3 ¼ 1.0623 x3 ¼ 0.0001
NOTE.—p: number of parameters of the model; l: log-likelihood value; r2 and r3: relative branch lengths.
Evolution of the Drosophila Hox gene labial 2051
FIG. 6.—Multiple alignment of proteins coded by lab orthologs from different organisms. D. buzzatii (Db) and D. virilis (Dv) from present work.
GenBank accession numbers: D. melanogaster (Dm, CAB57787), A. gambiae (Ag, AAF91398.1 and AAAB01008961), Tribolium castaneum (Tc,
AAF64147.1 and AAF64148.1), Thermobia domestica (Td, AAD50360.1), Homo sapiens HoxA1 (HsHoxA1, NP_005513.1). Dv, Ag, and Td are
partial sequences. Amino acids in lowercase correspond to degenerate primers. Alignments were performed with ClustalW with default parameters,
except for Ag and Td, which have been aligned individually with Db, increasing the gap extension penalty and included in the global alignment
2052 Negre et al.
conserved between Db lab and Dm lab are also more or less melanogaster revealed by sequence analysis. Genetics
conserved in all other sequences studied, including the 139:1359–1369.
human gene HoxA1. The domain conservation suggests that Cáceres, M., M. Puig, and A. Ruiz. 2001. Molecular character-
the protein functionality relies basically on a few domains, ization of two natural hotspots in the Drosophila buzzatii
leaving the rest of the sequence with more freedom to genome induced by transposon insertions. Genome Res.
11:1353–1364.
change. This fact could explain why minigenes carrying the
Chouinard, S., and T. C. Kaufman. 1991. Control of expression
chicken Hoxb1 gene are able to rescue null mutants for the of the homeotic labial (lab) locus of Drosophila melanogas-
lab gene in D. melanogaster (Lutz et al. 1996). Thus, Labial ter: evidence for positive and negative autogenous regulation.
proteins seem to be built on highly conserved blocks joined Development 113:1267–1280.
by rapidly evolving sequences. For example, the spacer or Cook, C. E., M. L. Smith, M. J. Telford, A. Bastianello, and M.
‘‘gene specific region’’ (Peterson et al. 1999) localized Akam. 2001. Hox genes and the phylogeny of arthropods.
between the YKWM motif and the homeodomain is much Curr. Biol. 11:759–763.
more variable, both in sequence and in length, in the lab gene DeRosa, R., J. K. Grenier, T. Andreeva, C. E. Cook, A. Adoutte,
than in other Hox genes even between Drosophila species. M. Akam, S. B. Carroll, and G. Balavoine. 1999. Hox genes
in brachiopods and priapulids and protostome evolution.
Nature 399:772–776.
Concluding Remarks Devenport, M. P., C. Blass, and P. Eggleston. 2000. Character-
Although Hox gene colinearity and cluster conserva- ization of the Hox gene cluster i the malaria vector mosquito,
tion represents a paradigm in Evo-Devo, in the genus Anopheles gambiae. Evol. Dev. 2:326–339.
Drosophila the Hox cluster has been split at least three Diederich, R. J., V. K. L. Merrill, M. A. Pultz, and T. C.
Kaufman. 1989. Isolation, structure, and expresion of labial
times. The lab-pb split described here is the most recent
a homeotic gene of the Antennapedia Complex involved in
example. Even after moving to a new location, labial Drosophila head development. Genes Dev. 3:399–414.
appears to evolve consistently, showing no alteration of its Dubchak, I., M. Brudno, G. G. Loots, L. S. Pachter, C. Mayor,
molecular evolutionary patterns. A paracentric inversion as E. M. Rubin, and K. A. Frazer. 2000. Active conservation of
well as a gene transposition event might, in principle, noncoding sequences revealed by three-way species compar-
account for the lab-pb disruption and the relocation of lab isons. Genome Res. 10:1304–1306.
near the genes abd-A and Abd-B observed in the repleta Duboule, D., and G. Morata. 1994. Colinearity and functional
group of Drosophila. However, the former possibility hierarchy among genes of the homeotic complexes. Trends
seems more likely given the extensive reorganization of the Genet. 10:358–364.
Drosophila genome produced by chromosomal inversions Duncan, I. 1987. The Bithorax complex. Annu. Rev. Genet.
(Powell 1997; Ranz, Casals, and Ruiz 2001). The molecular 21:285–319.
characterization of the chromosomal regions 59 of lab and Finnerty, J. R., and M. Q. Martindale. 1998. The evolution of the
Hox cluster: insights from outgroups. Curr. Opin. Genet. Dev.
39 of pb in D. buzzatii should help to determine which was
8:681–687.
the molecular mechanism that caused such a split and the FlyBase 2002. The FlyBase database of the Drosophila genome
precise location of the breakpoint between lab and pb. This projects and community literature. Nucleic Acids Res.
study would also elucidate the exact physical relationship, 30:106–108. (http//flybase.bio.indiana.edu/)
if any, between lab and the genes abd-A and Abd-B. Gehring, W. J., M. Affolter, and T. Bürglin. 1993. Homeodomain
Furthermore, although the rearrangement has not altered the proteins. Annu. Rev. Biochem. 63:487–526.
lab coding region, and although some of its regulatory Gehring, W. J., and Y. Hiromi. 1986. Homeotic genes and the
sequences seem to be in place, whether the relocation of lab homeobox. Annu. Rev. Genet. 20:147–173.
had any consequences upon the level or spatiotemporal Gubenko, I. S., and M. B. Evgen’ev. 1984. Cytological and
pattern of expression of this gene remains an open question, linkage maps of Drosophila virilis chromosomes. Genetics
one that we are already investigating. 65:127–139.
Hartl, D. L., and E. R. Lozovskaya. 1994. Genome evolution:
between the nucleosome and the chromosome. Pp 579–592 in
Acknowledgments B. Schierwater, B. Streit, G. P. Wagner, and R. DeSalle, eds.
We thank S. Celniker, J. E. Hooper, and E. Sánchez- Molecular ecology and evolution: approaches and applica-
tions, Birkhauser Verlag, Basel.
Herrero for probes; the National Drosophila Species
Hazelrigg, T., and T. C. Kaufman. 1982. Revertants of dominant
Resource Center (Bowling Green, Ohio, USA) for fly mutations associated with the Antennapedia gene complex of
strains; and J. González and E. Betrán for advice and helpful Drosophila melanogaster: cytology and genetics. Genetics
comments on the manuscript. This work was supported by 105:581–600.
grant PB98-0900-C02-01 from the Dirección General de Hellen, C. U. T., and P. Sarnow. 2001. Internal ribosome entry
Enseñanza Superior e Investigación Cientı́fica (MEC, sites in eukaryotic mRNA molecules. Genes Dev. 15:1593–
Spain) awarded to A.R. and by a doctoral FI/DGR fellow- 1612.
ship from the Generalitat de Catalunya awarded to B.N. Hong, Y. S., S. Y. Kim, A. Bhattacharya, D. R. Pratt, W. K.
Hong, and M. A. Tainsky. 1995. Structure and function of
Literature Cited the Hox A1 human homeobox gene cDNA. Gene 159:209–
214.
Adams, M. D., S. E. Celniker, R. A. Holt et al. (192 co-authors). Hooper, J. E., M. Pérez-Alonso, J. R. Bermingham, M. Prout,
2000. The genome sequence of Drosophila melanogaster. B. A. Rocklein, M. Wagenbach, J. E. Edstrom, R. de Frutos,
Science 287:2185–2195. and M. P. Scott. 1992. Comparative studies of Drosophila
Arkhipova, I. R. 1995. Promoter elements in Drosophila Antennapedia genes. Genetics 132:453–569.
Evolution of the Drosophila Hox gene labial 2053
Hughes, C., and T. C. Kaufman. 2002. Hox genes and the Powell, J. R. 1997. Progress and prospects in evolutionary biology:
evolution of the arthropod body plan. Evol. Dev. 4:459–499. the Drosophila model. Oxford University Press, New York.
Johnson, F. B., E. Parker, and M. A. Krasnow. 1995. Extradenticle Powers T. P., J. Hogan, Z. Ke, K. Dymbrowski, X. Wang, F. H.
protein is a selective cofactor for the Drosophila hometics: role Collins, and T. C. Kaufman. 2000. Characterization of the
of the homeodomain and YPWM amino acid motif in the Hox cluster from the mosquito Anopheles gambiae (Diptera:
interaction. Proc. Natl. Acad. Sci. USA 92:739–743. Culicidae). Evol. Dev. 2:311–325.
Kaufman, T. C., M. A. Seeger, and G. Olsen. 1990. Molecular Randazzo, F. M., M. A. Seeger, C. A. Huss, M. A. Sweeney,
and genetic organization of the Antennapedia gene complex J. K. Cecil, and T. C. Kaufman. 1994. Structural changes
of Drosophila melanogaster. Adv. Gen. 27:309–362. in the Antennapedia complex of Drosophila pseudoobscura.
Kutach A. K., and J. T. Kadonaga. 2000. The downstream Genetics 133:319–330.
promoter element DPE appears to be as widely used as the Ranz, J. M., F. Casals, and A. Ruiz. 2001. How maleable is the
TATA Box in Drosophila core promoters. Mol. Cell Biol. eukaryotic genome? Extreme rate of cromosomal rearrange-
20:4754–4764. ment in the genus Drosophila. Genome Res. 11:
Kwiatowski, J., and F. J. Ayala. 1999. Phylogeny of Drosophila 230–239.
and related genera: conflict between molecular and anatomical Ranz, J. M., J. González, F. Casals, and A. Ruiz. 2003. Low
analyses. Mol. Phylogenet. Evol. 13:319–328. occurrence of gene transposition events during the evolution
LeCalvez. 1953. Carte cytologique des chromosomes geants des of the genus Drosophila. Evolution 57:1325–1335.
glandes salivaires de Drosophila immigrans Sturt. Chromo- Ranz, J. M., C. Segarra, and A. Ruiz. 1997. Chromosomal
soma 6:170–174. homology and molecular organization of Muller’s elements D
Lewis, E. B. 1978. A gene complex controlling segmentation in and E in the Drosophila repleta species group. Genetics
Drosophila. Nature 276:565–570. 145:281–295.
Lutz, B., H. C. Lu, G. Miller, and T. C. Kaufman. 1996. Rescue Rauskolb C., and E. Wieschaus. 1994. Coordinate regulation of
of the Drosophila labial null mutant by chicken ortholog downstream genes by extradenticle and the honeotic selector
Hoxb-1 demonstrates that the function of Hox genes is proteins. EMBO J. 13:3561–3569.
phylogenetically conserved. Genes Dev. 10:176–184. Russo, C. A. M., N. Takezaki, and M. Nei. 1995. Molecular
Maier, D., A. Preiss, and J. R. Powell. 1990. Regulation of the phylogeny and divergence times of Drosophilid species. Mol.
segmentation gene fushi tarazu has been functionally Biol. Evol. 12:391–404.
conserved in Drosophila. EMBO J. 9:3957–3966. Ruvkun, G., and O. Hobert. 1998. The taxonomy of de-
Maier, D., D. Sperlich, and J. R. Powell. 1993. Conservation and velopmental control in Caenorhabditis elegans. Science
change of the developmentally crucial fushi tarazu gene in 282:2033–2041.
Drosophila. J. Mol. Evol. 36:315–326. Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular
Martin, C. H., C. A. Mayeda, C. A. Davis, C. L. Ericsson, J. D. cloning: a laboratory manual, 2nd edition. Cold Spring Harbor
Knafels, D. R. Mathog, S. E. Celniker, E. B. Lewis, and M. J. Laboratory Press, Cold Spring Harbor, N.Y.
Palazzolo. 1995. Complete sequence of the bithorax complex Shippy, T. D., S. J. Brown, and R. E. Denell. 1998. Molecular
of Drosophila. Proc. Natl. Acad. Sci. USA 92:8398–8402. characterization of the Tribolium abdominal-A ortholog and
Marty, T., M. A. Vigano, C. Ribeiro, U. Nussbaumer, N. C. implications for the products of the Drosophila gene. Dev.
Grieder, and M. Affolter. 2001. A HOX complex, a repressor Genes Evol. 207:446–452.
element and a 50 bp sequence confer regional specificity to Spicer, G. S. 1988. Molecular evolution among some Drosophila
a DPP-responsive enhancer. Development 128:2833–2845. species groups as indicated by two-dimensional electropho-
Mayor, C., M. Brudno, J. R. Schwartz, A. Poliakov, E. M. Rubin, resis. J. Mol. Evol. 27:250–260.
K. A. Frazer, L. S. Pachter, and I. Dubchak. 2000. VISTA: Stone, W. S. 1962. The dominance of natural selection and the
visualizing global DNA sequence alignments of arbitrary reality of superspecies (species groups) in the evolution of
length. Bioinformatics 16:1046–1047. Drosophila. Univ. Texas Publ. 6205:507–537.
McGinnis, W., and R. Krumlauf. 1992. Homeobox genes and Struhl, G. 1984. Splitting the bithorax complex of Drosophila.
axial patterning. Cell 68:283–302. Nature 308:454–457.
Mlodzik, M., A. Fjose, and W. J. Gehring. 1988. Molecular Terol, J., M. Perez-Alonso, and R. Frutos. 1991. Molecular
structure and spatial expression of a homeobox gene from the characterization of the zerknüllt region of the Antennapedia
labial region of the Antennapedia-complex. EMBO J. complex of D. subobscura. Hereditas 114:131–139.
7:2569–2578. Tiong, S. Y., J. R. Whittle, and M. C. Gribbin. 1987.
Nie, W., B. Stronach, G. Panganiban, T. Shippy, S. Brown, and Chromosomal continuity in the abdominal region of the
R. Denell. 2001. Molecular characterization of Tclabial and bithorax complex of Drosophila is not essential for its
the 39 end of the Tribolium homeotic complex. Dev. Genes contribution to metameric identity. Development 101:
Evol. 211:244–251. 135–142.
O’Neil, M. T., and J. M. Belote. 1992. Interspecific comparison Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994.
of the transformer gene of Drosophila reveals an unusually CLUSTAL W: improving the sensitivity of progressive
high degree of evolutionary divergence. Genetics 131:113– multiple sequence alignment through sequence weighting,
128. position-specific gap penalties and weight matrix choice.
Pesole, G., and S. Liuni. 1999. Internet resources for the Nucleic Acids Res. 22:4673–4680.
functional analysis of 59 and 39 untranslated regions of Von Allmen, G., I. Hogga, A. Spierer, F. Karch, W. Bender, H.
eukaryotic mRNAs. Trends Genet. 15:378. Gyurkovics, and E. Lewis. 1996. Splits in fruitfly Hox gene
Peterson, M. D., B. T. Rogers, A. Popadić, and T. C. Kaufman. complexes. Nature 380:116.
1999. The embryonic expression pattern of labial, posterior Wasserman, M. 1992. Cytological evolution of the Drosophila
homeotic complex genes and the teashirt homologue in an repleta species group. Pp. 455–552 in C. B. Krimbas and J .R.
apterygot insect. Dev. Genes Evol. 209:77–90. Powell, eds. Drosophila inversion polymorphism. CRC Press,
Piñol, J., A. Francino, A. Fontdevila, and O. Cabré. 1988. Rapid Boca Raton, Fla.
isolation of Drosophila high molecular weight DNA to obtain Wharton, L. T. 1942. Analysis of the repleta group of
genomic libraries. Nucleic Acids Res. 16:2736–2737. Drosophila. Univ. Texas Publ. 4228:23–59.
2054 Negre et al.
Wilder, J., and H. Hollocher. 2001. Mobile elements and the selective pressures among site classes. Mol. Biol. Evol.
genesis of microsatellites in Dipterans. Mol. Biol. Evol. 19:49–57.
18:384–392. Ye, X., P. Fong, N. Iizuka, D. Choate, and D. R. Cavener. 1997.
Yang, Z. 1997. PAML: a program package for phylogenetic Ultrabithorax and Antennapedia 59 untranslated regions
analysis by maximum likelihood. Comp. Appl. Biosci. promote developmentally regulated internal translation initi-
13:555–556. ation. Mol. Cell Biol. 17:1714–1721.
Yang, Z., and J. P. Bielawski. 2000. Statistical methods for
detecting molecular adaptation. Trends Evol. 15:496–503. William R. Jeffery, Associate Editor
Yang, Z., and W. J. Swanson. 2002. Codon substitution models
to detect adaptive evolution that account for heterogeneous Accepted July 14, 2003
RESULTADOS
64
RESULTADOS
2.2. Artículo 2
Este segundo trabajo va dirigido a resolver dos cuestiones. La primera se centra
en deteminar de forma precisa los puntos en los que se rompió el complejo de genes
Hox y el mecanismo que produjo las roturas, y la segunda en determinar si estas roturas
han tenido algún efecto funcional sobre los genes Hox adyacentes a ellas.
Primero, hemos secuenciado dos regiones del genoma de D. buzzatii. Una
contiene los genes lab y abdA y la otra incluye el gen pb, con un total de 257 kb. Se han
anotado los genes presentes en estas regiones y se ha comparado su organización con la
de D. melanogaster y D. pseudoobscura para localizar con precisión los puntos de
rotura. Se ha observado la presencia de bloques conservados en la secuencia no
codificante y se ha comparado con las regiones reguladoras conocidas. Las roturas se
han producido respetando los genes completos, es decir, la unidad de transcripción más
sus regiones cis-reguladoras. La posición de las roturas se ha reducido hasta un margen
de unas 3 kb. Todo indica que el mecanismo que produjo las roturas són la inversiones
paracéntricas.
La elevada presencia de bloques conservados en la secuencia no codificante y la
conservación de las regiones reguladoras apuntan a que los patrones de expresión de
estos genes deberían estar conservados. No parece haber habido pérdida de regiones
reguladoras con las roturas, pero se podrían haber adquirido nuevas secuencias que
modificaran el patrón de expresión. Para contrastar esta posibilidad, se ha analizado la
expresión de los tres genes Hox, lab, pb y abdA, en embriones y discos imaginales de
cuatro especies de Drosophila portadoras de diferentes organizaciones del complejo.
Los resultados muestran que las roturas del HOM-C no han alterado los patrones de
expresión de los genes Hox adyacentes a ellas. Concluímos que, en Drosophila, la
agrupación de los genes Hox no es imprescindible para su correcto funcionamiento.
Mas bién la organización de los genes Hox parece ser modular y su agrupación
genómica el resultado de la inercia evolutiva.
65
RESULTADOS
66
Letter
Homeotic (Hox) genes are usually clustered and arranged in the same order as they are expressed along the
anteroposterior body axis of metazoans. The mechanistic explanation for this colinearity has been elusive, and it may
well be that a single and universal cause does not exist. The Hox-gene complex (HOM-C) has been rearranged
differently in several Drosophila species, producing a striking diversity of Hox gene organizations. We investigated the
genomic and functional consequences of the two HOM-C splits present in Drosophila buzzatii. Firstly, we sequenced two
regions of the D. buzzatii genome, one containing the genes labial and abdominal A, and another one including
proboscipedia, and compared their organization with that of D. melanogaster and D. pseudoobscura in order to map
precisely the two splits. Then, a plethora of conserved noncoding sequences, which are putative enhancers, were
identified around the three Hox genes closer to the splits. The position and order of these enhancers are conserved,
with minor exceptions, between the three Drosophila species. Finally, we analyzed the expression patterns of the same
three genes in embryos and imaginal discs of four Drosophila species with different Hox-gene organizations. The results
show that their expression patterns are conserved despite the HOM-C splits. We conclude that, in Drosophila, Hox-gene
clustering is not an absolute requirement for proper function. Rather, the organization of Hox genes is modular, and
their clustering seems the result of phylogenetic inertia more than functional necessity.
[Supplemental material is available online at www.genome.org. The sequence data from this study have been
submitted to GenBank under accession nos. AY900631–AY900632 and AY897430–AY897434.]
Homeotic (Hox) genes were discovered in Drosophila melanogaster cestral HOM-C of arthropods (Cook et al. 2001; Hughes and Kauf-
as mutations that transform one body part into another. Lewis man 2002; Hughes et al. 2004). In winged insects, including Dro-
(1978) and Kaufman et al. (1980) found that these genes are sophila, the genes Hox3 and fushi tarazu (ftz) lost their homeotic
clustered and arranged in the chromosome in the same order as function, and thus, only eight truly homeotic genes remain.
their domains of action in the body of flies. Homologous Hox Three different splits of the ancestral HOM-C have been found so
genes were subsequently found in many other animals and their far in the Drosophila genus (Fig. 1A). In D. melanogaster, the com-
arrangement in complexes (HOM-C) shown to be the general plex is split between the genes Antennapedia (Antp) and Ultrabi-
rule (McGinnis and Krumlauf 1992; Ruddle et al. 1994). Hox thorax (Ubx), leaving two separate gene clusters as follows: the
genes encode transcription factors involved in the determination Antennapedia complex, ANT-C (Kaufman et al. 1990) that speci-
of segment identity along the anteroposterior body axis, and fies the identity of the mouth parts and anterior thorax, and the
thus, play a fundamental role in animal development. The con- Bithorax complex, BX-C (Duncan 1987; Martin et al. 1995) in-
served colinearity between Hox gene chromosomal arrangement volved in the development of the posterior thorax and abdomen.
and expression domain is a basic notion of developmental biol- In D. pseudoobscura, the HOM-C is also similarly divided in the
ogy, yet this is an enigmatic phenomenon for which no single ANT-C and BX-C complexes (Lewis et al. 2003). A different split
satisfactory explanation exists (Kmita and Duboule 2003). Fur- between Ubx and abdominal A (abdA) occurs in D. virilis (Von
thermore, HOM-C splits have been observed in Drosophila (Von Allmen et al. 1996), D. repleta (Ranz et al. 2001), D. buzzatii, and
Allmen et al. 1996; Lewis et al. 2003; Negre et al. 2003), Bombyx other species of the Drosophila subgenus (Negre et al. 2003; Fig.
(Yasukochi et al. 2004), nematodes (Aboobaker and Blaxter 1B). Finally, an additional split, between labial (lab) and probosci-
2003), and tunicates (Ikuta et al. 2004; Seo et al. 2004). pedia (pb), is present in D. buzzatii and other species of the repleta
Ten genes arranged in a single complex comprised the an- group (Negre et al. 2003). This third split separated the gene lab
far from pb and the anterior genes of the Hox complex and relo-
5
Corresponding author. cated it near the posterior genes abdA and Abdominal B (AbdB) in
E-mail [email protected]; fax 0034-93-581-23-87.
Article and publication are at http://www.genome.org/cgi/doi/10.1101/ a flagrant violation of the colinearity rule. The functional conse-
gr.3468605. quences of these splits are unknown.
692 Genome Research 15:692–700 ©2005 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/05; www.genome.org
www.genome.org
Drosophila Hox gene complex splits
Results
Molecular characterization of Hox-gene complex breakpoints
To characterize the two HOM-C splits present in D. buzzatii,
we isolated and sequenced two BAC clones, one (5H14, 124,024
bp) containing the lab-abdA region, and another (40C11,
132,938 bp) including the pb region (see Methods). The organi-
zation of the two regions of D. buzzatii chromosome 2 is shown
in Figure 2 along with the homologous regions of D. melanogaster
and D. pseudoobscura for comparison. D. melanogaster and D.
pseudoobscura are homosequential in the analyzed regions,
except where indicated. The sequenced pb region (Fig. 2A)
contains 16 ORFs including Dbuz\pb, Dbuz\zerknüllt
(Dbuz\zen), Dbuz\zerknüllt-related (Dbuz\zen2), and
Dbuz\bicoid (Dbuz\bcd). These four genes are present in the
ANT-C of D. melanogaster and also in the homologous region of
D. pseudoobscura (Fig. 2B,C). The orientation of Dbuz\zen2 is
Figure 1. Genomic (A) and phylogenetic (B) localization of the three
Hox gene complex splits observed in the Drosophila genus. (A) Ancestral the same as that of Dpse\zen2, but inverted with regard to
arrangement of the eight Hox genes within the insects is as follows: labial Dmel\zen2. The remaining 12 genes in this region are ortholo-
(lab), proboscipedia (pb), Deformed (Dfd), Sex combs reduced (Scr), Anten- gous to D. melanogaster genes from four different regions (84D1–
napedia (Antp), Ultrabithorax (Ubx), abdominal A (abdA), and Abdominal B 2, 89D2, 84E5, and 91D4–5) of chromosomal arm 3R. One of the
(AbdB). (B) Phylogenetic relationships and divergence times for the five
Drosophila species included in this study. D. melanogaster and D. pseudo- genes, CG14609, is represented by six copies, in contrast to the
obscura belong to the Sophophora subgenus. D. repleta and D. buzzatii single copy present in D. melanogaster or D. pseudoobscura. A
(both in the repleta species group) and D. virilis (virilis species group) total of four breakpoints are fixed in this region between D. buz-
belong to the Drosophila subgenus (see Negre et al. 2003 for details). zatii and D. melanogaster beside the zen2 microrearrangement.
That corresponding to the lab-pb split is located in the ∼3-kb
intergenic segment between Dbuz\pb and Dbuz\CG17836
In order to ascertain the consequences of Drosophila HOM-C (Fig. 2A).
splits, we have carried out a genomic and functional character- The sequenced lab-abdA region contains 11 ORFs, including
ization of the two splits present in D. buzzatii. We isolated and Dbuz\lab, the cuticular cluster genes (Dbuz\Ccp), and Dbuz\abdA
sequenced two BAC clones containing
the lab-abdA and pb chromosomal re-
gions of D. buzzatii. The gene organiza-
tion in these regions is compared with
that of the homologous regions in D.
melanogaster and D. pseudoobscura to
map the precise site of the two splits.
None of the two splits has altered the
coding regions of Hox genes. We then
searched for Conserved Noncoding Se-
quences (CNS), which are putative regu-
latory sequences, around the genes lab,
pb, and abdA, to find out whether the
splits removed or altered any Hox-gene
enhancer. The position of CNS around
Hox genes is compared with experimen-
tally identified Hox-gene enhancers, and
the arrangement of CNS is compared
between Hox and non-Hox genes. Fi-
nally, we analyzed the expression pat-
terns of three Hox genes, lab, pb, and
abdA, in four Drosophila species with dif- Figure 2. Gene organization of the lab-abdA and pb genomic regions of D. buzzatii compared with
ferent Hox-gene organizations (with and the homologous regions of D. melanogaster and D. pseudoobscura. The localization of the lab-pb split
without the splits) in whole-mount em- (arrow) and the Ubx-abdA split (large arrowhead) are indicated. (A) Sequence of D. buzzatii BAC 40C11
containing the pb region. (B) Organization of the lab-pb region in D. melanogaster. (C) Idem in D.
bryos and imaginal discs. The results
pseudoobscura. (D) Sequence of D. buzzatii BAC 5H14 containing the lab-abdA region. (E) Organization
show that, in Drosophila species, Hox of the abdA region in D. melanogaster and D. pseudoobscura. Genes are represented as open (UTRs) and
genes, as well as their regulatory regions filled boxes (coding sequences) with arrows indicating the sense of transcription. Hox genes are colored
and expression patterns, are conserved, in dark blue, Hox-derived genes in light blue, non-Hox genes in red, noncoding RNA genes in orange,
despite the Hox complex breaks. Thus, and the BcDNA:LP03188 and orthologous sequences in green. Transposable element insertions (usu-
ally ISBu elements, see Negre et al. 2003) are shown as yellow boxes. Large shaded rectangles include
the functional significance of the Hox- homologous Hox-gene regions in different species. Ochre triangles denote small inversions and inser-
gene clustering in Drosophila is question- tions or deletions. Small arrowheads show breakpoints between D. buzzatii and D. melanogaster in
able. non-Hox regions.
(Fig. 2D). The number of Ccp copies (including the gene Edg) is Conserved noncoding sequences in Hox gene regions
eight in the three species, but there is a small inversion encom-
passing two copies (plus the cDNA BcDNA:LP03188) in D. mela- We analyzed the conservation of noncoding sequences around
nogaster in comparison to D. buzzatii or D. pseudoobscura, as well the three Hox genes lab, pb, and abdA by comparing the se-
as one gain and one loss (Fig. 2B–D). These 11 genes come from quences of the three species D. buzzatii, D. melanogaster, and D.
three different regions (84A2–5, 86E11–13, and 89E2) of D. me- pseudoobscura as done previously by other authors (Bergman and
lanogaster chromosomal arm 3R, which means two fixed break- Kreitman 2001; Bergman et al. 2002) (see Methods). Figure 3
points between D. buzzatii and D. melanogaster, beside the small shows the VISTA graph, where the conservation between the
inversion of Ccp genes. One breakpoint corresponds to the lab-pb aligned sequences is plotted (when higher than 50%) and the
split and is found ∼40 kb upstream of Dbuz\lab, in the 5-kb regions that meet the selected criteria (75% identity in a 25-bp
between the sequence similar to BcDNA:LP03188 and the gene window) are highlighted for both coding and noncoding se-
Dbuz\CG31363. The second breakpoint is that of the Ubx-abdA quences. A preliminary analysis showed no differences between
split and is located between 11 and 15 kb downstream of intergenic and intronic regions, in agreement with previous stud-
Dbuz\abdA. The two breakpoints are separated by a DNA seg- ies (Bergman and Kreitman 2001). Thus, CNS are defined as in-
ment of only ∼22 kb encoding a single gene, Dbuz\CG31363 tergenic (excluding UTRs) or intronic sequences that meet the
(Fig. 2D). above criteria. The characteristics of observed CNS are given in
Figure 3. Nucleotide sequence conservation in the lab-abdA and pb regions between Drosophila species. The three panels in each VISTA plot represent
pairwise comparisons between D. melanogaster and D. pseudoobscura (mel/pse), D. melanogaster and D. buzzatii (mel/buz) and D. pseudoobscura and D.
buzzatii (pse/buz). The x-axis represents D. melanogaster coordinates, and y-axis sequence identity (50%–100%). Gray arrows show the location and
orientation of genes. Conservation in exons and UTRs is shown in dark and light blue, respectively. Pink regions represent CNS. Experimentally identified
regulatory sequences (solid purple bars) or segments with negative results (empty bars) are indicated on top of each plot. Five microinversions detected
in the lab or pb regions are enclosed in blue frames, and the VISTA graphs generated with the inverted sequences shown to the right of the main plots.
VISTA plots for the CG17836-CG14290, CG31363, and CG1288-CG14609-CG2520 regions (adjacent to Hox genes) are shown at the bottom of the figure
for comparison.
Tables 1 and 2, and the results of statistical analysis are shown in the CG31363 gene region, between lab and abdA, and the
Supplemental Tables S1 and S2. CG17836-CG14290 and CG1288-CG2520 regions, near pb (Fig.
When D. buzzatii is compared with D. melanogaster or D. 2). These regions include one, two, and three genes, respectively.
pseudoobscura, 395 and 440 CNS are found, respectively, around The pattern of CNS detected is shown in Figure 3 and summa-
the three Hox genes (Table 1). This gives a density of 4.5 and 5 rized in Table 2. In the comparisons with D. buzzatii, we found
CNS per kilobase, respectively. These conserved blocks show a around 100 CNS (∼2/kb), which represents <8% of noncoding
mean size of 44 bp with 86.5% nucleotide identity and represent sequence. Thus, in these non-Hox regions, a much smaller num-
20%–22% of the analyzed noncoding sequence. When D. mela- ber of CNS is observed and the proportion of sequence in CNS is
nogaster and D. pseudoobscura are compared, 563 CNS are de- also significantly lower than in Hox-gene regions (Supplemental
tected (6.5/kb) with a mean size of 55 bp and an average identity Table S1). In the D. melanogaster–D. pseudoobscura comparison,
of 87.4%. In this comparison, the sequence in CNS represents there are 326 CNS (5.7/kb) which represents a 23% of noncoding
36% of noncoding sequence. In all three comparisons, the three sequence. Thus, in this case, the density is similar between Hox
regions around the Hox genes lab, pb, and abdA are homogeneous and non-Hox-gene regions, but the size of CNS is significantly
with little variation either in CNS density, size, or nucleotide smaller in the latter regions (Supplemental Table S1). Conse-
identity (Supplemental Table S1). It is worth noting that CNS are quently, the proportion of sequence in CNS is also significantly
coincident in all three comparisons (Fig. 3), which means that all lower in the non-Hox-gene regions. It should be noted that non-
CNS detected when comparing D. buzzatii with either D. mela- Hox regions show a significant variation for CNS density and also
nogaster or D. pseudoobscura are also found in the comparison for the proportion of sequence in CNS that is not observed in
between the latter two species. Although most CNS keep colin- Hox-gene regions (Supplemental Table S1). The higher variation
earity (relative position and orientation), we could identify four observed between non-Hox regions is probably due to the het-
microinversions, around 1–2 kb in size. One is located within the erogeneity of the sample from a functional point of view. There
large intron of lab and the other three in introns 2 and 3 of pb is little information available on the function and expression
(Fig. 3). pattern of the six non-Hox genes analyzed, which probably rep-
D. buzzatii is equally distant phylogenetically from either D. resent a mixture of genes with different regulatory needs and
melanogaster or D. pseudoobscura (Fig. 1). The latter two species number of enhancers.
belong to the same subgenus and are phylogenetically closer. We
compared the characteristics of the CNS found in the three pair- Conservation of known regulatory sequences
wise comparisons. As expected, there are no statistical differences
Regulatory sequences of the genes lab, pb, and abdA have been
between the CNS found when comparing D. buzzatii with either
experimentally identified in D. melanogaster (Karch et al. 1985;
D. melanogaster or D. pseudoobscura (Supplemental Table S2). The
Chouinard and Kaufman 1991; Kapoun and Kaufman 1995; Mar-
CNS density and the proportion of sequence in CNS are signifi-
tin et al. 1995). We compared their position with the pattern of
cantly higher when comparing the phylogenetically closer spe-
CNS found around Hox genes. As shown in Figure 3, the regula-
cies D. melanogaster and D. pseudoobscura. Increasing divergence
tory sequences identified in D. melanogaster generally contain or
time does not seem to affect the nucleotide identity of the CNS,
correspond to CNS in D. buzzatii. For instance, CNS are found in
although the size of the CNS detected in the Hox-gene regions
the sites corresponding to the iab2 PRE and iab2(1.7) enhancers
shows a significant decrease (Supplemental Table S2).
of abdA (Shimell et al. 1994, 2000). Similarly, a prominent con-
servation peak is observed at the site of the lab550 enhancer,
Conserved noncoding sequences in non-Hox gene regions which directs the expression of lab in the embryo midgut (Marty
To find out whether the observed pattern of CNS is a particular et al. 2001). Also, the inverted segment found in the large intron
feature of Hox genes, we also analyzed the presence of CNS in of the lab gene roughly corresponds to the segment responsible
regions of the sequenced BACs adjacent, but unrelated, to Hox for lab expression in the posterior midgut. Sequence details of the
genes. We used the three microsyntenic regions between D. buz- lab550 and iab2(1.7) enhancer and binding site conservation are
zatii, D. melanogaster, and D. pseudoobscura longer than 10 kb, i.e., shown in Supplemental Figure S1. The Homeotic Response Ele-
Table 1. Characteristics of conserved noncoding sequences (CNS) detected with mVISTA in comparisons of Hox gene regions between
D. melanogaster (mel), D. pseudoobscura (pse), and D. buzzatii (buz)
lab 19,227 mel/pse 129 6.71 (0.59) 53.20 (34.03) 87.69 35.69
mel/buz 73 3.80 (0.44) 46.70 (24.35) 87.33 17.73
pse/buz 84 4.37 (0.48) 44.54 (27.05) 88.08 19.46
pb 42,056 mel/pse 265 6.30 (0.39) 55.04 (35.79) 87.29 34.68
mel/buz 196 4.66 (0.33) 41.88 (22.38) 87.09 19.52
pse/buz 215 5.11 (0.35) 42.92 (24.65) 86.38 21.94
abdA 26,043 mel/pse 169 6.49 (0.50) 59.29 (38.43) 87.44 38.45
mel/buz 126 4.84 (0.43) 45.98 (26.15) 86.18 22.25
pse/buz 141 5.41 (0.46) 46.11 (24.97) 85.44 24.97
Total 87,326 mel/pse 563 6.45 (0.27) 55.89 (36.23) 87.42 36.03
Hox gene regions mel/buz 395 4.52 (0.23) 44.08 (24.04) 86.83 19.94
pse/buz 440 5.04 (0.24) 44.25 (25.53) 86.39 22.30
a
Density = number of CNS per kilobase.
Table 2. Characteristics of conserved noncoding sequences (CNS) detected with mVISTA in comparisons of non-Hox gene regions
between D. melanogaster (mel), D. pseudoobscura (pse), and D. buzzatii (buz)
CG1288-CG2520 18,333 mel/pse 127 6.93 (0.61) 43.48 (28.16) 86.31 30.12
mel/buz 65 3.55 (0.44) 45.26 (31.51) 86.30 16.05
pse/buz 67 3.65 (0.45) 42.44 (29.40) 86.81 15.59
CG17836-CG14290 10,921 mel/pse 46 4.21 (0.62) 45.02 (33.27) 82.67 18.96
mel/buz 18 1.65 (0.39) 39.61 (22.67) 82.88 6.53
pse/buz 22 2.01 (0.43) 42.09 (24.34) 84.34 8.48
CG31363 27,510 mel/pse 153 5.56 (0.45) 35.51 (14.95) 87.17 19.75
mel/buz 22 0.80 (0.17) 26.09 (3.94) 82.93 2.09
pse/buz 23 0.84 (0.17) 28.78 (7.70) 83.23 2.41
Total 56,764 mel/pse 326 5.74 (0.32) 39.96 (24.15) 86.09 22.95
non-Hox gene regions mel/buz 105 1.84 (0.18) 40.28 (27.50) 85.27 7.45
pse/buz 112 1.97 (0.19) 39.59 (25.88) 85.76 7.83
a
Density = number of CNS per kilobase.
ment (HOMRE) of the lab550 enhancer contains four binding Hox gene expression patterns
sites; all of them are conserved in the three species. In the
iab2(1.7) enhancer, there are five Hunchback (HB)-binding sites, The conservation of regulatory sequences suggests that splits of
three of which are conserved in the three species, whereas the the HOM-C had no consequences on Hox-gene expression. To
other two vary in position between species. This enhancer also test this prediction, we compared the expression patterns of the
contains a unique Krüppel (KR)-binding site, where point muta- Hox genes lab, pb, and abdA between D. melanogaster, D. virilis, D.
tions in D. melanogaster cause gain-of-expression mutants (Hab1 buzzatii, and D. repleta. These four Drosophila species represent
and Hab2) (Shimell et al. 1994). This binding site is conserved in three different Hox-gene organizations (Figs. 1,2). D. melanogaster
all three species (Supplemental Fig. S1). The conservation be- possess the Antp-Ubx split only, whereas D. virilis has the Ubx-
tween D. melanogaster and D. buzzatii around abdA ends 9 kb (in abdA split instead. Both D. buzzatii and D. repleta present the
D. melanogaster) and 11 kb (in D. buzzatii) downstream of this Ubx-abdA and lab-pb splits. We used in situ hybridization and
gene (Fig. 3). This boundary lies between the iab2 and pbx regu- antibody staining to whole-mount embryos and to imaginal
latory sequences, which control the expression of abdA and Ubx, discs from third instar larvae and prepupae (see Methods). De-
respectively (Karch et al. 1985). We have shown that the iab2 tailed results are given in Figure 4 and Supplemental Figures S2–
region downstream of abdA is conserved in D. buzzatii. We have S6. The expression patterns of the four species closely follow
not sequenced the Ubx region in D. buzzatii, but we assume that those described for D. melanogaster (for review, see Hughes and
the pbx regulatory sequence will conserve its position upstream Kaufman 2002). Interspecific variation was detected only in the
of Ubx, i.e., there are no rearrangements between Ubx and its pb gene, which in D. virilis presents an extra domain in the em-
regulatory sequences (see below). bryo mesoderm (Fig. 4). As this expression domain is not shared
It is worth noting though, that CNS were also found in by D. melanogaster, it is seemingly not related with the lab-pb
fragments not experimentally tested or described as with no ef- split. Although our analysis is qualitative, and slightly quantita-
fect on expression (Fig. 3). This observation suggests that the tive changes or domain changes of a few cells may remain un-
regulation of these genes may be even more complex than cur- detected, it shows that the reorganization of the HOM-C caused
rently envisaged, and that more regulatory modules may be op- no major alterations of the expression patterns of the three Hox
erative in nature than those experimentally identified in the genes adjacent to the splits, in good agreement with the conser-
laboratory. vation of regulatory sequences (see above). Likewise, Bomze and
Figure 4. Expression pattern of pb in embryos. (A–D) stage 11 embryos, (E–H) stage 17 embryos. (A,E) D. melanogaster, (B,F) D. virilis, (C,G) D. buzzatii,
and (D,H) D. repleta. (A–D) Expression on the ectoderm of the maxillary and labial lobes. Later in development (E–H) pb is detected in the derivatives
of the maxillary (white arrowhead) and labial (black arrowhead) lobes, and in the ventral nervous system (boxed area). (A–D) pb expression is detected
in the mesodermal layer of the mandibular segment (black arrow) in all four species. In D. virilis only (B), pb is also expressed in the mesodermal layer
of the maxillary segment (white arrow). The mandibular (mn), maxillary (mx), and labial (lab) segments are shown in A.
López (1994) found that the expression pattern of Ubx in em- a higher degree of conservation in noncoding regions than more
bryos is conserved between D. melanogaster, D. pseudoobscura, D. simple genes with metabolic or housekeeping functions (Berg-
virilis, and D. hydei (a species of the repleta group), despite their man and Kreitman 2001; Bergman et al. 2002; Halligan et al.
different Hox-gene organization (Figure 1). 2004). Moreover, Hox genes are associated with larger noncoding
regions. Hox genes harbor some of the longest introns of any
Drosophila gene (Moriyama et al. 1998) and mean intron size is
Discussion significantly greater in the Hox than in the non-Hox genes ana-
lyzed here (F = 4.69, df = 1, P < 0.05). This observation also fits
zen2 predates the Drosophila radiation
with the notion that the amount of noncoding DNA must be
The zen and bcd genes come from a duplication of Hox3 in the larger in those genes with complex developmental functions in
ancestor of Cyclorraphan flies (Stauber et al. 2002). A second order to harbor the required CRM (Nelson et al. 2004).
duplication of zen gave birth to zen2, which was thought to be a
recent event in D. melanogaster (Randazzo et al. 1993), where it
has no discernible function. However, the existence of Dpse\zen2 HOM-C evolution in Drosophila
and Dbuz\zen2 shows that the zen–zen2 duplication must predate In Drosophila, Hox genes are arranged in the same 5⬘→3⬘ orien-
the divergence of the Sophophora and Drosophila subgenus, and tation (with only one exception, the Deformed gene in D. mela-
that this gene has been kept during at least 40–60 Myr of evolu- nogaster). Their regulatory sequences are usually located up-
tion. Whether this gene is also present in other flies outside of stream of each gene and in the introns. If we look at the three
the Drosophila genus is still unknown. HOM-C splits known in Drosophila, a common pattern arises. As
can be seen in Figure 2, the lab–pb split took place close to the 3⬘
Patterns of conserved noncoding sequence evolution end of pb and far from the lab 5⬘ end. Likewise, the split between
Cis-Regulatory Modules (CRM) are transcription regulatory DNA the genes Ubx and abdA took place near the abdA 3⬘ end and far
segments (from a few hundred base pair to 1 kb in size) that from the Ubx 5⬘ end, in the short space between their respective
control gene expression in higher eukaryotes (Wray et al. 2003). regulatory sequences pbx and iab2. This is approximately the
CRM have a complex structure still not fully understood. They same position where an experimental break that does not affect
contain one or several binding sites for different transcription development has been observed (Struhl 1984), although the de-
factors, which act cooperatively to activate or repress transcrip- ficiencies used in the complementation tests both carry a fraction
tion of the target gene. As CRM are functionally constrained to of the pbx and iab2 regions. Finally, sequence comparison be-
maintain the expression of the target gene, they evolve slower tween D. melanogaster and D. virilis (Lewis et al. 2003) show that
than nonfunctional sequences. Therefore, the conservation of both the insertion of the CG31217 gene and the Antp–Ubx split
noncoding sequences between phylogenetically distant species took place close to the Ubx 3⬘ end, and far from the Antp 5⬘ end
may be used as a guide for identification of regulatory sequences. (results not shown). Thus, all three splits seem to have occurred
Several recent studies (Bergman and Kreitman 2001; Bergman et far from the 5⬘ end of one gene and much closer to the 3⬘ end of
al. 2002; Cooper and Sidow 2003; Nobrega et al. 2003; Santini et the next one, in such a way as to keep in place the regulatory
al. 2003) support the use of comparative sequence analysis and sequences of both genes. In this way, rearrangements did not
characterization of CNS as a useful approach to detect putative alter any of the known regulatory sequences of these Hox genes;
CRM in Drosophila and other organisms. The clustering of previ- this would explain the absence of gene expression changes.
ously characterized transcription-factor binding sites may be also In the repleta group species, the anterior gene lab is located
used for detection of CRM (Berman et al. 2004). However, the near the posterior genes abdA and AbdB. The sequence analysis
absence of high-quality binding data for most Drosophila tran- shows that lab and abdA are only 75 kb apart and show the same
scription factors represent a great current limitation in the wide- orientation. The breakpoint of the lab–pb split occurred at ∼22 kb
spread application of this method. from that of the Ubx–abdA split. None of those splits seem to
We exhaustively searched for CNS around lab, pb, and abdA have affected the regulatory regions of the Hox genes, because the
and around adjacent non-Hox genes by comparing three species expression patterns of lab and abdA are unaffected. Although it is
pairs. A plethora of highly conserved blocks was found surround- intriguing, the proximity between these genes in the D. buzzatii
ing the three Hox genes in the comparison between the phylo- genome seems purely accidental and lacking any functional sig-
genetically distant species D. buzzatii and D. melanogaster or D. nificance.
pseudoobscura (Fig. 1). The proportion of noncoding sequence The most likely mechanisms for the generation of the
included in CNS was 20%–22%. In most cases, these CNS keep HOM-C splits are paracentric inversions (Ranz et al. 2001; Gonza-
their relative position and colinearity, although a few microrear- lez et al. 2002). A plausible reconstruction of HOM-C evolution
rangements were found. The interpretation of these CNS as regu- in the Drosophila subgenus that accounts for the current organi-
latory sequences is supported by the high neutral substitution zation of Hox genes in D. buzzatii is shown in Figure 5. In lower
rate (Moriyama and Gojobori 1992) and intrinsic rate of DNA Dipterans, such as Anopheles gambiae, the eight Hox genes, plus
loss (Petrov et al. 1996; Singh and Petrov 2004) in Drosophila. Hox3 and ftz, are arranged as a single cluster (Powers et al. 2000).
Noncoding sequences are not expected to be conserved between Before the radiation of the Drosophila genus, two transpositions
such distantly related species unless they are functionally con- occurred as follows: the Ccp gene cluster between lab and pb, and
strained. The coincidence between CNS and known enhancers the gene CG31217 between Antp and Ubx (Lewis et al. 2003).
such as iab2 PRE or lab550 (Supplemental Fig. S1) further sup- Also, zen, zen2, and bcd evolved from the Hox3 gene (see above).
ports this interpretation. In the lineage of the Drosophila subgenus, an inversion took place
A lower CNS density was observed around non-Hox genes. with one breakpoint between Ubx and abdA (split 2 in Fig. 1) and
This result fits well with previous observations showing that the other one between CG31363 and an unknown ORF (X). This
genes with complex developmentally regulated expression show HOM-C structure is now present in species of the Drosophila sub-
Figure 5. Reconstruction of the Hox gene complex evolution in the Drosophila subgenus. Genes are shown as arrows when the orientation 5⬘→3⬘ is
known, and as rectangles otherwise. Hox genes are in black, Hox-related genes in gray, and non-Hox genes in white. (A) Lower Dipterans. (B) Before the
radiation of the Drosophila genus. (C) Drosophila subgenus after its separation from that of the Sophophora subgenus. (D) Ancestor of the repleta group.
(E) Present arrangement of Hox genes in Drosophila buzzatii (cf. Fig. 2A,D).
genus outside the repleta group, such as D. virilis (see Fig. 1). A guillon 2003). In nematodes, the pattern of Hox-gene evolution
second inversion, in the ancestor of the repleta group, split the seems indicative of the move to a deterministic developmental
HOM-C between lab and pb (split 3 in Fig. 1). This inversion, mode (Aboobaker and Blaxter 2003). Bombyx embryogenesis,
which relocated lab close to abdA, had one breakpoint between which is difficult to assign to a short or a long germ insect, is
pb and the Ccp cluster genes and the second breakpoint between characterized by a quick development (Davis and Patel 2002).
CG31363 and CG17836. These two genes are not adjacent in the Drosophila is a long germ insect, where all Hox genes are activated
D. melanogaster genome, but we infer that they were so in the almost simultaneously during the cellular blastoderm stage.
ancestor of the Drosophila subgenus. Thus, none of these organisms seems to show temporal colinear-
ity. A common feature between all organisms shown so far to
Do flies have a Hox gene complex? have a split Hox complex seems to be a derived mode of embryo-
Despite the striking conservation of Hox-gene clustering in meta- genesis characterized by a fast early development.
zoans, if we compare two of the most deeply studied organisms, The loss of temporal progression in the activation of Hox
Drosophila and vertebrates, important differences arise (Ferrier genes in a very rapid mode of embryogenesis could be the ulti-
and Minguillon 2003; Santini et al. 2003; Wagner et al. 2003). mate cause for the modular organization of those Hox “clusters,”
Drosophila Hox-gene regions (1) are much larger than those of where modules can be taken apart without loss of function.
vertebrates, e.g., the human HoxA cluster is only 110 kb long, Given the high rate of chromosomal rearrangement in the genus
whereas the D. melanogaster HOM-C spans 665 kb; (2) contain Drosophila (Ranz et al. 2001; Gonzalez et al. 2002), we anticipate
transposable element insertions, which are remarkably absent in that an even greater variety of Hox-gene organizations will be
those of vertebrates; (3) contain also non-Hox genes that are in- discovered when more species are investigated. It is ironical that
serted between the Hox genes, and tandem duplications within Hox-gene colinearity was discovered in Drosophila, an organism
the complex, such as those of the zen-related genes; (4) allow for with a partially disassembled complex, which may be the by-
small inversions of Hox genes, such as Dfd (Randazzo et al. 1993), product of phylogenetic inertia more than that of functional ne-
and non-Hox genes, such as zen2 (Fig. 2); and (5) are split in three cessity.
ways in different lineages, apparently without consequences on
gene expression. These observations suggest a highly dynamic Methods
evolution in Drosophila that contrasts with the compact structure
seen in vertebrates. Thus, the splits of HOM-C in Drosophila in- Flies
dicate a release of functional requirements present in other meta- D. buzzatii stock st-1 was used for construction of a genomic BAC
zoan. library (González et al. 2005). The following species and stocks
Moreover, Drosophila is not the only organism known to were used for gene expression experiments: D. buzzatii (j19), D.
have a split HOM-C. Split Hox-gene complexes were also known repleta (1611.2), D. virilis (Tokyo-Japan), and D. melanogaster
in nematodes, and recently have been described in Bombyx and (Canton S and Oregon R).
tunicates. What do those organisms have in common in addition
to the split HOM-C? Vertebrate development follows a rostral- BAC sequencing
to-caudal temporal progression, and the colinearity of Hox genes The genomic BAC library was screened with probes from the lab,
is not only spatial, but also temporal (the Hox clock) (Kmita and pb, and abdA genes (González et al. 2005). Positive clones were
Duboule 2003). In the tunicate Oikopleura, Hox gene expression used to build physical maps for the lab-abdA and pb chromo-
still evokes spatial colinearity but not temporal (Seo et al. 2004), somal regions, and one BAC clone from each region was chosen
which favors the argument that the constraining force of HOM-C for sequencing. Shotgun sublibraries were constructed for each
structure conservation is temporal colinearity (Ferrier and Min- BAC using the vector TOPO, and enough plasmid clones were se-
quenced by both ends to reach an ∼6⳯ redundancy. Reads were Josefa González, Roger Hoskins, Barret Pfeiffer, Theresa Ren, Ar-
assembled with the PHRED-PHRAD-CONSED software (Ewing mand Sánchez, Chung Li Shu, and three anonymous referees for
and Green 1998; Ewing et al. 1998; Gordon et al. 1998) and help or comments on the manuscript. This work was supported
sequences finished with one round of AUTOFINISH (Gordon et by grants BMC2002-01708 awarded to A.R. and BMC2002-00300
al. 2001), followed by PCR to bridge the remaining gaps. A con- to E.S.-H. from the Dirección General de Enseñanza Superior e
tinuous high-quality sequence (PHRED score >40) was obtained Investigación Cientı́fica (MEC, Spain), an Institutional Grant
for BAC clones 5H14 (124,024 bp), and 40C11 (132,938 bp). Sta- from the Fundación Ramón Areces, a doctoral FI/DGR fellowship
tistic details of the sequencing process are given in Supplemental from the Generalitat de Catalunya awarded to B.N., and a doc-
Table S3. toral FPI fellowship from the Ministerio de Ciencia y Tecnologı́a
(BES-2003–0416) awarded to S.C.
Sequence annotation
Nucleotide sequences were annotated with the aid of GENE-
SCRIPT (Hudek et al. 2003) and ARTEMIS (Berriman and Ruther- References
ford 2003). Predicted ORFs were corroborated with GOFIGURE
Aboobaker, A. and Blaxter, M. 2003. Hox gene evolution in nematodes:
(Khan et al. 2003) for automatic Gene Ontology (Harris et al. Novelty conserved. Curr. Opin. Genet. Dev. 13: 593–598.
2004) annotation, and BLAST (McGinnis and Madden 2004) for Alonso, C.R. and Akam, M. 2003. A Hox gene mutation that triggers
similarity searches. D. buzzatii sequences were compared with nonsense-mediated RNA decay and affects alternative splicing during
those of homologous regions in D. melanogaster (Celniker et al. Drosophila development. Nucleic Acids Res. 31: 3873–3880.
Bergman, C.M. and Kreitman, M. 2001. Analysis of conserved
2002) and D. pseudoobscura (Richards et al. 2005) genomes. D. noncoding DNA in Drosophila reveals similar constraints in
melanogaster sequences used were as follows: AE001572 (ANT-C), intergenic and intronic sequences. Genome Res. 11: 1335–1345.
DMU31961 (BX-C), and AE003692, AE003672, AE003713, Bergman, C., Pfeiffer, B., Rincon-Limas, D., Hoskins, R., Gnirke, A.,
AE003676, and AE003724 (other regions). D. pseudoobscura con- Mungall, C., Wang, A., Kronmiller, B., Pacleb, J., Park, S., et al. 2002.
Assessing the impact of comparative genomic sequence data on the
tigs AADE01000437 (lab), AADE01000149 (pb), AADE01000036 functional annotation of the Drosophila genome. Genome Biol.
(abdA), and AADE01000014, AADE000175, AADE01002495, 3: research0086.
AADE01000322 (non-Hox genes) were identified with Genome Berman, B.P., Pfeiffer, B.D., Laverty, T.R., Salzberg, S.L., Rubin, G.M.,
VISTA (Dubchak et al. 2000) and the regions of interest anno- Eisen, M.B., and Celniker, S.E. 2004. Computational identification of
developmental enhancers: Conservation and function of
tated. transcription factor binding-site clusters in Drosophila melanogaster
and Drosophila pseudoobscura. Genome Biol. 5: R61.
Analysis of regulatory sequences Berriman, M. and Rutherford, K. 2003. Viewing and annotating
sequence data with Artemis. Brief Bioinform. 4: 124–132.
Pairwise alignments of six homologous genomic regions between Bomze, H.M. and López, A.J. 1994. Evolutionary conservation of the
D. buzzatii, D. melanogaster, and D. pseudoobscura were performed structure and expression of alternatively spliced Ultrabithorax
with the AVID global-alignment tool using default parameters isoforms from Drosophila. Genetics 136: 965–977.
(Bray et al. 2003). CNS were identified in the alignments with Bray, N., Dubchak, I., and Pachter, L. 2003. AVID: A global alignment
program. Genome Res. 13: 97–102.
mVISTA (Mayor et al. 2000) using a window size of 25 bp and a Celniker, S.E., Wheeler, D.A., Kronmiller, B., Carlson, J.W., Halpern, A.,
minimum identity of 75%. Statistical tests were carried out to Patel, S., Adams, M., Champe, M., Dugan, S.P., Frise, E., et al. 2002.
compare the characteristics of the CNS found in the different Finishing a whole-genome shotgun: Release 3 of the Drosophila
regions. Comparisons of CNS size distributions, which depart melanogaster euchromatic genome sequence. Genome Biol.
3: research0079.
significantly from normality, were conducted using the G-test Chouinard, S. and Kaufman, T.C. 1991. Control of expression of the
(Sokal and Rohlf 1995). The number of CNS and the proportion homeotic labial (lab) locus of Drosophila melanogaster: Evidence for
of sequence within CNS was scored for 1-kb windows along the both positive and negative autogenous regulation. Development
analyzed regions (masking out exons). The resulting variables 113: 1267–1280.
Cook, C.E., Smith, M.L., Telford, M.J., Bastianello, A., and Akam, M.
(density and percent sequence in CNS) as well as the nucleotide 2001. Hox genes and the phylogeny of the arthropods. Curr. Biol.
identity (per CNS) were tested using ANOVA (Sokal and Rohlf 11: 759–763.
1995). A complete list of CNS detected is provided in Supplemen- Cooper, G.M. and Sidow, A. 2003. Genomic regulatory regions: Insights
tal Table S5. from comparative sequence analysis. Curr. Opin. Genet. Dev.
13: 604–610.
Davis, G.K. and Patel, N.H. 2002. Short, long, and beyond: Molecular
Gene-expression experiments and embryological approaches to insect segmentation. Annu. Rev.
In situ hybridizations and antibody staining were performed to Entomol. 47: 669–699.
Dubchak, I., Brudno, M., Loots, G.G., Pachter, L., Mayor, C., Rubin,
whole-mount embryos and to imaginal discs from third-instar E.M., and Frazer, K.A. 2000. Active conservation of noncoding
larvae and prepupae as described (Alonso and Akam 2003; Su- sequences revealed by three-way species comparisons. Genome Res.
zanne et al. 2003). cDNA clones were obtained for lab from the 10: 1304–1306.
four species, pb from D. buzzatii and D. melanogaster and abdA Duncan, I. 1987. The bithorax complex. Annu. Rev. Genet. 21: 285–319.
Ewing, B. and Green, P. 1998. Base-calling of automated sequencer
from D. buzzatii, D. repleta, and D. virilis as described (Negre et al. traces using phred. II. Error probabilities. Genome Res. 8: 186–194.
2003) (for primers see Supplemental Table S4). Sense and anti- Ewing, B., Hillier, L., Wendl, M.C., and Green, P. 1998. Base-calling of
sense RNA probes were produced as described (Suzanne et al. automated sequencer traces using phred. I. Accuracy assessment.
2003). When no species-specific probe was available, at least two Genome Res. 8: 175–185.
Ferrier, D.E. and Minguillon, C. 2003. Evolution of the Hox/ParaHox
different ones were used in independent experiments, and the gene clusters. Int. J. Dev. Biol. 47: 605–611.
results were always consistent. Specific antibodies against the González, J., Ranz, J.M., and Ruiz, A. 2002. Chromosomal elements
protein were used for abdA (Macias et al. 1990). evolve at different rates in the Drosophila genome. Genetics
161: 1137–1154.
González, J., Nefedov, M., Bosdet, I., Casals, F., Calvete, O., Delprat, A.,
Shin, H., Chiu, R., Mathewson, C., Wye, N., et al. 2005. A
Acknowledgments BAC-based physical map of the Drosophila buzzatii genome. Genome
Res. (in press).
We thank Claudio Alonso, Joaquin Ariño, Anna Barceló, James Gordon, D., Abajian, C., and Green, P. 1998. Consed: A graphical tool
Castelli-Gair, Jose Felix de Celis, Susan Celniker, Ariel Chipman, for sequence finishing. Genome Res. 8: 195–202.
Gordon, D., Desmarais, C., and Green, P. 2001. Automated finishing split of the Hox gene complex in Drosophila: Relocation and
with autofinish. Genome Res. 11: 614–625. evolution of the gene labial. Mol. Biol. Evol. 20: 2042–2054.
Halligan, D.L., Eyre-Walker, A., Andolfatto, P., and Keightley, P.D. 2004. Nelson, C.E., Hersh, B.M., and Carroll, S.B. 2004. The regulatory content
Patterns of evolutionary constraints in intronic and intergenic DNA of intergenic DNA shapes genome architecture. Genome Biol. 5: R25.
of Drosophila. Genome Res. 14: 273–279. Nobrega, M.A., Ovcharenko, I., Afzal, V., and Rubin, E.M. 2003.
Harris, M.A., Clark, J., Ireland, A., Lomax, J., Ashburner, M., Foulger, R., Scanning human gene deserts for long-range enhancers. Science
Eilbeck, K., Lewis, S., Marshall, B., Mungall, C., et al. 2004. The Gene 302: 413.
Ontology (GO) database and informatics resource. Nucleic Acids Res. Petrov, D.A., Lozovskaya, E.R., and Hartl, D.L. 1996. High intrinsic rate
32: D258–D261. of DNA loss in Drosophila. Nature 384: 346–349.
Hudek, A.K., Cheung, J., Boright, A.P., and Scherer, S.W. 2003. Powers, T.P., Hogan, J., Ke, Z., Dymbrowski, K., Wang, X., Collins, F.H.,
Genescript: DNA sequence annotation pipeline. Bioinformatics and Kaufman, T.C. 2000. Characterization of the Hox cluster from
19: 1177–1178. the mosquito Anopheles gambiae (Diptera: Culicidae). Evol. Dev.
Hughes, C.L. and Kaufman, T.C. 2002. Hox genes and the evolution of 2: 311–325.
the arthropod body plan. Evol. Dev. 4: 459–499. Randazzo, F.M., Seeger, M.A., Huss, C.A., Sweeney, M.A., Cecil, J.K., and
Hughes, C.L., Liu, P.Z., and Kaufman, T.C. 2004. Expression patterns of Kaufman, T.C. 1993. Structural changes in the antennapedia
the rogue Hox genes Hox3/zen and fushi tarazu in the apterygote complex of Drosophila pseudoobscura. Genetics 134: 319–330.
insect Thermobia domestica. Evol. Dev. 6: 393–401. Ranz, J.M., Casals, F., and Ruiz, A. 2001. How malleable is the
Ikuta, T., Yoshida, N., Satoh, N., and Saiga, H. 2004. Ciona intestinalis eukaryotic genome? Extreme rate of chromosomal rearrangement in
Hox gene cluster: Its dispersed structure and residual colinear the genus Drosophila. Genome Res. 11: 230–239.
expression in development. Proc. Natl. Acad. Sci. 101: 15118–15123. Richards, S., Liu, Y., Bettencourt, B.R., Hradecky, P., Letovsky, S.,
Kapoun, A.M. and Kaufman, T.C. 1995. A functional analysis of 5⬘, Nielsen, R., Thornton, K., Hubisz, M.J., Chen, R., Meisel, R.P., et al.
intronic and promoter regions of the homeotic gene proboscipedia in 2005. Comparative genome sequencing of Drosophila pseudoobscura:
Drosophila melanogaster. Development 121: 2127–2141. Chromosomal, gene, and cis-element evolution. Genome Res.
Karch, F., Weiffenbach, B., Peifer, M., Bender, W., Duncan, I., Celniker, 15: 1–18.
S., Crosby, M., and Lewis, E.B. 1985. The abdominal region of the Ruddle, F.H., Bartels, J.L., Bentley, K.L., Kappen, C., Murtha, M.T., and
bithorax complex. Cell 43: 81–96. Pendleton, J.W. 1994. Evolution of Hox genes. Annu. Rev. Genet.
Kaufman, T.C., Lewis, R., and Wakimoto, B. 1980. Cytogenetic analysis 28: 423–442.
of chromosome 3 in Drosophila melanogaster: The homoeotic gene Santini, S., Boore, J.L., and Meyer, A. 2003. Evolutionary conservation of
complex in polytene chromosome interval 84A-B. Genetics regulatory elements in vertebrate Hox gene clusters. Genome Res.
94: 115–133. 13: 1111–1122.
Kaufman, T.C., Seeger, M.A., and Olsen, G. 1990. Molecular and genetic Seo, H.C., Edvardsen, R.B., Maeland, A.D., Bjordal, M., Jensen, M.F.,
organization of the antennapedia gene complex of Drosophila Hansen, A., Flaat, M., Weissenbach, J., Lehrach, H., Wincker, P., et
melanogaster. Adv. Genet. 27: 309–362. al. 2004. Hox cluster disintegration with persistent anteroposterior
Khan, S., Situ, G., Decker, K., and Schmidt, C.J. 2003. GoFigure: order of expression in Oikopleura dioica. Nature 431: 67–71.
Automated Gene Ontology annotation. Bioinformatics Shimell, M.J., Simon, J., Bender, W., and O’Connor, M.B. 1994.
19: 2484–2485. Enhancer point mutation results in a homeotic transformation in
Kmita, M. and Duboule, D. 2003. Organizing axes in time and space; 25 Drosophila. Science 264: 968–971.
years of colinear tinkering. Science 301: 331–333. Shimell, M.J., Peterson, A.J., Burr, J., Simon, J.A., and O’Connor, M.B.
Lewis, E.B. 1978. A gene complex controlling segmentation in 2000. Functional analysis of repressor binding sites in the iab-2
Drosophila. Nature 276: 565–570. regulatory region of the abdominal-A homeotic gene. Dev. Biol.
Lewis, E.B., Pfeiffer, B.D., Mathog, D.R., and Celniker, S.E. 2003. 218: 38–52.
Evolution of the homeobox complex in the Diptera. Curr. Biol. Singh, N.D. and Petrov, D.A. 2004. Rapid sequence turnover at an
13: R587–R588. intergenic locus in Drosophila. Mol. Biol. Evol. 21: 670–680.
Macias, A., Casanova, J., and Morata, G. 1990. Expression and Sokal, R.R. and Rohlf, F.J. 1995. Biometry: The principles and practice of
regulation of the abd-A gene of Drosophila. Development statistics in biological research. W.H. Freeman and Co., New York.
110: 1197–1207. Stauber, M., Prell, A., and Schmidt-Ott, U. 2002. A single Hox3 gene
Martin, C.H., Mayeda, C.A., Davis, C.A., Ericsson, C.L., Knafels, J.D., with composite bicoid and zerknullt expression characteristics in
Mathog, D.R., Celniker, S.E., Lewis, E.B., and Palazzolo, M.J. 1995. non-Cyclorrhaphan flies. Proc. Natl. Acad. Sci. 99: 274–279.
Complete sequence of the bithorax complex of Drosophila. Proc. Natl. Struhl, G. 1984. Splitting the bithorax complex of Drosophila. Nature
Acad. Sci. 92: 8398–8402. 308: 454–457.
Marty, T., Vigano, M.A., Ribeiro, C., Nussbaumer, U., Grieder, N.C., and Suzanne, M., Estella, C., Calleja, M., and Sánchez-Herrero, E. 2003. The
Affolter, M. 2001. A Hox complex, a repressor element and a 50 bp hernandez and fernandez genes of Drosophila specify eye and antenna.
sequence confer regional specificity to a DPP-responsive enhancer. Dev. Biol. 260: 465–483.
Development 128: 2833–2845. Von Allmen, G., Hogga, I., Spierer, A., Karch, F., Bender, W.,
Mayor, C., Brudno, M., Schwartz, J.R., Poliakov, A., Rubin, E.M., Frazer, Gyurkovics, H., and Lewis, E. 1996. Splits in fruitfly Hox gene
K.A., Pachter, L.S., and Dubchak, I. 2000. VISTA: Visualizing global complexes. Nature 380: 116.
DNA sequence alignments of arbitrary length. Bioinformatics Wagner, G.P., Amemiya, C., and Ruddle, F. 2003. Hox cluster
16: 1046–1047. duplications and the opportunity for evolutionary novelties. Proc.
McGinnis, W. and Krumlauf, R. 1992. Homeobox genes and axial Natl. Acad. Sci. 100: 14603–14606.
patterning. Cell 68: 283–302. Wray, G.A., Hahn, M.W., Abouheif, E., Balhoff, J.P., Pizer, M., Rockman,
McGinnis, S. and Madden, T.L. 2004. BLAST: At the core of a powerful M.V., and Romano, L.A. 2003. The evolution of transcriptional
and diverse set of sequence analysis tools. Nucleic Acids Res. regulation in eukaryotes. Mol. Biol. Evol. 20: 1377–1419.
32: W20–W25. Yasukochi, Y., Ashakumary, L., Wu, C., Yoshido, A., Nohata, J., Mita, K.,
Moriyama, E.N. and Gojobori, T. 1992. Rates of synonymous and Sahara, K. 2004. Organization of the Hox gene cluster of the
substitution and base composition of nuclear genes in Drosophila. silkworm, Bombyx mori: A split of the Hox cluster in a non-Drosophila
Genetics 130: 855–864. insect. Dev. Genes Evol. 214: 606–614.
Moriyama, E.N., Petrov, D.A., and Hartl, D.L. 1998. Genome size and
intron size in Drosophila. Mol. Biol. Evol. 15: 770–773.
Negre, B., Ranz, J.M., Casals, F., Cáceres, M., and Ruiz, A. 2003. A new Received November 15, 2004; accepted in revised form January 26, 2005.
76
RESULTADOS
SUPPLEMENTAL MATERIAL
Table S1. Statistical analysis of CNS characteristics between regions for the three pair-wise
species comparisons.
Species Nucleotide Sequence in
(1)
Source of variation pair Density Size distribution identity (%) CNS (%)
d.f. F d.f. G d.f. F d.f. F
Between the six regions mel-pse 5 1.74 60 104.34*** 5 4.84*** 5 7.17***
mel-buz 5 15.31*** 60 58.19 5 2.13 5 16.16***
pse-buz 5 18.75*** 60 64.82 5 3.30* 5 18.76***
Between Hox regions mel-pse 2 0.23 24 30.63 2 0.13 2 1.06
mel-buz 2 1.40 24 26.52 2 0.49 2 1.19
pse-buz 2 1.18 24 38.04* 2 4.55* 2 1.56
Between non-Hox regions mel-pse 2 3.34* 24 24.90 2 9.71*** 2 3.64*
mel-buz 2 9.44*** 24 16.68 2 2.69 2 11.71***
pse-buz 2 11.02*** 24 9.65 2 0.40 2 10.52***
Hox vs non-Hox regions mel-pse 1 1.58 12 60.36*** 1 4.54* 1 26.48***
mel-buz 1 54.87*** 12 17.41 1 4.29* 1 55.00 ***
pse-buz 1 69.34*** 12 19.59 1 6.63* 1 69.65***
* P < = 0.05; ** P < 0.01; *** P < 0.001.
(1)
Variation between the six regions (df = 5) can be partitioned into three components: variation between
the three Hox gene regions (df = 2), variation between the three non-Hox gene regions (df = 2) and
variation between Hox and non-Hox gene regions (df = 1).
Table S2. Statistical analysis of CNS characteristics between species pairs for the different
regions.
Nucleotide Sequence in
(1)
Source of variation Region Density Size distribution identity (%) CNS (%)
d.f. F d.f. G d.f. F d.f. F
Between the three Hox genes 2 15.58*** 24 28.23 2 0.70 2 43.97***
species pairs Non-Hox genes 2 49.44*** 24 9.59 2 1.89 2 34.10***
Mel/buz versus pse/buz Hox gene 1 2.09 12 2.96 1 0.08 1 1.69
Non-Hox gene 1 0.10 12 3.60 1 0.19 1 0.03
Mel/pse versus Hox gene 1 29.06*** 12 27.69** 1 1.32 1 86.24***
Mel/buz or pse/buz Non-Hox gene 1 98.78*** 12 5.28 1 3.59 1 68.17***
** P < 0.01; *** P < 0.001.
(1)
Variation between the three species pairs, mel/pse, mel/buz and pse/buz (df = 2) can be partitioned into
variation between the equally distant species pairs mel/buz and pse/buz (df = 1) and variation between the
closely related mel/pse and the other two species pairs (df = 1).
77
RESULTADOS
Table S3. Statistics of the sequencing process of two BAC clones from a Drosophila buzzatii
genomic library.
BAC clone 5H14 40C11
Number of clones in plasmid library (insert size 1.5 kb) 960 1,056
Number of clones sequenced by one end 77 -
Number of clones sequenced by both ends 712 864
Number of reads in shotgun sequencing 1,501 1,728
Coverage 7.7x 6x
Number of contigs (> 2kb) after automatic assembly 7 7
Number of contigs (> 2kb) after manual editing 4 6
Number of contigs (< 2kb) after manual editing 2 0
Number of extra reads from plasmid library (clone ends) 47 43
Number of extra reads from plasmid library (custom primers) 22 31
Number of PCR fragments sequenced (physical gaps) 1 3
Number of reads from PCR fragments 4 11
Number of reads for finishing 73 85
Total number of reads produced 1,574 1,813
Total number of reads in final assembly 1,391 1,396
Average read size (bp) 829 890
Size of BAC clone insert (bp) 124,024 132,938
Error rate (per 10 kb) < 0.01 0.17
78
RESULTADOS
79
RESULTADOS
80
RESULTADOS
81
RESULTADOS
82
RESULTADOS
83
RESULTADOS
84
RESULTADOS
85
RESULTADOS
86
RESULTADOS
87
RESULTADOS
88
RESULTADOS
89
RESULTADOS
labial
D. buzzatii 5H14 D.melanogaster D. pseudoobscura
AE001572 AADE01000437
inicio fin tamaño % inicio fin inicio fin
93374 93409 = 36 al 77.8 pse 58888 58853
93856 93884 = 29 al 75.9 mel 374162 374190
94584 94606 = 23 al 78.3 pse 57871 57849
94627 94656 = 30 al 80.0 mel 374573 374601
94627 94664 = 38 al 84.2 pse 57813 57776
94675 94692 = 18 al 94.4 pse 57771 57754
94683 94699 = 17 al 94.1 mel 374622 374637
94749 94764 = 16 al 93.8 pse 57699 57684
94835 94865 = 31 al 87.1 pse 57636 57606
94891 94915 = 25 al 100.0 mel 374807 374831
94891 94922 = 32 al 96.9 pse 57558 57527
95039 95101 = 66 al 84.8 mel 374880 374945
95039 95109 = 72 al 84.7 pse 57475 57404
95273 95292 = 20 al 90.0 pse 57242 57224
95274 95296 = 23 al 82.6 mel 375131 375152
96975 96997 = 23 al 82.6 pse 56750 56729
97580 97604 = 25 al 76.0 mel 375584 375607
100621 100694 = 75 al 89.3 pse 56160 56086
100635 100728 = 97 al 88.7 mel 376536 376632
100849 100888 = 40 al 87.5 mel 376678 376717
100849 100892 = 44 al 88.6 pse 55975 55932
100943 100981 = 39 al 94.9 pse 55885 55847
100947 100981 = 35 al 97.1 mel 376771 376805
101045 101087 = 43 al 100.0 pse 55755 55713
101050 101087 = 38 al 100.0 mel 376871 376908
101149 101189 = 43 al 83.7 mel 376973 377014
101167 101197 = 31 al 83.9 pse 55625 55596
101200 101227 = 29 al 82.8 mel 377019 377047
101267 101297 = 31 al 87.1 pse 55442 55413
101277 101292 = 16 al 100.0 mel 377163 377178
101416 101434 = 20 al 90.0 pse 55359 55340
101438 101462 = 26 al 88.5 pse 55325 55300
101545 101581 = 38 al 86.8 mel 377931 377968
101604 101627 = 26 al 76.9 mel 378013 378037
101692 101714 = 24 al 83.3 mel 378096 378119
90
RESULTADOS
91
RESULTADOS
92
RESULTADOS
93
RESULTADOS
94
RESULTADOS
abd-A
D. buzzatii 5H14 D.melanogaster D. pseudoobscura
DMU31961 AADE01000036
inicio fin tamaño % inicio fin inicio fin
300 346 = 47 al 80.9 pse 77999 78045
317 343 = 27 al 77.8 mel 152630 152656
404 464 = 61 al 88.5 pse 78115 78175
407 456 = 50 al 94.0 mel 152739 152788
514 562 = 49 al 87.8 pse 78297 78344
518 572 = 56 al 76.8 mel 152855 152909
737 774 = 39 al 79.5 pse 78399 78437
753 780 = 28 al 78.6 mel 152948 152975
970 1052 = 83 al 84.3 mel 152998 153079
1211 1235 = 25 al 96.0 pse 79197 79221
1216 1243 = 28 al 85.7 mel 153102 153129
1363 1389 = 28 al 78.6 pse 79311 79338
1710 1727 = 18 al 100.0 pse 79528 79545
3529 3544 = 17 al 88.2 pse 81189 81205
4112 4168 = 57 al 78.9 mel 155222 155277
4125 4162 = 38 al 84.2 pse 81627 81664
4321 4344 = 24 al 79.2 pse 81880 81903
5232 5242 = 11 al 100.0 mel 156597 156607
5964 6007 = 44 al 86.4 mel 156698 156741
5974 6007 = 34 al 88.2 pse 83110 83143
6068 6129 = 62 al 82.3 pse 83226 83287
6073 6147 = 76 al 84.2 mel 156832 156907
6136 6152 = 17 al 100.0 pse 83316 83332
7240 7259 = 20 al 100.0 mel 158026 158045
7241 7259 = 19 al 94.7 pse 83955 83973
7288 7341 = 55 al 80.0 pse 84033 84087
7294 7340 = 48 al 83.3 mel 158116 158163
7433 7492 = 61 al 83.6 mel 158185 158243
7443 7480 = 38 al 86.8 pse 84203 84238
7558 7618 = 61 al 88.5 pse 84301 84361
7563 7611 = 49 al 87.8 mel 158270 158318
7689 7708 = 20 al 95.0 mel 158355 158373
7727 7763 = 37 al 75.7 mel 158376 158412
7863 7919 = 57 al 78.9 pse 84467 84523
7870 7919 = 52 al 82.7 mel 158454 158505
95
RESULTADOS
96
RESULTADOS
97
RESULTADOS
98
RESULTADOS
99
RESULTADOS
100
RESULTADOS
101
RESULTADOS
CG17836-CG14290
D. buzzatii 40C11 D.melanogaster D. pseudoobscura
AE003724 AADE001000175
inicio fin tamaño % inicio fin inicio fin
69953 69874 = 83 al 81.9 pse 313 395
69794 69762 = 33 al 78.8 pse 475 507
69632 69607 = 26 al 92.3 pse 641 666
69555 69525 = 31 al 80.6 pse 732 762
69512 69448 = 67 al 76.1 pse 769 835
69137 69109 = 29 al 79.3 mel 54961 54989
68556 68531 = 26 al 76.9 pse 1167 1192
67779 67755 = 25 al 76.0 mel 57456 57479
67719 67687 = 33 al 90.9 pse 1919 1949
67717 67662 = 56 al 82.1 mel 57532 57587
67658 67622 = 37 al 78.4 pse 1909 1942
67637 67580 = 62 al 80.6 mel 57597 57658
67597 67550 = 48 al 83.3 pse 2007 2054
67575 67532 = 48 al 77.1 mel 57663 57710
67457 67421 = 37 al 75.7 pse 2182 2217
67161 67047 = 117 al 89.7 pse 2502 2616
67148 67047 = 103 al 87.4 mel 58061 58163
67026 66997 = 30 al 86.7 mel 58187 58216
67026 66997 = 30 al 93.3 pse 2640 2669
66966 66915 = 56 al 78.6 mel 58240 58295
66966 66929 = 39 al 89.7 pse 2706 2744
65976 65954 = 24 al 95.8 mel 59263 59286
65976 65924 = 57 al 82.5 pse 3829 3885
65945 65924 = 22 al 90.9 mel 59300 59321
65885 65851 = 35 al 97.1 pse 3925 3959
65884 65847 = 38 al 94.7 mel 59370 59407
65797 65722 = 78 al 84.6 pse 4005 4081
65791 65723 = 71 al 83.1 mel 59444 59514
64307 64285 = 25 al 76.0 pse 4631 4655
63430 63413 = 18 al 88.9 pse 5181 5198
61461 61450 = 12 al 100.0 mel 62855 62866
61240 61216 = 26 al 76.9 mel 62991 63016
61237 61212 = 27 al 88.9 pse 7929 7955
59356 59327 = 30 al 80.0 pse 9792 9821
59264 59247 = 19 al 94.7 pse 9859 9877
102
RESULTADOS
103
RESULTADOS
CG31363
D. buzzatii 5H14 D.melanogaster D. pseudoobscura
AE003692
inicio fin tamaño % inicio fin inicio fin
59041 59017 = 25 al 80.0 mel 70059 70083 AADE01002495
59041 59008 = 38 al 78.9 pse 11971 11934
58917 58891 = 27 al 77.8 mel 70206 70232
58114 58095 = 21 al 85.7 pse 10251 10231
57090 57064 = 30 al 80.0 mel 75953 75982
56143 56123 = 21 al 100.0 mel 76122 76142
56113 56073 = 41 al 90.2 pse 6029 5989
56091 56067 = 25 al 96.0 mel 76215 76239
56007 55976 = 32 al 81.3 pse 5908 5877
55994 55957 = 38 al 84.2 mel 76302 76339
55607 55585 = 25 al 76.0 mel 76904 76927
55332 55299 = 34 al 91.2 pse 4920 4888
55320 55298 = 23 al 100.0 mel 77264 77286
55297 55275 = 23 al 95.7 pse 4870 4848
55293 55270 = 28 al 78.6 mel 77299 77326
51694 51682 = 13 al 100.0 pse 2602 2590
51394 51371 = 25 al 76.0 mel 81016 81040
50844 50812 = 33 al 87.9 mel 81456 81488
50844 50812 = 33 al 87.9 pse 1595 1563
50778 50759 = 20 al 90.0 mel 81500 81519
50778 50759 = 20 al 90.0 pse 1507 1488
50561 50534 = 28 al 89.3 mel 81610 81636
50561 50533 = 31 al 83.9 pse 1293 1263
50510 50484 = 28 al 78.6 mel 81672 81699
49361 49335 = 27 al 81.5 pse 889 866
35316 35291 = 26 al 80.8 mel 84411 84435 AADE01000322
29277 29253 = 25 al 88.0 pse 884 860
27772 27750 = 23 al 78.3 mel 85276 85298
26246 26222 = 25 al 76.0 pse 738 762
25320 25295 = 26 al 80.8 mel 90906 90931
22691 22654 = 38 al 76.3 pse 1026 1063
22568 22545 = 26 al 80.8 pse 1257 1282
22123 22091 = 33 al 78.8 pse 1341 1373
19528 19502 = 28 al 75.0 pse 8256 8283
19390 19368 = 23 al 87.0 pse 8395 8417
104
RESULTADOS
105
RESULTADOS
CG1288-CG14609-CG2520
D. buzzatii 40C11 D.melanogaster D. pseudoobscura
AE003672 AADE000014
inicio fin tamaño % inicio fin inicio fin
1087 1167 = 82 al 86.6 mel 125136 125056
1095 1156 = 62 al 93.5 pse 204478 204539
1228 1249 = 22 al 95.5 mel 125042 125021
1233 1266 = 34 al 82.4 pse 204564 204594
1288 1306 = 19 al 100.0 mel 124977 124959
1303 1368 = 71 al 81.7 pse 204680 204749
1317 1377 = 67 al 80.6 mel 124936 124870
2386 2416 = 31 al 93.5 mel 124735 124706
2673 2702 = 30 al 83.3 pse 205258 205286
2854 2886 = 33 al 84.8 pse 205482 205514
2857 2886 = 30 al 90.0 mel 123983 123954
2982 3022 = 42 al 85.7 mel 123858 123817
2983 3022 = 43 al 79.1 pse 205672 205714
3054 3079 = 26 al 96.2 mel 123398 123373
3054 3081 = 28 al 96.4 pse 205765 205792
3170 3235 = 68 al 91.2 pse 205816 205883
3175 3268 = 103 al 81.6 mel 123342 123241
3236 3258 = 23 al 91.3 pse 205893 205915
3297 3308 = 12 al 100.0 pse 205965 205976
4378 4398 = 21 al 85.7 pse 206342 206361
4388 4411 = 26 al 84.6 mel 122681 122656
4488 4529 = 42 al 78.6 pse 206379 206417
4499 4522 = 24 al 91.7 mel 122622 122599
4720 4779 = 60 al 81.7 mel 122404 122345
4731 4781 = 51 al 90.2 pse 206679 206729
4792 4818 = 28 al 78.6 pse 206739 206764
8582 8609 = 28 al 78.6 pse 208343 208367
8626 8661 = 36 al 83.3 pse 208396 208431
8628 8672 = 45 al 88.9 mel 120781 120737
8704 8742 = 39 al 89.7 mel 120713 120675
8708 8736 = 29 al 96.6 pse 208466 208494
8965 8989 = 25 al 80.0 mel 120566 120542
9021 9048 = 28 al 82.1 mel 120518 120495
9024 9048 = 25 al 88.0 pse 208671 208695
10255 10286 = 33 al 90.9 mel 120408 120376
106
RESULTADOS
107
RESULTADOS
108
RESULTADOS
109
RESULTADOS
expression
110
RESULTADOS
mel actgttcagt ctcataaaaa atgcagcgga aactagaaaa aatcacaatg attta----c aactgcagcc gagcaacata tacaattctg cacctttca-
pse ....ccg.a. a......... .....ta..t ttga----.. .cc.ga.... .....----. .g.....a.t ...a.cgca. cgg.gaatc. ..aga.atgc
buz ---------- ---------- ---------- ---------- ...t.tg..t ta.g.tatg. ....a....g .----.a..c .g..caaa.. .........-
mel --------ta atattacgta aaaattgtcc cttttttatt ttgtctgctt gaacagt--- -taccgttgg aaagtcgcaa cttgca-aac cataaaaaat
pse acctttcac. g.......c. .......... t......... ----.g...c ...tgt.gtt g..gt.--.. .......... ....a.-... ..........
buz --------a. ........c. ...tc..... ...gca.-.g .....g..-- ---------- ---------- -.ct...... .....cc... ..........
mel –tgccgttga catgaaaacc ttttgccag- -aaaactcgg cacaaaggca gcagtgcaaa aaaga----- --agagacat gaaa------ -gccataaaa
pse –...t.cc.. ..c......a ....a.ag.a a........a .......ca. .a.ca.ag.g .g...ggcac ag........ a...------ -.........
buz a.....c... .g.a...... ...ataa.a- -..------- ---------. -------... ...a.----- --c..a.a.. tg..aggtgc t.t.......
mel cggcaagcgg gtat----aa gtga-----g aggctgttag gctgggact- ---gggacag caatttgtag atagcgaatt tttgatttcc gtaacataaa
pse .......... tcgaccca.. ....-----. .....c.g.t c.----.gag gat.c.g... .......... ........a. .......... .a........
buz .a........ .-.a----.. c...tgcct. .cca..ca.t tt.----.-- ---------- --------.. .....t..a. .a........ aa........
mel cacggagcct gtgccattcc g--------- ---------- ---------- -agcgcatct –ccgcacacc gtcgcgcaaa aaattcgcac aggtagcatg
pse .......... ...------- ---------- ---------- ---------- ---------- ------.... .......... .......... .......g..
buz ....a.a.a. aaa.gca.t. ccatttctac tattctgctg caattccctt tgctg..... gc..c.c... a.....t... .......t.. .......g..
mel gaaaattgct aaaaataaga aacaaaaa-- ---------- -actagacaa g----aaagt agagaaagga atgcggg--- gaaaatcggg gaaaag-cgc
pse ......c.g. .....a.ga. ..t..c.c-- ---------- -.ag.at.gc ctggc.t..g ..c.t.g.a. .---ac.--- .....c..a. tc....a.a.
buz c.g....taa .....a..a. .ga.g...ta aaaacagcca g..a.at... c----....a .a.tc...a. .a.t.a.aca ....gg.at. .....c-gtt
mel tgg------c tgctggccac cagggcgaga cctgccgctg ggcgagtcgg tgtcggtatt ggta------ --gcttcacc Cttt-cacgc actgcgtgcg
pse a..gccata. ......g... .gc....... ........a. c.g...c... ...------- -..t------ --t....... ....-..t.. .........a
buz .c.------. .------..a .cca...g.. gt..gt...a caaa..ct.c .tct.c.tc. .c.tctgctt ct........ ....a..t.. cg...agt..
Figure S1. Detail of two CRM. Fragments of the AVID alignments (mel/pse and mel/buz) used
to generate the VISTA plots of Figure 3. (A) lab550-HOMRE (Homeotic Response Element):
The four binding sites of this CRM are conserved between all three species in one CNS. (B)
iab2(1.7) (fragment): All species show five Hunchback (HB) binding sites (three are conserved
between all three species; one is shared by D. melanogaster and D. pseudoobscura; another one
is shared by D. pseudoobscura and D. buzzatii; and there is one unique in D. melanogaster and
another one in D. buzzatii); a unique Krüppel (KR) binding site, where point mutations (red
uppercase) in D. melanogaster cause gain of expression mutants (Hab1 and Hab2) (Shimell et
al. 1994), is conserved in all three species. This region contains 9 CNS in the D. melanogaster-
D. pseudoobscura comparison and 6 CNS between D. melanogaster-D. buzzatii. PATCH
(http://www.gene-regulation.com/cgi-bin/pub/programs/patch/bin/patch.cgi) was used for
Hunchback binding site prediction in the iab2(1.7) enhancer.
111
RESULTADOS
Figure S2. Expression pattern of pb in imaginal discs. Labial disc: pb is expressed in the
labial disc of late third instar larvae and prepupae of the four tested species. Expression is
detected over the whole disc, which will give rise to the proboscis. Eye-antennal disc: pb is also
expressed in the presumptive palpus (arrowhead) of the eye-antennal disc of prepupa and white
pupa, which will give rise to the maxillary palps. In overstained discs some pb expression can
be detected in the presumptive ocellus. Some staining has also been observed in the lateral side
of the stalk in white pupae eye-antennal discs of D. repleta. However, it has not been
determined whether this expression domain is shared by the other species, as this region is
easily broken.
Figure S3. Expression pattern of lab in embryos. Cephalic region: lab is expressed in a
ventral stripe just anterior to the cephalic furrow, which corresponds to the intercalary segment
(black arrow). This domain of expression was first detected at stage 6 in D. buzzatii and at stage
8 in the other species. Later on lab is expressed in the dorsal ridge (white arrow) (stage 11-13)
and in sensory structures derived from the intercalary segment (boxed area). Midgut: at stage 8
lab expression is detected weakly in the anterior and posterior midgut primordial (arrowheads).
When the two primordia fuse (stage 12-13) lab expression in this region suddenly increases.
Expression is kept at high levels in the second midgut constriction and finally in the second
midgut chamber (stage 17).
112
RESULTADOS
Figure S4. Expression pattern of lab in the eye-antennal disc. Expression of lab was detected
in the peripodial membrane of the eye-antennal disc from late third instar larvae and prepupae.
This region will give rise to the lateral and posterior head capsule. No expression has been
detected in other imaginal discs.
Figure S6. Expression pattern of abdA in imaginal discs. Female genital disc: ABD-A
protein is observed in the presumptive internal genitalia region (derived from the A8 segment)
of the female genital disc. Proper staining with anti-ABD-A was only obtained for D.
melanogaster and D. buzzatii. All imaginal discs were tested for abdA expression. Central
nervous system of the larva: ABD-A protein is always detected in seven stripes in the larval
central nervous system of the four species, although the morphology of the larval central
nervous system is slightly different.
113
RESULTADOS
114
3. DISCUSIÓN
DISCUSIÓN
117
DISCUSIÓN
118
DISCUSIÓN
contener agrupaciones de CNS, donde cada uno contiene uno o varios lugares de unión.
Habría que estudiar más casos para ver cuán general es el patrón observado, ya que la
muestra es excesivamente pequeña. En el caso del módulo de la banda 2 de even-
skipped (eve S2E – even-skipped Stripe 2 Element) existen múltiples lugares de unión
para cada uno de los factores de transcripción involucrados (LUDWIG et al. 2000) lo que
concuerda con la observación anterior de elevada variación en la secuencia de algunos
lugares de unión.
De todos modos, esta dinámica y variabilidad en los lugares de unión no será
siempre neutra. Aunque la conservación de las secuencias no codificantes esté
demostrando ser un buén método para la detección de los módulos reguladores, no hay
que olvidar que la variabilidad fenotípica se debe en muchos casos a la variabilidad de
estas regiones (DE MEAUX et al. 2005). De manera que el uso de la conservación para la
detección de las regiones debería ir seguido de la disección de la variabilidad
nucleotídica y su posible relación con la variabilidad en la expresión génica y en último
término la variabilidad fenotípica. Debemos centrarnos en los cambios funcionales en
vez de en la conservación para entender la dinamica evolutiva a corto plazo de las
regiones cis-reguladoras.
119
DISCUSIÓN
anteroposterior (BERLETH et al. 1988), zen es un gen zigótico que se expresa en el tejido
extraembrionario del embrión, la amnioserosa (RUSHLOW et al. 1987), y zen2 es un gen
con patrón de expresión similar a zen pero sin función conocida ya que la ausencia de
este gen no presenta alteraciones en el fenotipo (RUSHLOW et al. 1987). Los genes bcd y
zen no son miembros del sistema Hox de determinación del patrón anteroposterior sino
que representan genes con homeodominio de evolución rápida (DE ROSA et al. 1999).
La clonación de los genes de tipo Hox3 en S. gregaria y T. castaneum mostró la
homología de este gen con zen. Mientras el homeodominio de estos genes era
claramente de la familia Hox3, tanto motivos de secuencia fuera del homeobox como su
patrón de expresión en los tejidos extraembrionarios mostraban que eran homólogos del
gen zen de D. melanogaster (FALCIANI et al. 1996). La descripción, en el ácaro
Archegozetes longisetosus, de un gen homólogo a zen pero con patrón de expresión tipo
Hox confirmó que en realidad zen era el gen Hox3 modificado (TELFORD y THOMAS
1998). En este organismo Hox3 presenta un patrón de expresión ampliamente solapado
con pb, coincidiendo en su límite anterior, lo que llevó a la hipótesis de que deben
presentar una cierta redundancia funcional. Esta redundancia entre pb y Hox3 habría
sido aprovechada por los insectos para utilizar Hox3 para una nueva función (FALCIANI
et al. 1996; TELFORD y THOMAS 1998). Un estudio reciente en el insecto áptero
Thermobia domestica muestra que Hox3 todavía retiene el patrón de expresión Hox,
pero ha extendido su expresión a los tejidos extraembrionarios (HUGHES et al. 2004). En
el ancestro de los insectos alados, como S. gregaria o T. castaneum, este gen perdió su
dominio de expresión Hox y reforzó su expresión en los tejidos extraembrionarios
(FALCIANI et al. 1996). El gen Hox3 parece haber evolucionado directamente a zen
durante la diversificación de los insectos (HUGHES et al. 2004).
Posteriormente, durante la diversificación de los dípteros, zen adquirió un nuevo
dominio de expresión que provee un aporte materno en el embrión (STAUBER et al.
2002). Una duplicación de este gen en la línea de las moscas ciclorrafas permitió la
separación de ambos dominios de expresión, dando lugar a dos genes con funciones
diferenciadas: el morfógeno materno bcd y el gen zigótico zen (STAUBER et al. 1999;
STAUBER et al. 2002).
Una segunda duplicación de zen originó el gen zen2. Un análisis comparativo de
la estructura del ANT-C entre D. pseudoobscura y D. melanogaster fue incapaz de
detectar un homólogo de zen2 ni en la región estudiada ni en el resto del genoma, de
manera que se pensó que zen2 era una duplicación muy reciente en el linaje de D.
120
DISCUSIÓN
121
DISCUSIÓN
En las tres especies de Drosophila analizadas, el gen zen2 está situado dentro de
las regiones reguladoras de pb, aunque en orientaciones diferentes. Esta situación se
podría interpretar como que la orientación ancestral es la común a D. pseudoobscura y
D. buzzatii, y que el gen se ha invertido en D. melanogaster. Pero el hecho que las
duplicaciones en tándem suelen dar lugar a dos genes situados en la misma orientación
y la posición de zen2 situado entre el promotor basal de pb y sus regiones cis-
reguladoras, apuntan a una evolución más compleja.
Después de las dos duplicaciones que dieron lugar a zen, bcd y zen2 a partir del
gen Hox3 tuvieron lugar dos inversiones (Figura 15). Una primera inversión con un
punto de rotura entre el promotor basal de pb y sus regiones reguladoras 5’ y el otro
entre zen y zen2 se habría producido en el ancestro de los subgéneros Drosophila y
Sophophora, dando lugar a la estructura que observamos en D. pseudoobscura y D.
buzzatii. Una segunda inversión, en el linaje de D. melanogaster, afectó únicamente al
gen zen2. Esta segunda inversión habría reutilizado el punto de rotura entre zen2 y el
promotor basal de pb, y tendría el segundo punto de rotura entre zen2 y la región
reguladora 5’ de pb. Este modelo explica, mediante solo dos inversiones, la orientación
del gen zen2 y la posición de las regiones reguladoras de pb en las tres especies donde
se ha analizado la región.
3.2.1.2. lab
labial (lab) es el gen más proximal del ANT-C de D. melanogaster (WAKIMOTO
y KAUFMAN 1981). Corresponde al primer grupo parálogo (PG1). Los mutantes para
122
DISCUSIÓN
123
DISCUSIÓN
del embrión en desarrollo en las diferentes especies analizadas muestra que esta función
fue adquirida por lab por lo menos antes de la diversificación del género Drosophila.
La comparación de la secuencia codificante de lab realizada entre tres especies
de Drosophila (ver Artículo 1) muestra una elevada tasa de substitución nucleotídica. Se
contrastó la posibilidad de una variación en la tasa de substitución nucleotídica como
consecuencia del cambio de posición. Los resultados mostraron una tasa homogénea a
lo largo de la divergencia del género, aunque heterogénea a lo largo de la región
codificante. La región más conservada corresponde al homeobox, pudiendose observar
la presencia de otras pequeñas regiones que se encuentran también conservadas, a nivel
de proteína, en especies alejadas como T. castaneum o incluso en Humanos. Esta
conservación a nivel de dominio proteicos explica como minigenes de Hoxb-1 de pollo,
expresados bajo el control de regiones reguladoras de lab de D. melanogaster, eran
capaces de rescatar mutantes nulos para lab (LUTZ et al. 1996). Aunque para algunas
funciones del gen los cambios en la proteína pueden ser significativos, para otras
funciones es más importante el momento y lugar de expresión, es decir, las regiones cis-
reguladoras.
3.2.1.3. pb
proboscipedia (pb) es el segundo gen del HOM-C de Drosophila. El primer
mutante de pb fue descrito en 1933 por Bridges y Dobzhansky. Las mutaciones nulas de
pb producen transformación homeótica específica del adulto, donde los palpos labiales
124
DISCUSIÓN
125
DISCUSIÓN
Figura 16. Expresión del gen pb en embriones (estadío 11) de cinco especies de Drosophila. La
flecha negra muestra la expresión en el mesodermo ventral del segmento mandibular, la flecha
blanca señala la expresión en el mesodermo ventral del segmento maxilar. El segmento
mandibular (mn), maxilar (mx) y labial (lab) se muestran en el primer recuadro.
Figure 16. Expression of pb in embryos (stage 11) from five Drosophila species. The black
arrow shows the expression in the ventral mesoderm of the mandibular segment, the white
arrow shows the expression in the ventral mesoderm of the maxillary segment. The mandibular
(mn), maxilar (mx) and labial (lab) segments are shown on the first embryo.
126
DISCUSIÓN
3.2.1.4. abdA
Cuando Edward Lewis (1978) describió el complejo Bithorax, indicó la
presencia de al menos ocho genes, responsable cada uno de dar identidad a un
segmento. Posteriormente, el aislamiento y caracterización de mutaciones letales en
heterocigosis con defiencias del BX-C mostraron la existencia de solo tres grupos de
complementación (SÁNCHEZ-HERRERO et al. 1985; TIONG et al. 1985). Estudios
moleculares confirmaron la existencia de solo tres genes codificantes: Ubx, abdA y
AbdB (CELNIKER et al. 1989; KARCH et al. 1990). Los “genes” identificados por LEWIS
(1978) eran en realidad regiones cis-reguladoras de estos tres genes. Cada una de estas
regiones cis-reguladoras se encarga de activar y mantener la expresión de uno de estos
genes en un segmento (o parasegmento) determinado, y están situadas de forma
colineal.
127
DISCUSIÓN
De los genes Hox estudiados en este trabajo, abdA parece ser el que presenta
más constricciones funcionales. La región abdA no presenta ninguna microinversión, a
diferencia de lab y pb, ni inserción de elementos móviles, los cuales son muy
abundantes en el gen lab situado a tan solo 75 kb (Artículo 2 Figura 2). Además los
elementos reguladores de abdA, a diferencia de lab y pb, presentan colinealidad.
128
DISCUSIÓN
129
DISCUSIÓN
citológico en solo dos regiones genómicas, ya que el gen lab está a tan solo 75 kb de
abdA.
130
DISCUSIÓN
Figura 17. Modelo de evolución del HOM-C en Drosophila. Arriba, filogenia de las especies de
Drosophila estudiadas en este trabajo. Se muestra el lugar de la filogenia en que tuvieron lugar
los diferentes eventos. El recuadro con el número (1) muestra los eventos sucedidos antes de la
diversificación del género. Los otros dos recuadros muestran el estado del complejo antes de la
diversificación, punto de partida donde se produjeron el resto de reorganizaciones. Símbolos y
colores igual que en la Figura 15.
Figure 17. Model of the evolution of the HOM-C in Drosophila. Up, phylogeny of the
Drosophila species studied in this thesis. The numbers show where in the phylogeny each event
took place. The rectangle with number (1) shows those events that happened before the
diversification of the Drosophila genus. The other rectangles show the structure of the complex
before the diversification, starting point where the rearrangements took place. Symbol and
colours as in Figure 15.
131
DISCUSIÓN
132
DISCUSIÓN
133
DISCUSIÓN
134
DISCUSIÓN
Hay que destacar que a pesar de la importancia del complejo de genes Hox, en
Drosophila es posible romper este complejo y tener una mosca viable. Trabajos de
roturas inducidas o rescates por minigenes parecían indicar la posibilidad de estas
roturas en condiciones de laboratorio, pero cuestionaban si serían capaces de sobrevivir
a las presiones selectivas de la naturaleza (RANDAZZO et al. 1991). Hay que recalcar que
las roturas del HOM-C descritas en el apartado 3.2.2.1 son naturales. No se trata de
roturas experimentales cuyos portadores sobreviven en el laboratorio sin aparentes
malformaciones, sino que se trata de inversiones fijadas en grupos de especies con
elevado éxito evolutivo. Incluso la rotura más reciente, entre lab y pb, está fijada en
todo el grupo repleta, uno de los grupos más exitosos con más de 100 especies.
135
DISCUSIÓN
136
DISCUSIÓN
137
DISCUSIÓN
138
DISCUSIÓN
139
DISCUSIÓN
Figura 18. Estructura del HOM-C en los Metazoos. Árbol filogenético con la información
disponible sobre la organización genómica de los genes Hox. Información obtenida de: S.
gregaria (FERRIER y AKAM 1996); T. castaneum (DENELL et al. 1996; BROWN et al. 2002); A.
gambiae (POWERS et al. 2000; LEWIS et al. 2003); B. mori (YASUKOCHI et al. 2004); C. elegans
(DE ROSA et al. 1999); D. rerio (AMORES et al. 2004); H. francisci de (KIM et al. 2000; CHIU et
al. 2002); B. floridae (MINGUILLÓN et al. 2004); C. intestinalis (SPAGNUOLO et al. 2003; IKUTA
et al. 2004); O. Dioica (SEO et al. 2004), y S. purpuratus (MARTINEZ et al. 1999).
Figure 18. HOM-C structure in the Metazoan. Phylogenetic tree with the available information
on the genomic organization of Hox genes. Information from: S. gregaria (FERRIER y AKAM
1996); T. castaneum (DENELL et al. 1996; BROWN et al. 2002); A. gambiae (POWERS et al.
2000; LEWIS et al. 2003); B. mori (YASUKOCHI et al. 2004); C. elegans (DE ROSA et al. 1999); D.
rerio (AMORES et al. 2004); H. francisci (KIM et al. 2000; CHIU et al. 2002); B. floridae
(MINGUILLÓN et al. 2004); C. intestinalis (SPAGNUOLO et al. 2003; IKUTA et al. 2004); O. Dioica
(SEO et al. 2004), y S. purpuratus (MARTINEZ et al. 1999).
140
DISCUSIÓN
que S. gregaria contiene extensos elementos reguladores específicos de cada gen, más
que elementos compactos y sobrepuestos como en vertebrados. (FERRIER y AKAM 1996)
El escarabajo T. castaneum también posee un HOM-C (BEEMAN 1987). En
algunos aspectos, este complejo es más parecido al HOM-C de vertebrados que al de
Drosophila, por ejemplo en que incluye unidades de transcripción grandes y distancias
intergénicas largas importantes para la elaborada regulación de los genes (DENELL et al.
1996). El HOM-C de Tribolium forma un bloque contiguo y no presenta genes
accesorios (MCGINNIS y KRUMLAUF 1992). Se ha clonado y secuenciado la región
equivalente al ANT-C (BROWN et al. 2002).
En el 2000 dos trabajos estudiaron la organización genómica de los genes Hox
en A. gambiae (DEVENPORT et al. 2000; POWERS et al. 2000). Posteriormente, la
secuencia genómica ha completado los datos restantes confirmando la presencia de un
complejo único y entero. Las transposiciones en el ancestro de Drosophila, así como las
duplicaciones génicas (Hox3), se produjeron después de la separación de los linajes de
Drosophila y Anopheles (LEWIS et al. 2003).
141
DISCUSIÓN
142
4. CONCLUSIONES
CONCLUSIONES
145
CONCLUSIONES
8. Los patrones de expresión de los genes Hox lab y abdA están conservados
entre D. melanogaster, D. virilis, D. repleta y D. buzzatii. El gen pb presenta un
dominio extra de expresión en los embriones de D. virilis, pero esta diferencia no parece
estar relacionada con la rotura lab-pb. Por lo tanto las roturas del HOM-C no han tenido
efecto sobre el patrón de expresión de los genes Hox adyacentes (lab, pb o abdA).
146
CONCLUSIONS
1. The nucleotide sequences from the regions lab-abdA and pb of D. buzzatii and
their comparison with the sequences of D. melanogaster and D. pseudoobcura shows
that the split of the HOM-C between the genes lab and pb was produced at less than 3
kb from the 3’ end of pb, between this gene and the cluster of cuticular genes (Ccp).
2. In D. buzzatii the genes lab and abdA, located at the same cytological
position, are at 75 kb only. In the intervening region there is only one gene, CG31363,
which is not related with the Hox genes. Despite their proximity, the position of these
two genes seems fortuitous and without functional meaning.
3. The sequences from the lab and pb regions of D. buzzatii and their
comparison with those of D. melanogaster and D. pseudoobscura show that the lab-pb
split was produced through a paracentric inversion and not by transposition. The
second breakpoint was located between the CG31363 and CG17836 genes. However,
there is no firm evidence about the mechanism that produced this inversion.
4. The comparative analysis of the coding regions of the lab gene between D.
buzzatii, D. virilis and D. melanogaster shows that the evolution rate is heterogeneous
along the gene, but homogeneous along the phylogeny. Thus, the change in position had
no effect on the nucleotide substitution rate of lab.
6. The comparison between the conserved blocks and the known regulatory
regions in D. melanogaster shows that the position and orientation of the regulatory
regions is conserved between the three species, except small inverted regions that
probably correspond to individual CRM.
147
CONCLUSIONS
7. The three HOM-C splits known in the Drosophila genus have been produced
between the regulatory regions of the adjacent genes. In all three cases the splits took
place near the 3’ end of one gene and far from the 5’ end of the other gene.
8. The expression patterns of the Hox genes lab and abdA are conserved
between D. melanogaster, D. virilis, D. repleta and D. buzzatii. The gene pb presents an
extra domain of expression in the embryos of D. virilis, but this difference seems
unrelated to the lab-pb split. Thus, the HOM-C splits had no effect on the expression
pattern of the adjacents Hox genes (lab, pb or abdA).
10. Besides Drosophila, four species are known to have a split HOM-C: the
nematode C. elegans, the lepidopteran B. mori and the tunicates C. intestinalis and O.
dioica. These organisms seem to keep spatial but not temporal colinearity. The presence
of a derived and quick mode of embryogenesis could be the ultimate cause of the loss of
temporal colinearity and the HOM-C disassembly.
148
5. REFERENCIAS
REFERENCIAS
151
REFERENCIAS
152
REFERENCIAS
153
REFERENCIAS
154
REFERENCIAS
155
REFERENCIAS
156
REFERENCIAS
157
REFERENCIAS
158
REFERENCIAS
159
REFERENCIAS
160
REFERENCIAS
161
REFERENCIAS
162
REFERENCIAS
163
REFERENCIAS
164
REFERENCIAS
165
REFERENCIAS
166
REFERENCIAS
167
REFERENCIAS
168
REFERENCIAS
169
REFERENCIAS
170
REFERENCIAS
171
REFERENCIAS
172
REFERENCIAS
173
REFERENCIAS
174
ANEXO 1. PROTOCOLOS
ANEXO 1
PROTOCOLOS
1. PROCEDIMIENTOS BÁSICOS 181
1.1. Medios de Cultivo 181
1.2. Soluciones Estoc 182
1.3. Electroforesis en Gel de Agarosa 183
1.4. Digestión con Enzimas de Restricción 184
177
ANEXO 1
178
ANEXO 1
179
ANEXO 1
180
ANEXO 1
1. PROCEDIMIENTOS BÁSICOS
TB (Terrific Broth)
Preparación 1,2 g triptona 0,231g KH2PO4
2,4 g extracto de levadura 1,254 g K2HPO4
0,4 mL glicerol 10 mL agua
90 mL agua destilada Autoclavar
Autoclavar
Dejar enfriar y mezclar las dos soluciones
X-gal 20 mg/mL
Preparación 200 mg 5-bromo-4-cloro-3-indolil-β-galactopiranosido
10 mL dimetilformamida (en tubo de polipropileno)
Hacer alícuotas de 1 mL
Guardar a –20ºC.
Usar 2 µL por mL de cultivo
Ampicilina 25 mg/mL
Preparación 250 mg Ampicilina
10 mL agua destilada
Esterilizar por filtración
Guardar a 4ºC.
Usar 2 µL por mL de cultivo
181
ANEXO 1
EDTA 0,5 M
Preparación 186,1 g etilendiamintetraacetatox2H2O
800 mL agua destilada
Añadir unos 20 g de NaOH
Ajustar a pH 8
Enrasar a 1 L, Autoclavar
TrisHCl 1M
Preparación 121,1 g Tris
800 mL agua destilada
Añadir unos 65 mL de HCl 35%
Ajustar al pH deseado (7, 7,5 ó 8)
Enrasar a 1 L, Autoclavar
SDS 20%
Preparación 20 g SDS
100 mL agua destilada
Calentar a 68ºC
182
ANEXO 1
183
ANEXO 1
Tiempo de incubación:
DNA genómico Toda la noche
DNA plasmídicos De 2 horas a toda la noche
DNA fagos lambda 1-2 horas
DNA BACs 3 horas
184
ANEXO 1
2.1. Puesta
1. Preparar los medios de puesta.
(Las placas con el medio se pueden guardar a 4ºC varios días).
2. Precalentar las placas y poner la levadura justo antes de usar.
3. Poner la placa de medio con levadura en la Caja o Falcon de puesta (Figura A1).
Dejar a la temperatura óptima de la cepa.
4. Cambiar los medios de puesta después del tiempo deseado.
5. Los embriones se pueden recolectar inmediatamente o dejar las placas hasta un par
de días en la nevera (adecuado cuando las puestas contienen muy pocos embriones).
# Es conveniente empezar siempre con dos puestas de 1 hora, descartar estas puestas.
Medios de puesta:
D. melanogaster Agar al 2% con zumo de manzana, y una pizca de levadura.
D. pseudoobscura
D. virilis Agar con zumo de manzana y una capa fina de levadura (para
evitar que los huevos queden demasiado undidos en el medio).
D. buzzatii Agar al 2% con ácido acético y una capa fina de levadura
D. repleta (sólo ponen embriones sobre la levadura).
185
ANEXO 1
bjhjkhfkfuk
2.3. Decorionización
Decorionizar los embriones mediante inmersión en lejía.
1. Introducir el filtro con los embriones en un vaso de precipitados
pequeño con lejía 50%.
2. Mantener entre 2 y 5 minutos.
3. Controlar el proceso bajo la lupa.
4. Lavar con agua del grifo o con solución de lavado.
5. Secar ligeramente.
186
ANEXO 1
2.4. Almacenaje
Para extracción de RNA:
1. Colocar el filtro boca abajo sobre una placa de Petri.
2. Hacer caer los embriones en la placa de petri con solución de lavado.
3. Rotar lentamente la placa de forma que los embriones se acumulen en el centro.
4. Recoger los embriones con una pipeta Pasteur y pasarlos a un eppendorf
(poner la pipeta vertical para que los embriones no se enganchen a la pared).
5. Centrifugar brevemente y eliminar el líquido.
6. Guardar los embriones a –20ºC.
Para hibridación in situ o tinción con anticuerpo fijar según protocolo (ver
Apartados 9.2 y 11.1.1).
187
ANEXO 1
188
ANEXO 1
Solución III
Preparación 11,5 mL ácido acético glacial
60 mL acetato potásico
28,5 mL agua destilada
189
ANEXO 1
190
ANEXO 1
191
ANEXO 1
192
ANEXO 1
193
ANEXO 1
194
ANEXO 1
195
ANEXO 1
196
ANEXO 1
197
ANEXO 1
198
ANEXO 1
199
ANEXO 1
5.2. Ligación
En pBSK: 1. Poner en un eppendorf: 7 µL DNA inserto
1 µL DNA vector linealizado
2. Incubar 5 minutos a 65ºC.
3. 10 minutos en hielo
4. Añadir: 1 µL de tampón de Ligasa 10x
1 µL de T4 DNA ligasa (1 unidad/ µL)
5. Incubar toda la noche a 14-16ªC,
ó 3-4 horas a temperatura ambiente.
En pGEM-T o pGEM-Teasy:
Tampón 10x 1. Poner en un eppendorf:
7 µL DNA inserto
1 µL vector
1 µL tampón ligasa 10x
1 µL ligasa
2. Incubar a 4ºC durante toda la noche.
200
ANEXO 1
201
ANEXO 1
5.4. Transformación
La tranformación se ha realizado por el método de choque térmico.
Rutinariament se han utilitzado 5 µL de ligación para transformar una alíquota de
células competentes, pero en casos en que es partía de poco DNA o de mala calidad se
utilizan los 10 µL.
202
ANEXO 1
Preparación de la Marcaje de
muestra la sonda
Desnaturalización de la sonda
y de la muestra
Hibridación
Lavados
Revelado
Visualización
203
ANEXO 1
7.2. Transferencia
Se utiliza el protocolo de transferencia por capilaridad a membrana de nylon en
condicions neutras descrito en Sambrook et al. (1989).
204
ANEXO 1
pes (500-1000 g)
placa de vidre
paper absorbente
3 papers Whatman 3MM
FILTRE gel d’electroforesis
marc impermeable
paper Whatman 3MM
10 x SSC
205
ANEXO 1
206
ANEXO 1
7.4. Hibridación
7.4.1. Prehibridación e Hibridación
1. Incubar el filtro 1-2 horas en 10 mL de solución de hibridación,
a la temperatura de hibridación (42ºC para homólogas y 37ºC para heterólogas)
2. Mientras desnaturalizar la sonda: Hervir 5 minutos
Poner en hielo
3. Poner la sonda desnaturalizada en 10 mL de solución de hibridación
precalentada a temperatura de hibridación.
# Las sondas reutilizadas se desnaturalizan 10 minutos a 65ºC.
4. Descartar la solución de prehibridación y
añadir la solución de hibridación con sonda.
5. Incubar a la temperatura de hibridación durante mínimo 16 horas para DNA
genómico ó 12 horas para DNA de clones.
207
ANEXO 1
208
ANEXO 1
Tampón 3
Preparación 100 mL TrisHCl 1 M
5,844 g NaCl
10,17 g MgCl2
Ajustar a pH 9,5 con NaOH (26-30 lentejas)
Enrasar a 1 L con agua destilada
Solución de bloqueo
Preparación 100 mL Tampón 1
1 g leche en polvo
ó 90 mL Tampón 1
10 mL Reactivo de bloqueo 10%
Solución de anticuerpo
Preparación 20 mL Solución de bloqueo
4 µL anticuerpo (Antidigoxigenin-AP)
Solución de coloración
Preparación 10 mL Tampón 3
200 µL NBT/BCIP Stock Solution
209
ANEXO 1
210
ANEXO 1
6. Poner los portaobjetos en una cubeta con solución SSC Denhart en una baño a 65°C.
Dejarlos en esta solución entre 2 y 3 horas.
7. Poner los portaobjetos en una cubeta con etanol:ácido acético (3:1)
durante 20 minutos.
8. Sumergir los portaobjetos durante 2 segundos unas 10 veces en una cubeta
con etanol al 95%.
9. Dejar secar y guardarlos a 4°C.
SSC Denhart
Preparación 0,012 g polivinilpirrolidona
0,012 g ficoll
0,012 g albúmina de suero bovino
90 ml 20xSSC
600 ml H2O destilada
211
ANEXO 1
9. Picar sobre el cubreobjetos con una punta de pipeta para romper los núcleos y
separar bien los cromosomas.
10. Poner un pliegue de papel absorbente sobre el cubreobjetos y presionar con el
pulgar. Hay que evitar el desplazamiento del cubreobjetos para que los
cromosomas no se rompan.
11. Observar las preparaciones al microscopio, seleccionar las mejores y colocarlas
en una carpeta plana que se dejará en posición horizontal toda una noche.
12. Coger un portaobjetos y sumergirlo en nitrógeno líquido hasta que este deje
de burbujear.
13. Colocar el filo de una cuchilla entre una esquina del cubreobjetos y el portaobjetos
y tirar hacia arriba hasta que salte el cubreobjetos.
14. Colocar el portaobjetos en una cubeta con etanol al 95% un mínimo de 15 minutos.
15. Dejar secar los portaobjetos y mirar al microscopio para seleccionar las mejores
preparaciones que se guardarán a 4°C hasta su hibridación.
212
ANEXO 1
213
ANEXO 1
214
ANEXO 1
8.3. Hibridación
8.3.1. Prehibridación
1. Incubar las preparaciones cromosómicas 30 minutos en 2xSSC en un baño a 65°C.
2. Poner las preparaciones en una cubeta con etanol al 70% durante 5 minutos.
3. Pasar las preparaciones a otra cubeta de etanol 70% donde se dejarán otros 5 minutos.
4. Pasar las preparaciones a otra cubeta con etanol 95% durante 5 minutos.
5. Dejar secar las preparaciones cromosómicas (mínimo 15 minutos).
8.3.2. Hibridación
Se han realizado dos tipos de hibridaciones: homólogas, cuando el DNA de la
sonda proviene de la misma especie que el material donde se hibrida, y heterólogas,
cuando el DNA de la sonda proviene de una especie diferente.
215
ANEXO 1
8.3.3. Posthibridación
1. Poner las preparaciones cromosómicas en una cubeta con 2xSSC a 37°C,
incubar durante 10 minutos.
2. Incubar otros 10 minutos en una segunda con 2xSSC a 37°C.
3. Incubar 10 minutos en una cubeta de 2xSSC a temperatura ambiente.
Repetir la operación.
4. Poner las preparaciones en una cubeta con 1xPBS durante 5 minutos.
8.3.4. Revelado
1. Añadir a 210 µL de solución I 3,5 µL de reactivo A y
3,5 µL de reactivo B del kit Elite.
2. Añadir a cada preparación 20 µL de la solución anterior.
3. Poner un cubreobjetos sobre la zona donde se ha añadido la solución.
Si quedan burbujas eliminarlas con la punta de la pipeta.
4. Incubar las preparaciones a 37°C durante 45 minutos en una cámara húmeda.
5. Poner las preparaciones en 1xPBS durante 10 minutos.
Repetir dos veces.
6. Añadir a cada preparación 500 µL de solución III.
7. Incubar las preparaciones a 37°C durante 45 minutos en una cámara húmeda.
8. Poner las preparaciones en 1xPBS de 5 minutos a 1 hora.
9. Preparar la solución de tinción añadiendo 5 mL de Giemsa a 100 mL de tampón P.
10. Poner las preparaciones en la solución de tinción durante 1 minuto.
11. Poner las preparaciones en una cubeta con H2O destilada.
Repetir en otra cubeta de H2O destilada.
12. Dejar secar las preparaciones.
13. Añadir a cada preparación una gota de montador de preparaciones biológicas
y poner un cubreobjetos.
216
ANEXO 1
solución III
Preparación 5 µg de diaminobencidina (DAB)
1 mL de Tris HCl 50 mM pH 7,6
1 µL de H2O2
tampón P
Preparación 40% NaH2PO4H2O 50 mM
60% Na2HPO4 12 H2O 50 mM
Ajustar a pH 7 con las soluciones anteriores
217
ANEXO 1
218
ANEXO 1
219
ANEXO 1
PAF 4%
Preparación 10 mL 10xPBS
70 mL H2O destilada
Calentar a 55ºC.
Añadir 4 g de paraformaldehído.
(Al disolverse queda transparente)
Poner 1 gota de NaOH 1M.
Ajustar el pH a 7,3 a Temperatura ambiente.
Enrasar a 100 mL con agua destilada.
Alicuotar y guardar a –20ºC.
paraformaldehído Fluka-Chemika 76240 (CH2O)n Mr 30.03n
9.3. Hibridación
1. Rehidratar los embriones con baños de 5 minutos em 1 mL de metanol/PBT:
70:50, 50:50, 25:75, 10:90.
2. Lavar 3 veces con 1 mL de PBT 5 minutos.
3. Incubar 30 minutos en i mL de PBT.
4. Reemplazar 0,5 mL de PBT con 0,5 mL de solución de hibridación.
5. Incubar 10 minutos.
6. Prehibridar en 1 mL de solución de hibridación 1 hora a 60ºC.
7. Mientras desnaturalizar la sonda: Hervir 10 minutos
Poner en hielo
8. Poner la sonda desnaturalizada en 300 µL de solución de hibridación
precalentada a temperatura de hibridación.
# Las sondas reutilizadas se desnaturalizan 10 minutos a 80ºC.
9. Descartar la solución de prehibridación y añadir la solución de hibridación con sonda.
10. Incubar toda la noche a 60ºC.
220
ANEXO 1
Solución de Hibridación
Preparación 50 mL Formamida
25 mL 20x SSC
200 µL tRNA (50 mg/ml)
1 mL Heparina (50 mg/ml)
1 mL Tween 10%
23 mL H2O destilada
9.4. Lavados
1. Retirar la solución de hibridación con sonda.(Guardar la sonda a –20ºC).
2. Añadir 1,5 mL de solución de hibridación.
3. Incubar 1 hora a 60ºC.
4. Retirar 1 mL de soluciuón y 750 µL de PBT prescalentado a 60ºC.
5. Incubar 5 minutos a 60ºC.
6. Retirar el líquido. Añadir 1 mL de solución de hibridación/PBT 1/1 a 60ºC.
7. 15 minutos en agitación a temperatura ambiente.
8. Lavar 5 veces con PBT durante 5 minutos.
9.5. Revelado
1. Poner 100 µL de anti-DIG-AP 1/2000 en PBT.
2. Incubar 2 horas a temperatura ambiente o toda la noche a 4ºC.
3. Lavar 5 veces con PBT durante 5 minutos.
4. Lavar 3 veces con APB durante 5 minutos.
5. Colorear en un pocillo con 1 mL de APB + 10 µL NBT + 5 µL BCIP.
6. Parar la coloración con 3 lavados de 15 minutos con PBT.
221
ANEXO 1
9.6. Montaje
1. Lavar 2 veces con PBS.
2. Deshidratar los embriones con lavados sucesivos en etanol/PBS:
30:70, 50:50, 70:30 y etanol.
3. Lavar 3 veces con etanol 100%.
4. Montar con DURCUPAN
a. Poner una gota de medio DURCUPAN en el centro de un portaobjetos.
b. Colocar a cada lado medio cubreobjetos (del nº cero).
c. coger los embriones (con un poquito de etanol)
con una pipeta de 200 µL con la punta cortada.
d. poner los embriones sobre la gota de medio.
e. colocar un cubreobjetos de 24x, con cuidado de que no quede aire.
f. moviendo con cuidado el cubreobjetos podemos girar los embriones
para verlos en diferente orientación.
Guardar a –20ºC mientras queramos poder girar los embriones. Siempre en horizontal.
Para endurecer el medio dejar un par de días a temperatura ambiente.
222
ANEXO 1
1. Poner una placa sobre hielo para mantener el material frío durante la disección.
2. Poner 1xPBS en tres pocillos de cada placa.
3. Colocar las larvas del estadio deseado o las pupas a diseccionar en el primer pocillo
de la placa en hielo.
4. Coger una larva o pupa y situarla en el primer pocillo de la segunda placa o en un
porta para realizar la disección.
Bajo la lupa:
5. Eliminar los restos de medio adheridos a la larva o pupa y pasarla a otro pocillo.
6. Partir la larva o pupa por la mitad con dos pinzas finas (larvas)
o con tijeras de disección (pupas).
7. Coger cada mitad y dar la vuelta a la cutícula como si fuera un guante, colocando
una pinza en el extremo y moviendo la parte abierta de la cutícula con la otra pinza.
(si solo se quiere un tipo determinado de disco, la otra mitad se desecha directamente)
8. Eliminar con cuidado la grasa y las estructuras internas de la larva.
9. Los discos quedan unidos a la cutícula por las tráqueas, separarlos un poco para que
queden accesibles pero sin soltarlos de la cutícula.
10. Pasar las mitades invertidas y limpias al tercer pocillo.
11. Eliminar los restos adheridos y pasarlas a un pocillo de la placa en hielo.
La disección debe durar un máximo de 20 minutos, para evitar la degradación del RNA.
12. Fijar según el protocolo a usar (hibridación in situ o tinción con anticuerpo).
(Si se ha obtenido poco material, después de fijadas se pueden juntar dos
o más muestras)
223
ANEXO 1
10.4. Hibridación
1. Prehibridar 1 hora en solución de hibridación a 55º (en bloque).
2. Desnaturalizar la sonda 10 minutos a 80ºC.
(5 µL de sonda en 50 µL de solución de hibridación)
(Sondas reutilizadas suelen dar mejor resultado, se reduce el fondo.)
3. Mantener la sonda en hielo.
4. Descartar la solución de prehibridación.
5. Poner la solución de hibridación con sonda desnaturalizada.
Incubar toda la noche a 55ºC (en bloque).
10.5. Revelado
1. Precalentar (5 minutos) 3 mL de solución de hibridación por cada muestra.
2. Quitar la sonda de cada muestra, sin que quede completamente seco
(guardar la sonda a –20ºC para reutilizar).
3. Lavar 2 minutos a 55ºC (en el bloque) con 1 mL de solución de hibridación
precalentada.
4. Dos lavados de 15 minutos a 55ºC con 1 mL de solución de hibridación precalentada.
5. Lavar 15 minutos con 50% solución de hibridación:50% PBT en balancín.
(Quitar 500 µL de sol hibridación y añadir 500 µL PBT).
6. Cinco lavados de 10 minutos en PBT (o más lavados más cortos p.e. 7 x 5 minutos).
7. Incubar con anti-DIG 1/2000 1,5 horas a Temperatura ambiente.
Usar 5 µL de anti-DIG (1/50) preabsorvido en 200 µL de PBT.
8. Cuatro lavados de 20 minutos en PBT.
9. Preparar el tampón de coloración.
224
ANEXO 1
Solución de Hibridación
Preparación 50 mL Formamida
25 mL 20x SSC
200 µL tRNA (50 mg/ml)
1 mL Heparina (50 mg/ml)
1 mL Tween 10%
23 mL H2O destilada
tRNA 109525 ROCHE
Tampón de coloración
Preparación 500 µL NaCl 1M (o 125 µL 4M)
250 µL MgCl2 1M
200 µL tRNA (50 mg/ml)
500 µL Tris 1M pH 9,5)
50 µL Tween 10%
3,75 mL H2O destilada
225
ANEXO 1
Solución de coloración
Preparación 1 mL de tampón coloración
9 µL de NBT
7 µL de BCIP (llamado también XPhos)
Bajo la lupa:
6. Diseccionar los discos uno a uno y colocarlos sobre una gota de glicerol en un
portaobjetos nuevo.
(Se pueden levantar con las pinzas o agujas de disección o con una pipeta de 20 µL
con la punta cortada, en este segundo caso tener cuidado de lo coger otros restos)
7. Una vez diseccionada toda la muestra, colocar un cubreobjetos (con cuidado de que
no queden burbujas) y sellar con pintauñas (mejor si es transparente).
8. Guardar a 4ºC.
226
ANEXO 1
227
ANEXO 1
11.2.2. Anticuerpos
1. Incubar con anticuerpo primario a 4ºC toda la noche (o 4h a temperatura ambiente).
Siempre en balancín.
2. Lavar 4 veces 20 minutos con PBT/BSA.
3. Incubar con anticuerpo secundario 2 h a Temperatura ambiente.
4. Lavar 3 minutos 20 minutos con PBT/BSA.
Al empezar el segundo lavado preincubar el AB complex
300 µL PBT + 3 µL A + 3 µL B
30 minutos a temperatura ambiente
5. Incubar los discos con el AB complex 45 minutos.
228
ANEXO 1
Anticuerpos secundarios
Anti-rabbit Ig, biotinylated species-specific whole antibody (from
donkey) 2 ml. RPN 1004. Amersham Pharmacia Biotech.
Anti-mouse Ig, biotinylated species-specific whole antibody (from
sheep) 2 ml. RPN 1001. Amersham.
229
ANEXO 1
LB TopAgar
Preparación 1 g Tryptona
0,5 g Yeast extract
0,5 g NaCl
0,6 g Agarosa
Enrasar a 100 mL con agua destilada
Autoclavar
Guardar a 45º
Antes de usar añadir 1 mL de MgSO4 1 M
230
ANEXO 1
SM buffer
Preparación 1 mL gelatina 10 mg/ml
5 mL Tris-HCl 1 M
0,58 g NaCl
0,8 mL MgSO4 1 M
Enrasar a 100 mL con agua destilada
Autoclavar
231
ANEXO 1
12.4. Rastreo
12.4.1. Infección y Transferencia a Filtro
Para tener una probabilidad de 0,99 de obtener un fago positivo, se rastrean unos
75.000 fagos distribuidos en 10 placas (~7.500 calvas/placa). Los cálculos se muestran
en la Figura A5.
ln (1 – 0’99)
N= = 75.045 clons
Ln (1 – 13’5 / 2’2x105)
Figura A5. Cálculo del numero de fagos que hay que rastrear para obtenir como mínimo un
positivo (Sambrook et al. 1989, Klug y Cummings 1999).
232
ANEXO 1
1. Sembrar los fagos a hibridar en placas LB, a baja densidad para poder identificar
luego los fagos positivos. Crecer a 37ºC.
2. Guardar las placas 1 hora a 4ºC (para que el agar se endurezca).
3. Cortar filtros circulares del tamaño de las placas, uno por placa.
4. Numerar las placas y los filtros (estos últimos a lápiz), y hacer tres marcas
asimétricas en la placa que servirán de referencia para posicionar los positivos.
5. Abrir la placa y colocar el filtro correspondiente con el número hacia arriba y
haciéndolo coincidir con el de la placa. Dejar caer el filtro poco a poco para que
no queden burbujas, haciendo que toque primero el centro.
6. Dejar transferir 3 minutos.
7. Pinchar el filtro con una aguja al lado de cada una de las tres marcas hechas en la
placa, aprovechando el último pinchazo para levantar el filtro.
8. Coger el filtro con las pinzas y colocarlo boca arriba sobre papel Whatmann.
9. Dejar secar.
10. En dos cápsulas de petri de vidrio de mayor tamaño que los filtros poner dos
papeles Whatmann impregnados uno con solución de Desnaturalización (recién
preparada) y el otro con solución de Neutralización.
11. En una fiambrera grande de vidrio poner otro papel Whatmann con 2xSSC.
12. Poner los filtros, siempre boca arriba y controlando que la solución no los cubra
ni que se seque el papel: 5 minutos en sol. de desnaturalización
5 minutos en sol. de neutralización
15 minutos en 2xSSC
13. Dejar secar.
14. Fijar el ADN irradiando 3 minutos por cara con UV.
233
ANEXO 1
234
ANEXO 1
13. SECUENCIACIÓN
235
ANEXO 1
236
ANEXO 2. MAPAS FÍSICOS
ANEXO 2
Estos mapas se han realizado a partir de los clones BAC obtenidos de cada
región al rastrear la genoteca CHORI-225 de D. buzzatii. Los mapas se han utilizado
para escoger un BAC de cada región para secuenciar.
Actualmente se dispone de un mapa físico detallado de D. buzzatii realizado a
partir de esta misma genoteca. Los clones de la región labial-abdominalA se encuentran
en el contig 976 de dicho mapa y la región proboscipedia en el 968.
239
240
ANEXO 2
Figura A5. Mapa de la región lab-abdA. Arriba se muestra la posición de las sondas utilizadas y el gen al que corresponden. Las sondas de abdA se muestran
en relación a la estructura y tamaño del gen de D. melanogaster, única información disponible al realizar el mapa. Abajo se muestran los clones BAC, se
agrupan aquellos que contienen los mismos marcadores. Para la región del gen lab se han utilizado como sondas 3 clones en pBSK: pGPE307 (región 3’),
pGEP 302.1.1 (fragmento del intrón 1) y pGPE320.5 (extremo 5’). Para el gen abdA se ha utilizado una oligoprove diseñada sobre secuencia de D. virilis
(3’abdA) y un fragmento del gen (exones 2-4) amplificado a partir de DNA genómico de D. melanogaster (abdA2).
Figura A6. Mapa de la región pb. Arriba se muestra, como referencia, la estructura de parte del complejo ANT-C de D. melanogaster. Sobre los clones BAC
se muestra el mapa de restricción de NotI y AvrII de la región. Los puntos negros muestran la posición de los STS. La sonda utilizada para el reastreo de la
genoteca es un clon de cDNA de pb de D. buzzatii (pGPE 332).
241
ANEXO 2
ANEXO 2
242
ANEXO 3. SECUENCIAS
ANEXO 3
245
ANEXO 3