Molecular insights into the colonization and chromosomal diversification of Madeiran house mice

The colonization history of Madeiran house mice was investigated by analysing the complete mitochondrial (mt) D-loop sequences of 156 mice from the island of Madeira and mainland Portugal, extending on previous studies. The numbers of mtDNA haplotypes from Madeira and mainland Portugal were substantially increased (17 and 14 new haplotypes respectively), and phylogenetic analysis confirmed the previously reported link between the Madeiran archipelago and northern Europe. Sequence analysis revealed the presence of four mtDNA lineages in mainland Portugal, of which one was particularly common and widespread (termed the 'Portugal Main Clade'). There was no support for population bottlenecks during the formation of the six Robertsonian chromosome races on the island of Madeira, and D-loop sequence variation was not found to be structured according to karyotype. The colonization time of the Madeiran archipelago by Mus musculus domesticus was approached using two molecular dating methods (mismatch distribution and Bayesian skyline plot). Time estimates based on D-loop sequence variation at mainland sites (including previously published data from France and Turkey) were evaluated in the context of the zooarchaeological record of M. m. domesticus. A range of values for mutation rate (μ) and number of mouse generations per year was considered in these analyses because of the uncertainty surrounding these two parameters. The colonization of Portugal and Madeira by house mice is discussed in the context of the best-supported parameter values. In keeping with recent studies, our results suggest that mutation rate estimates based on interspecific divergence lead to gross overestimates concerning the timing of recent within-species events.     

Is chromosomal speciation occurring in house mice in Tunisia?

Two chromosomal races of house mice are present in Tunisia, one represented by mice carrying the 40-acrocentric standard karyotype and the other by a Robertsonian race (2re = 22) homozygous for nine centric fusions (Rb). A comparative summary on allozyme divergence, geographical distribution and level of reproductive isolation in the Tunisian and European Rb races is presented, to which new data on mitochondrial DNA and morphological divergence are added. The Tunisian 22Rb race revealed unique features not matched by the European chromosomal races, such as a decrease in allozymic variability, a higher level of genetic and morphological differentiation and a mosaic geographical distribution. The mtDNA analysis argued in favour of a local origin of the chromosomal divergence suggesting that the higher level of differentiation between the Tunisian races resulted from the older age of the 22Rb race and/or from a severe botdeneck. The decrease in fertility of chromosomal hybrids between the Tunisian races was compatible with the limited genetic introgression between diem. Furthermore, data on the restricted distribution of hybrid populations suggested that premating reproductive barriers may be evolving. The Tunisian 22Rb race is thus an appropriate model to investigate a chromosomally-mediated speciation event.     

DEVELOPMENTAL INSTABILITY IN WILD CHROMOSOMAL HYBRIDS OF THE HOUSE MOUSE

In wild populations of the house mouse from Tunisia, fluctuating asymmetry and character size of tooth traits were compared between chromosomal races (2n = 40, all acrocentric standard karyotype, and 2n = 22, with nine fixed Robertsonian fusions) and their natural hybrids. Developmental stability was impaired in hybrids compared to both parental groups. Because genetic divergence measured by allozyme markers was low, genomic incompatibilities were not expected between the chromosomal races. This suggests that differentiation of gene systems specifically involved in development may have occurred between the chromosomal races. Support for the latter was found in the study of character size which showed that the 2n = 22 mice had smaller teeth than either the hybrid or the standard mice. The study of Tunisian chromosomal races thus shows that chromosomal evolution may lead to important changes in coadapted gene systems without involving extensive genic differentiation.     

Lack of congruence between morphometric evolution and genetic differentiation suggests a recent dispersal and local habitat adaptation of the Madeiran lizard Lacerta dugesii

Genetic differentiation among nine populations of the endemic lizard Lacerta dugesii Milne-Edwards 1829 (Lacertidae) from four groups of islands constituting the Archipelago of Madeira, was investigated by protein electrophoresis at 23 enzyme loci. Among twenty polymorphic loci, the total genetic diversity was due primarily to intra-population variation. The allele and genotypic frequencies among populations showed some heterogeneity, allowing the species to present a structuring pattern compatible with their geographical clustering. Some evidence suggests that selection acting on some loci in different ecological conditions may be responsible for the clustering of the populations studied. There was no apparent isolation effect expected under an "island" model of population divergence, and no correlation was found between genetic and geographic distances among populations. Morphological variation of the proposed three L. dugesii subspecies is not congruent with the allozyme analysis. This most probably suggests a rapid colonization of the islands followed by a strong effect of selection operating over the morphological characters used to define the subspecies.     

Mitochondrial portraits of the Madeira and Açores archipelagos witness different genetic pools of its settlers

We have studied the matrilineal genetic composition of the Madeira and Açores north Atlantic archipelagos, which were settled by the Portuguese in the 15th century. Both archipelagos, and particularly Madeira, were involved in a complex commercial network established by the Portuguese, which included the trading of slaves across the Atlantic. One hundred and fifty-five mtDNAs sampled from the Madeira and 179 from the Açores archipelagos were analysed for the hypervariable segment I (HVS-I), and for haplogroup-diagnostic coding-region RFLPs. The different settlement histories of both groups of islands are well reflected in their present day mtDNA pool. Although both archipelagos show identical diversity values, they are clearly different in their haplogroup content. Madeira displays a stronger sub-Saharan imprint, with haplogroups L1–L3 constituting about 13% of the lineages. Also, the relative frequencies of L sub-clusters in Madeira and mainland Portugal suggests that, at least in part, African presence in Madeira can be attributed to a direct gene flow from West Africa and not via Portugal. A comparison of the genetic composition of these two archipelagos with the Canary Islands, specially taking into account that their European source population was essentially from the Iberian Peninsula, testifies the stronger impact of the North African U6 cluster in the Canaries. This group is present in Madeira at a moderate frequency, but very reduced in the Açores. Nevertheless the recorded introduction of Canary native Guanches, who are characterized by the presence of particular sub-clade U6b1, has left no detectable imprints in the present day population of Madeira.     

Influence of physical environmental characteristics and anthropogenic factors on the position and structure of a contact zone between two chromosomal races of the house mouse on the island of Madeira (North Atlantic, Portugal)

Aim This study aimed to investigate if and how environmental characteristics (physical factors of the natural environment and the human impact on the landscape) influence the position and structure of a contact zone between two chromosomal races of the house mouse (Mus musculus domesticus Rutty 1712) from the island of Madeira. Location The western part of Madeira, a volcanic island in the North Atlantic. Methods Mice were sampled along a south/north‐western transect following the main road, in human‐modified outdoor habitats. Karyotypes of mice were determined using the yeast‐stimulated bone marrow cell method. Trapping sites were characterized in terms of their physical (altitude, temperature, precipitation and soil type) and habitat (human landscape use and occupancy) features. Demographic parameters of mouse populations, based on trapping‐with‐removal techniques, were also analysed (relative abundance, sex‐ratio, juvenile ratio and female fertility ratio), as well as body size (weight and length). Results Four chromosomal zones were identified on the basis of the frequency of two diagnostic rearrangements (Rb(6.7) in race E. Calheta and Rb(7.15) in race A. Cruz). E. Calheta was present in the two southern‐most zones, followed by the contact zone characterized by the presence of two inter‐racial hybrids and the co‐occurrence of mice belonging to the two races. The northern‐most part of the transect was occupied by A. Cruz. Environmental features differed leading us to split the transect into two parts. The southern part is characterized by lower altitude and precipitation, milder temperature, better soil quality supporting vegetable crops and vineyards, and more abundant and evenly distributed human habitats. This southern part is occupied by E. Calheta mice. The north‐western part presents characteristics opposite to those described above with cereals as the main cultivated crop, and it includes the contact zone as well as the zone inhabited by A. Cruz mice. The demographic parameters evaluated in this study did not differ significantly between chromosomal zones. Main conclusions This ecological survey highlights differences in climatic and edaphic features that have moulded the agricultural activities of humans, contributing to a differentiation of their spatial development, and hence the structure of potential habitats for mice. Results are interpreted within the source–sink framework of population dynamics, following which E. Calheta may function as a source and the areas where A. Cruz and the contact zone are located may function as a sink. Our study suggests that the position and chromosomal composition of the contact zone is influenced by the human component underlying broader environmental features. Similar characteristics were most likely present during the historical settlement of Madeira. They may have favoured the independent divergence of the two races and influenced the dynamics of the contact zone.     

Genetic structure of house mouse (Mus musculus Linnaeus 1758) populations in the Atlantic archipelago of the Azores: colonization and dispersal

We analyzed the genetic structure and relationships of house mouse (Mus musculus) populations in the remote Atlantic archipelago of the Azores using nuclear sequences and microsatellites. We typed Btk and Zfy2 to confirm that the subspecies Mus musculus domesticus was the predominant genome in the archipelago. Nineteen microsatellite loci (one per autosome) were typed in a total of 380 individuals from all nine Azorean islands, the neighbouring Madeiran archipelago (Madeira and Porto Santo islands), and mainland Portugal. Levels of heterozygosity were high on the islands, arguing against population bottlenecking. The Azorean house mouse populations were differentiated from the Portuguese and Madeiran populations and no evidence of recent migration between the three was obtained. Within the Azores, the Eastern, Western, and Central island groups tended to act as separate genetic units for house mice, with some exceptions. In particular, there was evidence of recent migration events among islands of the Central island group, whose populations were relatively undifferentiated. Santa Maria had genetically distinctive mice, which may relate to its colonization history. © 2013 The Linnean Society of London     

Evolution rampant: house mice on Madeira

House mice are extra-ordinary animals – extra -ordinary in the literal sense of that word. They are pests – but also a valued laboratory animal. They are generalized rodents – and successful in habitats from tundra to tropics and from sea-level to high altitudes. They have differentiated into a perplexity of taxa, yet differ little in their general morphology. They were long scorned by ecologists as recently arrived commensals, but are increasingly illuminating evolutionary processes as new techniques are applied to their study. Local forms, once valued only by taxonomists, are proving ever more interesting as their genetics are probed. In 1992, Mathias & Mira described the apparently unexciting characteristics of mice living on the two main islands of the Madeira group, 600 km west of continental Portugal. Then in 2000, Britton-Davidian et al. discovered that there were at least six chromosomal (Robertsonian) races on the main island. In the past decade, studies of molecular and mitochondrial genomes have shown an array of variables and posed questions about the origins and subsequent evolution of these island mice. In this issue of Molecular Ecology , Förster et al. report on the mtDNA haplotypes found on the island and in mainland Portugal, discuss the probable source of the island colonizers, and consider data which might give information about the timing of the colonizing event(s).     

Fertility assessment in hybrids between monobrachially homologous Rb races of the house mouse from the island of Madeira: implications for modes of chromosomal evolution

The speciation model of divergence by monobrachially homologous fusions (that is, with one arm in common) benefits from a wide conceptual acceptance, because heterozygotes between populations carrying such fusions suffer from high levels of meiotic dysfunction. The same meiotic configurations can also be generated by WART (whole-arm reciprocal translocation), rearrangements that are known to occur in mammals. Estimating the disadvantage of heterozygotes carrying monobrachially homologous fusions is required to evaluate the relevance of this mode of chromosomal evolution in diversification and speciation. House mice are an excellent study models because chromosomal races exist carrying monobrachially homologous fusions, and WARTs have been documented in this species. The fertility of heterozygote mice carrying the smallest number of monobrachially homologous fusions (that is, a chain of four chromosomes, C4) was investigated in laboratory-bred hybrids between two parapatric chromosomal races from the island of Madeira. Meiotic nondisjunction analyses and histological sections of testes showed that aneuploidy (16.7%) and germ cell death (50.9%) rates reached significantly higher mean values in hybrids than in homozygotes. In females, however, the histological analysis of ovarian follicle parameters revealed no significant differences between hybrid and homozygous individuals. Overall, the reproductive assays indicated that these C4-carrying hybrids were not sterile but showed an approximately 50% decrease in fertility compared to homozygous parental mice. Implications for modes of chromosomal evolution involving monobrachially homologous fusions are discussed.     

Phylogeography of the Marsh Fritillary Euphydryas aurinia (Lepidoptera: Nymphalidae) in the UK

We sampled extant (and extinct) populations of Euphydryas aurinia to examine the phylogeography of the species in the UK. We were interested in whether the genetic structuring of populations reflects anything other than a single recent post-glacial colonization event. Four hundred base pairs of the mitochondrial Cytochrome b gene were sequenced from individuals from populations throughout the UK, as well as populations in France and Portugal and seven polymorphic allozyme loci resolved. The mean number of allozyme alleles per locus was 4 and isolation by distance was shown not to be a factor in the geographic structuring of genetic diversity either with or without the inclusion of the French data (Mantel statistic Z = 0.015 and 0.112 respectively, P >0.5). Cytochrome b nucleotide diversity (average number of nucleotide differences per site between two sequences) was low overall (0.003, n = 63) but mean cytochrome b gene diversity over all populations was 0.77. The presence of the Portuguese and French haplotype in Scottish populations indicates that the Iberian peninsula was likely to be one glacial refugium for E. aurinia populations. The pattern of mitochondrial DNA found in the UK could be interpreted in one of two ways: (1) two separate colonization events or (2) a single slow colonization event. Allele frequency distributions followed a similar geographical pattern as mitochondrial DNA haplotypes. An AMOVA assigned just 2.68% of allozyme genetic variation to the grouping tested, providing more support for the single colonization event theory.     

Origin and evolution of the Robertsonian populations of the house mouse (Rodentia, Muridae) in Tunisia based on allozyme studies

The central coastal region of Tunisia harbours two chromosomal races of the house mouse (2 n  = 22, 22Rb; 2 n  = 40, 40Std), which are genetically differentiated and show a high level of reproductive isolation. This study presents an allozyme analysis of house mice from the inland city of Kairouan 75 km from the coastal region. Results showed that the 22Rb and 40Std mice from Kairouan shared the same high degree of variability, and were not genetically differentiated. This contrasts with the genetic data registered for the two chromosomal races from the coastal towns, in which a particularly low genic diversity was observed in the 22Rb populations. As the two races in Kairouan show the same low hybridization rates as the populations in the coastal region, these results argue in favour of a local differentiation of the Rb race in Tunisia, most likely originating in Kairouan. The data further suggest that the decrease in variability observed in the coastal 22Rb populations is related to a founder effect or a bottleneck following the introduction of mice from Kairouan into one of these localities, prior to their dispersion throughout the coastal region.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 515–521.     

Genic differentiation and origin of Robertsonian populations of the house mouse (Mus musculus domesticus Rutty)

This paper examines the relation between chromosomal and nuclear-gene divergence in 28 wild populations of the house mouse semi-species, Mus musculus domesticus, in Western Europe and North Africa. Besides describing the karyotypes of 15 of these populations and comparing them to those of 13 populations for which such information was already known, it reports the results of an electrophoretic survey of proteins encoded by 34 nuclear loci in all 28 populations. Karyotypic variation in this taxon involves only centric (or Robertsonian) fusions which often differ in arm combination and number between chromosomal races. The electrophoretic analysis showed that the amount of genic variation within Robertsonian (Rb) populations was similar to that for all-acrocentric populations, i.e. bearing the standard karyotype. Moreover, divergence between the two types of populations was extremely low. These results imply that centric fusions in mice have not modified either the level or the nature of genic variability. The genetic similarity between Rb and all-acrocentric populations is not attributed to the persistence of gene flow, since multiple fusions cause marked reproductive isolation. Rather, we attribute this extreme similarity to the very recent origin of chromosomal races in Europe. Furthermore, genic diversity measures suggest that geographically separated Rb populations have in situ and independent origins. Thus, Rb translocations are probably not unique events, but originated repeatedly. Two models are presented to explain how the rapid fixation of a series of chromosomal rearrangements can occur in a population without lowering variability in the nuclear genes. The first model assumes that chromosomal mutation rates are between 10(-3) and 10(-4) and that populations underwent a series of transient bottlenecks in which the effective population size did not fall below 35. In the second model, genic variability is restored following severe bottlenecks, through gene flow and recombination.     

Molecular differentiation at nuclear loci in French host races of the European corn borer (Ostrinia nubilalis)

French populations of the European corn borer consist of two sympatric and genetically differentiated host races. As such, they are well suited to study processes that could be involved in sympatric speciation, but the initial conditions of host-race divergence need to be elucidated. Gene genealogies can provide insight into the processes involved in speciation. We used DNA sequences of four nuclear genes to (1) document the genetic structure of the two French host races previously delineated with allozyme markers, (2) find genes directly or indirectly involved in reproductive isolation between host races, and (3) estimate the time since divergence of the two taxa and see whether this estimate is compatible with this divergence being the result of a host shift onto maize after its introduction into Europe approximately 500 years ago. Gene genealogies revealed extensive shared polymorphism, but confirmed the previously observed genetic differentiation between the two host races. Significant departures from the predictions of neutral molecular evolution models were detected at three loci but were apparently unrelated to reproductive isolation between host races. Estimates of time since divergence between French host races varied from approximately 75,000 to approximately 150,000 years, suggesting that the two taxa diverged recently but probably long before the introduction of maize into Europe.     

How to Use Genetic Data to Distinguish Between Natural and Human-Mediated Introduction of <Emphasis Type="Italic">Littorina littorea</Emphasis> to North America

The rapid range southward expansion of the periwinkle Littorina littorea from the Canadian maritimes has fueled a long-running debate over whether this species was introduced to North America by human activity. A reappraisal of the mitochondrial DNA sequence evidence finds considerable endemic allelic diversity in the American population. The degree of endemic genetic diversity is higher than expected from human-mediated colonization, but not so much to suggest that it survived the last glacial maximum in America. Coalescent estimates of population divergence agree that colonization of America preceded European contact. A reappraisal of the ITS nuclear sequence data finds extensive recombination. Taking this recombination into account strengthens the genetic case against human-mediated introduction. Finally, a reappraisal of conflicting allozyme studies from the 1970’s supports a claim of limited divergence between American and European populations. This is consistent with post-glacial colonization, but the allozyme data cannot distinguish between natural or human-mediated colonization. Taken as a whole, the DNA sequence data supports the many sub-fossil reports of an American L. littorea population in the Canadian maritimes that preceded even the first visits by the Vikings.     

Parapatric diversification after post-glacial range expansion in the gall fly Urophora cardui (Tephritidae)

Aim Primary and secondary genetic clines in post-glacial colonized regions have different implications for biogeographic distributions and the origin of species. Primary clines arise in situ after colonization as adaptive responses to environmental gradients, while secondary clines are caused by contact between vicariant lineages. Here we analyse primary versus secondary origin of a genetic cline in the tephritid fly Urophora cardui in Jutland, Denmark, in a post-glacial landscape. Location Western Palaearctic. Methods Phylogeographic and demographic analyses of U. cardui based on mitochondrial DNA (mtDNA) genealogies, hierarchical genetic variance tests based on allozymes and distribution analysis of a rare allele from the Jutland cline. Results There was no phylogeographic divergence between the Jutland population of U. cardui north of the cline and neighbouring western European regional populations, which all shared the common western European mtDNA haplotype H1. At nuclear loci, by contrast, the North Jutland population was diverged above the mean level of divergence among regional populations and had no loss of genetic variation. A rare allozyme allele that was frequent in the cline area (up to 45%) and was missing north of the cline also occurred at low frequency (0–14%) elsewhere in the sampling range. Shallow phylogeographic divergence was observed between Russian and western European populations and between English and continental populations. Main conclusions The genetic variation patterns support primary cline evolution and parapatric divergence in Jutland following a demographic expansion of a western European ancestral source population of U. cardui, and suggest cryptic refugia and/or selection in other European population assemblages. The patterns of intra-specific regional divergence are discussed with respect to the interpretation of cryptic refugia in Europe after the most recent ice age.     

Allozyme variation of the pygmy wood mouse Sylvaemus uralensis (Rodentia, Muridae) and estimation of the divergence of its chromosome forms

The genetic divergence between the eastern European, southern European, and Asian chromosome forms of the pygmy wood mouse Sylvaemus uralensis, whose karyotypes differ from one another in the amount of pericentromeric heterochromatin, has been reevaluated using allozyme analysis. In general, Asian S. uralensis living in eastern Kazakhstan, eastern Turkmenistan (the Kugitang Ridge), and Uzbekistan are more monomorphic than European populations of this species. However, the allozyme differences between all chromosome forms of the pygmy wood mouse is comparable with the interpopulation differences within each form and are an order of magnitude smaller than those between "good" species of the genus Sylvaemus. Thus, the chromosome forms of S. uralensis cannot be considered to be separate species. The concept of races as large population groups that have not diverged enough to regard them as species but differ from one another in some genetic characters is used to describe the differentiation of S. uralensis forms more adequately. The currently available evidence suggests the existence of two S. uralensis races, the Asian and the European ones, and two chromosome forms (eastern and western) of the European race. The possible historical factors that have determined the formation of the races of the pygmy wood mouse are considered. According to the most plausible hypothesis, the shift and fragmentation of the broad-leaved forest zone during the most recent glacial period (late Pleistocene) were the crucial factors of the formation of these races, because they resulted in a prolonged isolation of the European and Asian population groups of S. uralensis from each other.     

Allozyme variation and racial differentiation in Swedish Carex lepidocarpa s.l. (Cyperaceae)

Two morphological races have previously been recognized within the sedge, Carex lepidocarpa , in Sweden. These largely allopatric races are accorded specific status, as C. lepidocarpa s.s. and C. jemtlandica , in Scandinavian floras. A study of allozyme variation in populations from 80 Swedish sites supports the morphological evidence for racial differentiation within C. lepidocarpa. The two races differ from each other in terms of allele frequencies at polymorphic loci and also show different levels of within-population genetic diversity. Material that is morphologically referable to C. lepidocarpa s.s. is characterized by relatively high levels of allozyme variation, both within and between populations. Carex lepidocarpa s.s. is widespread in southern Sweden. In contrast, material that is morphologically assignable to C. jemtlandica shows low levels of within-population genetic diversity, and there is little differentiation between the geographically separated isolates of C. jemtlandica in northern Sweden and on the Baltic island of Gotland. The high degree of morphological similarity and moderate levels of genetic differentiation between the two races within C. lepidocarpa indicate that it is more appropriate to recognize the races as subspecies than as species. The low levels of genetic variation in C. jemtlandica , both within and between populations, suggest that C. jemtlandica may have arisen from C. lepidocarpa (or a near ancestor of C. lepidocarpa) as a result of population fragmentation and isolation in glacial refugia, or during the process of post-glacial colonization of Scandinavia. Lack of allozyme evidence for extensive hybridization between the two races of C. lepidocarpa , despite their ability to hybridize freely where their ranges overlap at present, supports the suggestion that the two races have had separate post-glacial histories.     

Colonization history of atlantic island common chaffinches (Fringilla coelebs) revealed by mitochondrial DNA

Common chaffinches (Fringilla coelebs) are thought to have colonized the Atlantic island archipelagoes (the Azores, Madeira, and the Canaries) from neighboring continental populations (Iberia and north Africa) within the last million years. However, colonization may have occurred separately from north Africa to the Canaries and from Iberia to the Azores (as would be predicted geographically) or in one wave from Iberia to the Azores and then to Madeira and the Canaries. These alternatives have different implications for the evolution of morphometric and plumage differentiation in island chaffinches. To determine the most likely colonization route, we estimated the phylogenetic relationships among island and continental subspecies of common chaffinch using sequences from four mtDNA genes (cytochrome b, ATPase 6, NADH 5, and the control region). The most strongly supported mtDNA phylogeny places the continental subspecies together as the sister group to a monophyletic clade containing the island subspecies. This is consistent with a single wave of colonization, and suggests that patterns of similarity among Atlantic island common chaffinches, such as blue pigmentation, short wings, and long tarsi, are due to common colonization history rather than to convergent evolution in a common island environment. However, spectral analysis of phylogenetic splits showed that although monophyly of island haplotypes is favored, there is also substantial support for their polyphyletic origin. We attribute the latter to the confounding effect of homoplasy at multistate sites and to the relatively rapid sequence of colonization events which provided insufficient time for the accumulation of strong phylogenetic signal. These problems are likely to be significant impediments in attempts to test hypotheses of phylogenetic histories of recently evolved populations and taxa.     

Genetic diversity of the Macaronesian leafy liverwort Porella canariensis inferred from RAPD markers

Plant colonization of the North Atlantic raises the intriguing question of the relationships between extant island species with their continental counterparts (European, African, and American), which may provide clues to past geographic distribution and colonization history. It has been suggested that during past glaciations, many plant species with typical Mediterranean distributions survived in the Atlantic islands that belong to what is today known as Macronesia. We used random amplified polymorphic DNA (RAPD) markers to study 12 populations of the liverwort Porella canariensis partly covering its present-day distribution (Azores, Madeira, Canary and Cape Verde Islands, and Iberian Peninsula). Unweighted pair-group (UPGMA) and principal component (PCO) analyses showed a similar geographical pattern that suggested a close relationship between Iberian populations and those from the Canaries and Cape Verde Islands. Populations from Madeira had more genetic variation than those from the Azores, a result from either a richer diversity of habitats in Madeira, which prompted more population diversification, successive colonization waves from different origins, or an older colonization of Madeira. The data show that continuous patches of liverworts are often comprised of more than one individual. Finally, RAPDs can be used to investigate intraspecific diversity within a comparatively large geographic area and, with utmost care, can be used to infer a historic context to explain the patterns observed.