Chromosome evolution in Iberolacerta, a genus that deviates from the standard karyotype formula of Lacertidae.

This paper describes the preparation of flow-sorted chromosome paints from the Iberian Rock lizard Iberolacerta monticola, exemplifying their subsequent use in cross-species comparisons of chromosome painting. We carried out comparative analyses of chromosome evolution in the congeneric species I. galani and I. bonnali, as well as in two other species of Lacertini (Lacerta schreiberi and Timon lepidus) whose sex chromosomes were also studied through comparative genomic hybridization. Most species of Lacertini possess a diplod number of 2n = 38, with 36 acrocentric macrochromosomes and 2 microchromosomes. However, the nine species included in the genus Iberolacerta do not possess microchromosomes. Furthermore, very conspicuous differences from the standard Lacertini karyotype were observed in the three Pyrenean species of this genus, which included several biarmed metacentrics and a Z1Z2W multiple sex-chromosome system. With the possible exception of L. schreiberi, all the species of the family Lacertidae described to date appear to share homologous Z chromosomes, which date back to the last common ancestor of the whole group. We provide conclusive evidence that L. schreiberi should no longer be considered an exception to this rule, and demonstrate that the loss of microchromosomes in Iberolacerta was produced by their fusion to a middle-sized chromosome. Furthermore, we show that the multiple sex-chromosome system of the Pyrenean species of Iberolacerta originated from the fusion of the ancestral W chromosome with one of the shortest autosomes, and provide additional evidence of the fast evolution of DNA sequences linked to the W chromosome in Lacertini.

Phylogeography of the Spanish Moon Moth Graellsia isabellae (Lepidoptera, Saturniidae).

BACKGROUND: Geographic and demographic factors as well as specialisation to a new host-plant may lead to host-associated differentiation in plant-feeding insects. We explored the phylogeography of a protected moth, Graellsia isabellae, and its two recognised host-plant species (Pinus sylvestris and P. nigra) in order to seek for any concordance useful to disentangle the evolutionary history of this iconic lepidopteran. RESULTS: DNA variation in one mitochondrial marker and nine nuclear microsatellite loci revealed a strong phylogeographic pattern across 28 populations of G. isabellae studied in Spain and France comprising six groups mostly distributed along different mountain ranges. Reanalysis of a previously published chloroplast microsatellite dataset revealed a three and two-group structure for Spanish P. sylvestris and P. nigra, respectively. Overall, the population groupings of this protected moth did not match the ones of P. sylvestris and P. nigra. CONCLUSIONS: There was no evidence of host-associated differentiation between populations using P. sylvestris and the ones inhabiting P. nigra. The two major mitochondrial clades of G. isabellae likely diverged before the Last Glacial Maximum and geographically separated the species into a "southern" (Central and Southern Iberian clusters) and a "northern" lineage (Eastern Iberian, Pyrenean and French Alpine clusters). The Eastern Iberian System, where this insect uses both host-plants, harboured the highest level of genetic diversity. Such a group independently colonised the West and East parts of the Pyrenees. Our results point to a native origin for the French populations occurring in the Alps, genetically related to the Eastern Iberian and Pyrenean sites. The Central Iberian group derived from Southern Iberian ancestors. Secondary contacts were inferred between the Southern/Central Iberian populations and Eastern Iberian cluster as well as between the two Pyrenean ones. The mito-nuclear discordance observed with regard to the Eastern Iberian cluster is congruent with a secondary contact after the evolution of mito-nuclear incompatibilities in geographically isolated areas.

Evolutionary dynamics of two satellite DNA families in rock lizards of the genus Iberolacerta (Squamata, Lacertidae): different histories but common traits.

Satellite DNAs compose a large portion of all higher eukaryotic genomes. The turnover of these highly repetitive sequences is an important element in genome organization and evolution. However, information about the structure and dynamics of reptilian satellite DNA is still scarce. Two satellite DNA families, HindIII and TaqI, have been previously characterized in four species of the genus Iberolacerta. These families showed different chromosomal locations, abundances, and evolutionary rates. Here, we extend the study of both satellite DNAs (satDNAs) to the remaining Iberolacerta species, with the aim to investigate the patterns of variability and factors influencing the evolution of these repetitive sequences. Our results revealed disparate patterns but also common traits in the evolutionary histories of these satellite families: (i) each satellite DNA is made up of a library of monomer variants or subfamilies shared by related species; (ii) species-specific profiles of satellite repeats are shaped by expansions and/or contractions of different variants from the library; (iii) different turnover rates, even among closely related species, result in great differences in overall sequence homogeneity and in concerted or non-concerted evolution patterns, which may not reflect the phylogenetic relationships among taxa. Contrasting turnover rates are possibly related to genomic constraints such as karyotype architecture and the interspersed organization of diverging repeat variants in satellite arrays. Moreover, rapid changes in copy number, especially in the centromeric HindIII satDNA, may have been associated with chromosomal rearrangements and even contributed to speciation within Iberolacerta.

Isolation and characterization of two satellite DNAs in some Iberian rock lizards (Squamata, Lacertidae).

Satellite DNAs represent a large portion of all high eukaryotic genomes. They consist of numerous very similar repeated sequences, tandemly arranged in large clusters up to 100 million base pairs in length, usually located in the heterochromatic parts of chromosomes. The biological significance of satDNAs is still under discussion, but most of their proposed functions are related to heterochromatin and/or centromere formation and function. Because information about the structure of reptilian satDNA is far from exhaustive, we present a molecular and cytogenetic characterization of two satDNA families in four lacertid species. Two families of tandemly repeated DNAs, namely TaqI and HindIII satDNAs, have been cloned and sequenced from four species belonging to the genus Iberolacerta. These satDNAs are characterized by a monomer length of 171-188 and 170-172 bp, and by an AT content of 60.5% and 58.1%, respectively. FISH experiments with TaqI satDNA probe produced bright signals in pericentromeric regions of a subset of chromosomes whereas all the centromeres were marked by HindIII probe. The results obtained in this study suggest that chromosome location and abundance of satDNAs influence the evolution of these elements, with centromeric families evolving tenfold faster than interstitial/pericentromeric ones. Such different rates render different satellites useful for phylogenetic investigation at different taxonomic ranks.

Evidence for the persistence of an active endogenous retrovirus (ERVE) in humans.

Transposable elements (TEs) account for nearly half (44 %) of the human genome. However, their overall activity has been steadily declining over the past 35-50 million years, so that <0.05 % of TEs are presumably still "alive" (potentially transposable) in human populations. All the active elements are retrotransposons, either autonomous (LINE-1 and possibly the endogenous retrovirus ERVK), or non-autonomous (Alu and SVA, whose transposition is dependent on the LINE-1 enzymatic machinery). Here we show that a lineage of the endogenous retrovirus ERVE was recently engaged in ectopic recombination events and may have at least one potentially fully functional representative, initially reported as a novel retrovirus isolated from blood cells of a Chinese patient with chronic myeloid leukemia, which bears signals of positive selection on its envelope region. Altogether, there is strong evidence that ERVE should be included in the short list of potentially active TEs, and we give clues on how to identify human specific insertions of this element that are likely to be segregating in some of our populations.

A new mountain lizard from Montes de León (NW Iberian Peninsula): Iberolacerta monticola astur ssp. nov. (Squamata: Lacertidae).

Iberolacerta populations from the Northern Montes de Leon (NML) were studied by means of external morphology (scalation and biometry), osteology and genetics (mtDNA and microsatellites), searching for their homogeneity ("intrazonalanalysis") and, once verified, comparing them with Iberolacerta monticola s. str. (from Central Cantabrian Mountains)and/. gal ani (from Southern Montes de Leon) ("extrazonal analysis") from neighboring areas.Our "intrazonal analysis" revealed discordances between the different approaches, especially the patterns of variation of nuclear microsatellites (congruent with external morphology) and mtDNA, namely a very low nuclear differentiation between relatively highly differentiated mtDNA lineages. The morphological approach was unable to discriminate any of the populations as significantly different from the others in the NML. Mitochondrial DNA revealed a haplotype lineage closely related to I. galani (MNL-II in our text) in some specimens of Sierra de Villabandfn and Suspiron, but these populations are morphologically indistinguishable from the main part of the other populations that belong to lineage NML-1,phylogenetically closer to/. monticola. After a separation from I. manti cola ca. 1.8 Mya, the populations in this geographic region must have suffered at least two different waves of gene flow from I. gal ani, the second one not much later than 0.5 Mya. Microsatellite results indicate that all the NML populations are genetically similar in terms of their nuclear genomes,independently of their mitochondrial differentiation (NML-I vs. NML-II haplotype groups). Since all the morphological and microsatellite evidences point towards the fact that, independently of the mitochondrial haplotypes that they bear (NML-1 or NML-II), there is only one taxon in the area, we describe it as: Iberolacerta monticola astur ssp. nov.Concerning the relationships of I. m. astur ssp. nov. with I. monticola s. str. and I. gal ani ("extra zonal analysis"), in the female analyses the new taxon centroid is closer to I. monticola s. str. than to I. gal ani (more similarity with I manticolas.str.), whereas in the male analyses the relationship is just the contrary (closer to I. gal ani, paralleling the direction of the hypothesized past hybridization). Moreover, in both sexes' ANOVA, I. m. astur ssp. nov. results more similar (lessP<0.05 differences) to I. galani than to I. monticola s. str. Osteologically, I. m. astur ssp. nov. is slightly more similar toI. monticola s. str. than to I. galani, especially in the squamosal bone, which is regularly arched (primitive shape). Genetically,as indicated above, the NML populations can be subdivided in two groups according to their mitochondrial DNA,namely NML-I (bearing clearly differentiated haplotypes, phylogenetically closer to I. monticola) and NML-II (whose haplotypes could have been mistaken for those of an I. gal ani population). This mitochondrial subdivision has at most a subtle nuclear correlate, however. According to the nuclear microsatellite markers, all the NML populations belong to a single group(/. m. astur ssp. nov.), which would be more similar to I. gal ani than to I monticola, with NML-II populations lying closer to I. galani than those from the NML-I group and, correspondingly, more distant from I. monticola. The discordant phylogenetic signal of mitochondrial and nuclear markers is discussed in terms of past introgression events and sex-biases in phylopatry and dispersion in these species. Iberolacerta manti cola astur ssp. nov., inhabits the Northern Montes de Leon (Sierra de Gistreo sensu latissimo ): Gistredo,Catoute, Tambaron, Nevadfn, Villabandfn (or Macizo del Alto de Ia Canada), Arcos del Agua (or Fernan Perez),Tiendas and Suspiron, mainly in quartzite and slate rock substrates. Its current distribution, cornered in the NW of theNorthern part of the Montes de Leon, suggests a possible competitive exclusion between this taxon and/. galani, as the galani haplotypes (NML-II) appear cornered in the most harsh and continental areas, speaking also about a, even in the past, very limited presence of this species in the area that probably was soon absorbed by I. m. astur ssp. nov. (with NMLI haplotypes). Variation in watershed limits (especially with l montico/a s. str. in the North) and Pleistocene climatic oscillations(with I. gal ani in the South) probably played a crucial role in isolation of the different Iberolacerta colonizationwaves in this zone. These changes in the boundaries among watersheds limited the contact between the NML and the main Cantabrian Mountains, restricting to narrow points (different along time) the contact between the two ranges, and thus,the areas for possible contact between I. m. astur ssp. nov. and I. monticola s. str. (see Fig. lB). The origin of this tax on dates back to the end of Pliocene or Lower Pleistocene (around 1.8 Mya), according to mtDNA divergence. On the other side, climatic oscillations allowed expansion and contact with the more continental harsh climate-dwelling I. gal ani.

Evolution of the mdg1 lineage of the Ty3/gypsy group of LTR retrotransposons in Anopheles gambiae.

So far, only a few retrovirus-like transposable elements (TEs) have been reported in Anopheles mosquitoes, although a large fraction of their genomes is made up of these middle repetitive sequences. By screening the A. gambiae genome databases, we have found 10 element families belonging to the mdg1 lineage of the Ty3/gypsy group of long terminal repeat (LTR) retrotransposons. These Anopheles families constitute a sister clade of the Drosophila representatives of this same lineage. According to the phylogenetic reconstruction of their open reading frame (ORF)2 enzymatic domains, the analysis of patterns of nucleotide substitution therein, and the estimation of the age of particular insertions, all these elements must have been active until quite recently, and some of them must be very young. On the other hand, the fact that all these element families are primarily composed of fragmentary copies (mostly solos) or full-length copies with inactivating mutations indicates that their turnover rate has been probably very low. Finally, incongruent phylogenies obtained from different regions of the elements strongly suggest that recombination has played a significant role in their evolutionary history.

Causes and evolutionary consequences of population subdivision of an Iberian mountain lizard, Iberolacerta monticola.

AIM: The study of the factors that influence population connectivity and spatial distribution of genetic variation is crucial for understanding speciation and for predicting the effects of landscape modification and habitat fragmentation, which are considered severe threats to global biodiversity. This dual perspective is obtained from analyses of subalpine mountain species, whose present distribution may have been shaped both by cyclical climate changes over ice ages and anthropogenic perturbations of their habitats. Here, we examine the phylogeography, population structure and genetic diversity of the lacertid lizard Iberolacerta monticola, an endemism considered to be facing a high risk of extinction in several populations. LOCATION: Northwestern quadrant of the Iberian Peninsula. METHODS: We analyzed the mtDNA variation at the control region (454 bp) and the cytochrome b (598 bp) loci, as well as at 10 nuclear microsatellite loci from 17 populations throughout the distribution range of the species. RESULTS: According to nuclear markers, most sampling sites are defined as distinct, genetically differentiated populations, and many of them show traces of recent bottlenecks. Mitochondrial data identify a relatively old, geographically restricted lineage, and four to six younger geographically vicariant sister clades, whose origin may be traced back to the mid-Pleistocene revolution, with several subclades possibly associated to the mid-Bruhnes transition. Geographic range fragmentation of one of these clades, which includes lowland sites, is very recent, and most likely due to the accelerated loss of Atlantic forests by human intervention. MAIN CONCLUSIONS: Altogether, the data fit a "refugia within refugia" model, some lack of pattern uniformity notwithstanding, and suggest that these mountains might be the cradles of new species of Iberolacerta. However, the changes operated during the Holocene severely compromise the long-term survival of those genetic lineages more exposed to the anthropogenic perturbations of their habitats.

Evolutionary dynamics of the Ty3/gypsy LTR retrotransposons in the genome of Anopheles gambiae.

Ty3/gypsy elements represent one of the most abundant and diverse LTR-retrotransposon (LTRr) groups in the Anopheles gambiae genome, but their evolutionary dynamics have not been explored in detail. Here, we conduct an in silico analysis of the distribution and abundance of the full complement of 1045 copies in the updated AgamP3 assembly. Chromosomal distribution of Ty3/gypsy elements is inversely related to arm length, with densities being greatest on the X, and greater on the short versus long arms of both autosomes. Taking into account the different heterochromatic and euchromatic compartments of the genome, our data suggest that the relative abundance of Ty3/gypsy LTRrs along each chromosome arm is determined mainly by the different proportions of heterochromatin, particularly pericentric heterochromatin, relative to total arm length. Additionally, the breakpoint regions of chromosomal inversion 2La appears to be a haven for LTRrs. These elements are underrepresented more than 7-fold in euchromatin, where 33% of the Ty3/gypsy copies are associated with genes. The euchromatin on chromosome 3R shows a faster turnover rate of Ty3/gypsy elements, characterized by a deficit of proviral sequences and the lowest average sequence divergence of any autosomal region analyzed in this study. This probably reflects a principal role of purifying selection against insertion for the preservation of longer conserved syntenyc blocks with adaptive importance located in 3R. Although some Ty3/gypsy LTRrs show evidence of recent activity, an important fraction are inactive remnants of relatively ancient insertions apparently subject to genetic drift. Consistent with these computational predictions, an analysis of the occupancy rate of putatively older insertions in natural populations suggested that the degenerate copies have been fixed across the species range in this mosquito, and also are shared with the sibling species Anopheles arabiensis.

Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics.

Culex quinquefasciatus (the southern house mosquito) is an important mosquito vector of viruses such as West Nile virus and St. Louis encephalitis virus, as well as of nematodes that cause lymphatic filariasis. C. quinquefasciatus is one species within the Culex pipiens species complex and can be found throughout tropical and temperate climates of the world. The ability of C. quinquefasciatus to take blood meals from birds, livestock, and humans contributes to its ability to vector pathogens between species. Here, we describe the genomic sequence of C. quinquefasciatus: Its repertoire of 18,883 protein-coding genes is 22% larger than that of Aedes aegypti and 52% larger than that of Anopheles gambiae with multiple gene-family expansions, including olfactory and gustatory receptors, salivary gland genes, and genes associated with xenobiotic detoxification.

Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae.

The recent availability of the genome of Anopheles gambiae offers an extraordinary opportunity for comparative studies of the diversity of transposable elements (TEs) and their evolutionary dynamics between two related species, taking advantage of the existing information from Drosophila melanogaster. To this goal, we screened the genome of A. gambiae for elements belonging to the Ty3/gypsy group of long-terminal repeat (LTR) retrotransposons. The A. gambiae genome displays a rich diversity of LTR retrotransposons, clearly greater than D. melanogaster. We have characterized in detail 63 families, belonging to five of the nine main lineages of the Ty3/gypsy group. The Mag lineage is the most diverse and abundant, with more than 30 families. In sharp contrast with this finding, a single family belonging to this lineage has been found in D. melanogaster, here reported for the first time in the literature, most probably consisting of old inactive elements. The CsRn1 lineage is also abundant in A. gambiae but almost absent from D. melanogaster. Conversely, the Osvaldo lineage has been detected in Drosophila but not in Anopheles. Comparison of structural characteristics of different families led to the identification of several lineage-specific features such as the primer-binding site (PBS), the gag-pol translational recoding signal (TRS), which is extraordinarily diverse within the Ty3/gypsy retrotransposons of A. gambiae, or the presence/absence of specific amino acid motifs. Interestingly, some of these characteristics, although in general well conserved within lineages, may have evolved independently in particular branches of the phylogenetic tree. We also show evidence of recent activity for around 75% of the families. Nevertheless, almost all families contain a high proportion of degenerate members and solitary LTRs (solo LTRs), indicative of a lower turnover rate of retrotransposons belonging to the Ty3/gypsy group in A. gambiae than in D. melanogaster. Finally, we have detected significant overrepresentations of insertions on the X chromosome versus autosomes and of putatively active insertions on euchromatin versus heterochromatin.

Paleogenomic record of the extinction of human endogenous retrovirus ERV9.

An outstanding question of genome evolution is what stops the invasion of a host genome by transposable elements (TEs). The human genome, harboring the remnants of many extinct TE families, offers an extraordinary opportunity to investigate this problem. ERV9 is an endogenous retrovirus repeatedly mobilized during primate evolution, 15 to 6 million years ago (MYA), which left a trace of over a hundred provirus-like copies and at least 4,000 solitary long terminal repeats (LTRs) in the human genome. Then, its proliferation ceased for unknown reasons, and the family went extinct. We have made a detailed reconstruction of its last active subfamily, ERV9_XII, by examining 115 solitary LTRs from it. These insertions were grouped into 11 sets according to shared nucleotide variants, which could be placed in a sequential order of 10 to 6 MYA. At least 75% of the subfamily was produced 8 to 6 MYA, during a stage of intense proliferation. With new analytical tools, we show that the youngest and most prolific sets may have been produced by effectively instantaneous expansions of corresponding single-sequence variants. The extinction of this family apparently was not a consequence of its slow gradual degeneration, but the outcome of the fixation of specific restrictive alleles in the human-chimpanzee ancestral population. Three species-specific insertions (two in humans and one in chimpanzees) were identified, further supporting that extinction took place when these two species were beginning to diverge. These are the only fixed differences of this kind so far observed between humans and chimpanzees, apart from those belonging to the human endogenous retrovirus K family.

On the relative roles of faster-X evolution and dominance in the establishment of intrinsic postzygotic isolating barriers.

The modern theory of speciation assigns a prominent role to the recessivity of genetic incompatibilities in the two rules of speciation, namely Haldane's rule and the 'large X effect', and considers that the contribution of faster evolution of the X versus the autosomes to those patterns is generally of relatively minor importance. By extending Turelli and Orr's previous analysis of the model of two-locus Dobzhansky-Muller incompatibilities, I first show that when the X and the autosomes evolve at the same rate, the two dominance parameters involved in that model are not equally important for the declaration of a large X effect, but that the degree of recessivity of homozygous-homozygous incompatibilities is the major determinant for such a declaration. When the X evolves faster than the autosomes, the model obviously predicts that the importance of both dominance parameters will progressively vanish. It is then of importance to obtain estimates of the relative evolutionary rate of X-linked incompatibility loci. Several different procedures to obtain such estimates from the perspective of the large X effect are suggested. The application of the appropriate test to the only suitable data from Drosophila hybridizations so far available leads to the conclusion that the X actually evolves at least 2.5 times faster than the autosomes, as far as hybrid male sterility determinants are concerned, thus making dominance considerations absolutely irrelevant. Notwithstanding the necessity of further tests, the relative roles currently assigned to faster-X evolution and dominance in the theory of speciation should be revised, giving due prominence to faster-X evolution, at least for hybrid male sterility in the genus Drosophila.