[en] Drosophila subobscura belongs to the obscura species group, which is one of the known ten species groups of the subgenus Sophophora of the genus Drosophila. Originally endemic from the Palearctic region, it has recently colonized North and South America. The chromosomal inversion system of D. subobscura represents one of the most interesting models to investigate the evolution of this type of genome rearrangement, because i) it shows extremely high levels of polymorphism, exhibiting inversions of all kinds regarding length and chromosomal location; and ii) it has been identified to be involved in the species' adaptation to contemporary global climate warming. The lack of a reference genome for the species has, however, stood as a major obstacle to its study. To overcome this limitation, here we have tackled a de novo assembly of the genome of D. subobscura using PacBio long-read technology. Raw PacBio reads were assembled using the Canu assembler. A semi-automated pipeline for assessing the quality of the contigs and scaffolding the genome has been developed that combined both synteny information from previously assembled genomes of D. melanogaster and D. pseudoobscura, and published data from 560 genetic markers derived from in situ hybridization experiments and genetic linkage analyses. Canu assembled contigs were then scaffolded resulting in 186 scaffolds with a N50 of approximately 6Mb. Chromosomal assignment, ordering and orientation of the scaffolds resulted in six pseudochromosomes, one for each of the six D. subobscura chromosomes. A total of 13,939 protein-coding genes were predicted, and 13% of the species genome was found to consist of repetitive sequences. Finally, the amounts of genome rearrangement between D. subobscura, D. guanche, D. melanogaster and D. pseudoobscura was assessed. The breakpoints of the fixed inversions between D. subobscura and D. guanche were determined and characterized. Here we illustrate that genome structure evolution in D. subobscura is driven indirectly, through the inversions' recombination-suppression effects in maintaining sets of adaptive alleles together in the face of gene flow. Presently, D. subobscura is experiencing a rapid replacement of high-latitude by low-latitude inversions associated with global warming. However, not all low-latitude inversions are correlated with the secular warming trend. The mixed behavior of O7 inversion across components of the ambient temperature suggests that it is driven by selective factors other than temperature alone. Research into this question has been hindered by lacking knowledge of the genomic breakpoint sequences of the inversion. To tackle this limitation, we generated a PacBio long read-based chromosome-scale genome assembly, from an O3+4+7 isogenic line. The complete continuous sequence of O7 was isolated using synteny analysis with the reference genome. Inversion O7 was shown to stretch 9.936 Mb, containing over 1,000 annotated genes. We illustrate that inversion O7 had a complex origin, involving multiple breaks associated with non-B DNA motifs, formation of a microinversion, and ectopic repair in trans with the two homologous chromosomes. Our findings support that inversion O7 breakpoints carry a pre-inversion record of fragility, including a sequence insertion, and transposition with later inverted duplication of an Attacin immunity gene. The O7 inversion was found to have relocated the major insulin signaling forkhead box subgroup O (foxo) gene bringing it in tight linkage with its antagonistic regulatory partner serine/threonine-protein kinase B (Akt1). Further, its distal breakpoint disrupted concerted evolution of the two inverted Attacin duplicates, reattaching them to dFOXO metabolic enhancers. We suggest that O7 exerts antagonistic pleiotropic effects on reproduction and immunity, setting a framework to understand its relationship with climate change. Our findings have general implications for current theories on the molecular mechanisms of formation of inversions and the contribution of breakage versus repair in shaping inversion-breakpoint junctions.

[es] Drosophila subobscura es una especie del grupo obscura, uno de los diez grupos de especies conocidos del subgénero Sophophora, género Drosophila. Originariamente endémica de la región paleártica, la especie colonizó recientemente América del Norte y del Sur. El sistema de inversiones cromosómicas de D. subobscura representa uno de los modelos más interesantes para investigar la evolución de este tipo de rearreglos estructurales del genoma, ya que i) muestra niveles inusualmente elevados de polimorfismo y ii) se ha detectado que está implicado en la adaptación de la especie al calentamiento global. Sin embargo la falta de un genoma de referencia para la especie ha supuesto un obstáculo importante para su estudio. Para superar esta limitación hemos abordado un ensamblado de novo del genoma de D. subobscura utilizando la tecnología de "long-reads" PacBio. Los "long-reads" se ensamblaron utilizando el ensamblador Canu. Para evaluar la calidad de los "contigs" y el "scaffolding" del genoma, se desarrolló un protocolo semiautomático. Los "contigs" generados por Canu obteniéndose 186 "scaffolds" con un N50 de aproximadamente 6Mb. La asignación cromosómica, el orden y la orientación de los "scaffolds" dieron como resultado seis pseudocromosomas, uno por cada uno de los seis cromosomas de D. subobscura. Se predijeron un total de 13.939 genes codificadores de proteínas, y se estimó que el 13% del genoma de la especie consiste en secuencias repetitivas. Finalmente, se determinaron y caracterizaron los puntos de rotura de las inversiones fijadas entre D. subobscura y D. guanche. Los resultados obtenidos indican que la evolución de la estructura del genoma en D. subobscura es una consecuencia indirecta de los efectos de supresión de la recombinación de las inversiones en el mantenimiento de conjuntos de alelos adaptativos frente al flujo génico. D. subobscura está experimentando actualmente un rápido reemplazamiento de inversiones de latitudes altas por inversiones de latitudes bajas asociado con el calentamiento global. No todas las inversiones de latitudes bajas siguen esta pauta, lo que sugiere que sus frecuencias no dependen solo de la temperatura. Uno de los ejemplos mejor documentados es la inversión O7, cuya investigación se ha visto obstaculizada por el desconocimiento de las secuencias de sus puntos de rotura genómica. Para abordar esta limitación, generamos un ensamblado del genoma a escala cromosómica a partir de una línea isogénica O3+4+7 utilizando tecnología de PacBio. Para aislar la secuencia completa de O7 se utilizó un análisis de conservación de sintenia con el genoma de referencia. La inversión O7 se originó mediante un mecanismo complejo que implicó múltiples roturas asociadas con motivos de ADN no-B, formación de una microinversión y reparación ectópica en trans con los dos cromosomas homólogos. Nuestros hallazgos respaldan que los puntos de rotura de la inversión O7 tenían una historia de fragilidad anterior al origen de la inversión, incluyendo una inserción seguida de una transposición con una duplicación invertida de un gen de inmunidad de la familia de las atacinas. Con respecto a alteraciones génicas, la inversión O7 reubicó el principal gen del subgrupo O (foxo) de señalización de la insulina, acercándolo a su regulador antagónico, la serina/treonina-proteína quinasa B (Akt1). Además, su punto de rotura distal interrumpió la evolución concertada de los dos duplicados invertidos de la atacina, conectándolos a elementos potenciadores de dFOXO. Estos descubrimientos sugieren que O7 ejerce efectos pleiotrópicos antagónicos sobre la reproducción y la inmunidad, lo cual aporta un marco para comprender su relación con el cambio climático. Nuestros hallazgos tienen implicaciones generales para las teorías actuales sobre los mecanismos moleculares de formación de las inversiones y el papel relativo de la rotura y la reparación en la formación de sus puntos de rotura.