Global genetic diversity of Aedes aegypti

Mosquitoes, especially Aedes aegypti, are becoming important models for studying invasion biology. We characterized genetic variation at 12 microsatellite loci in 79 populations of Ae. aegypti from 30 countries in six continents, and used them to infer historical and modern patterns of invasion. Our results support the two subspecies Ae. aegypti formosus and Ae. aegypti aegypti as genetically distinct units. Ae. aegypti aegypti populations outside Africa are derived from ancestral African populations and are monophyletic. The two subspecies co‐occur in both East Africa (Kenya) and West Africa (Senegal). In rural/forest settings (Rabai District of Kenya), the two subspecies remain genetically distinct, whereas in urban settings, they introgress freely. Populations outside Africa are highly genetically structured likely due to a combination of recent founder effects, discrete discontinuous habitats and low migration rates. Ancestral populations in sub‐Saharan Africa are less genetically structured, as are the populations in Asia. Introduction of Ae. aegypti to the New World coinciding with trans‐Atlantic shipping in the 16th to 18th centuries was followed by its introduction to Asia in the late 19th century from the New World or from now extinct populations in the Mediterranean Basin. Aedes mascarensis is a genetically distinct sister species to Ae. aegypti s.l. This study provides a reference database of genetic diversity that can be used to determine the likely origin of new introductions that occur regularly for this invasive species. The genetic uniqueness of many populations and regions has important implications for attempts to control Ae. aegypti, especially for the methods using genetic modification of populations.     

Historical environmental change in Africa drives divergence and admixture of Aedes aegypti mosquitoes: a precursor to successful worldwide colonization?

Increasing globalization has promoted the spread of exotic species, including disease vectors. Understanding the evolutionary processes involved in such colonizations is both of intrinsic biological interest and important to predict and mitigate future disease risks. The Aedes aegypti mosquito is a major vector of dengue, chikungunya and Zika, the worldwide spread of which has been facilitated by Ae. aegypti's adaption to human‐modified environments. Understanding the evolutionary processes involved in this invasion requires characterization of the genetic make‐up of the source population(s). The application of approximate Bayesian computation (ABC) to sequence data from four nuclear and one mitochondrial marker revealed that African populations of Ae. aegypti best fit a demographic model of lineage diversification, historical admixture and recent population structuring. As ancestral Ae. aegypti were dependent on forests, this population history is consistent with the effects of forest fragmentation and expansion driven by Pleistocene climatic change. Alternatively, or additionally, historical human movement across the continent may have facilitated their recent spread and mixing. ABC analysis and haplotype networks support earlier inferences of a single out‐of‐Africa colonization event, while a cline of decreasing genetic diversity indicates that Ae. aegypti moved first from Africa to the Americas and then to Asia. ABC analysis was unable to verify this colonization route, possibly because the genetic signal of admixture obscures the true colonization pathway. By increasing genetic diversity and forming novel allelic combinations, divergence and historical admixture within Africa could have provided the adaptive potential needed for the successful worldwide spread of Ae. aegypti.     

HUMAN IMPACTS HAVE SHAPED HISTORICAL AND RECENT EVOLUTION IN AEDES AEGYPTI,THE DENGUE AND YELLOW FEVER MOSQUITO

Although anthropogenic impacts are often considered harmful to species, human modifications to the landscape can actually create novel niches to which other species can adapt. These “domestication” processes are especially important in the context of arthropod disease vectors, where ecological overlap of vector and human populations may lead to epidemics. Here, we present results of a global genetic study of one such species, the dengue and yellow fever mosquito, Aedes aegypti, whose evolutionary history and current distribution have been profoundly shaped by humans. We used DNA sequences of four nuclear genes and 1504 single nucleotide polymorphism (SNP) markers developed with restriction‐site associated DNA (RAD) sequencing to test the hypothesis that Ae. aegypti originated in Africa, where a domestic form arose and spread throughout the tropical and subtropical world with human trade and movement. Results confirmed African ancestry of the species, and supported a single subspeciation event leading to the pantropical domestic form. In addition, genetic data strongly supported the hypothesis that human trade routes first moved domestic Ae. aegypti out of Africa into the New World, followed by a later invasion from the New World into Southeast Asia and the Pacific. These patterns of domestication and invasion are relevant to many species worldwide, as anthropogenic forces increasingly impact evolutionary processes.     

Rapid evolution and the genomic consequences of selection against interspecific mating

While few species introduced into a new environment become invasive, those that do provide critical information on ecological mechanisms that determine invasions success and the evolutionary responses that follow invasion. Aedes albopictus (the Asian tiger mosquito) was introduced into the naturalized range of Aedes aegypti (the yellow fever mosquito) in the United States in the mid‐1980s, resulting in the displacement of A. aegypti in much of the south‐eastern United States. The rapid displacement was likely due to the superior competitive ability of A. albopictus as larvae and asymmetric mating interference competition, in which male A. albopictus mate with and sterilize A. aegypti females, a process called “satyrization.” The goal of this study was to examine the genomic responses of a resident species to an invasive species in which the mechanism of character displacement is understood. We used double‐digest restriction enzyme DNA sequencing (ddRADseq) to analyse outlier loci between selected and control lines of laboratory‐reared A. aegypti females from two populations (Tucson, AZ and Key West, Florida, USA), and individual females classified as either “resisted” or “mated with” A. albopictus males via mating trials of wild‐derived females from four populations in Florida. We found significant outlier loci in comparing selected and control lines and between mated and nonmated A. aegypti females in the laboratory and wild‐derived populations, respectively. We found overlap in specific outlier loci between different source populations that support consistent genomic signatures of selection within A. aegypti. Our results point to regions of the A. aegypti genome and potential candidate genes that may be involved in mating behaviour, and specifically in avoiding interspecific mating choices.     

Ancestral origins and invasion pathways in a globally invasive bird correlate with climate and influences from bird trade

Invasive species present a major threat to global biodiversity. Understanding genetic patterns and evolutionary processes that reinforce successful establishment is paramount for elucidating mechanisms underlying biological invasions. Among birds, the ring‐necked parakeet (Psittacula krameri) is one of the most successful invasive species, established in over 35 countries. However, little is known about the evolutionary genetic origins of this species and what population genetic signatures tell us about patterns of invasion. We reveal the ancestral origins of populations across the invasive range and explore the potential influence of climate and propagule pressure from the pet trade on observed genetic patterns. Ring‐necked parakeet samples representing the ancestral native range (n = 96) were collected from museum specimens, and modern samples from the invasive range (n = 855) were gathered from across Europe, Mauritius and Seychelles, and sequenced for two mitochondrial DNA markers comprising 868 bp of cytochrome b and control region, and genotyped at 10 microsatellite loci. Invasive populations comprise birds that originate predominantly from Pakistan and northern areas of India. Haplotypes associated with more northerly distribution limits in the ancestral native range were more prevalent in invasive populations in Europe, and the predominance of Asian haplotypes in Europe is consistent with the higher number of Asian birds transported by the pet trade outside the native range. Successful establishment of invasive species is likely to be underpinned by a combination of environmental and anthropogenic influences.     

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo‐Pacific populations of Aedes aegypti (L.)

Nations throughout the Indo‐Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target‐site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo‐Pacific region are described and their geographical distributions investigated using genome‐wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed‐models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome‐wide variation. Mixed‐models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome‐wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo‐Pacific region, pointing to undocumented mosquito invasions between countries.     

Wolbachia pipientis occurs in Aedes aegypti populations in New Mexico and Florida,USA

The mosquitoes Aedes aegypti (L.) and Ae. albopictus Skuse are the major vectors of dengue, Zika, yellow fever, and chikungunya viruses worldwide. Wolbachia, an endosymbiotic bacterium present in many insects, is being utilized in novel vector control strategies to manipulate mosquito life history and vector competence to curb virus transmission. Earlier studies have found that Wolbachia is commonly detected in Ae. albopictus but rarely detected in Ae. aegypti. In this study, we used a two‐step PCR assay to detect Wolbachia in wild‐collected samples of Ae. aegypti. The PCR products were sequenced to validate amplicons and identify Wolbachia strains. A loop‐mediated isothermal amplification (LAMP) assay was developed and used for detecting Wolbachia in selected mosquito specimens as well. We found Wolbachia in 85/148 (57.4%) wild Ae. aegypti specimens from various cities in New Mexico, and in 2/46 (4.3%) from St. Augustine, Florida. Wolbachia was not detected in 94 samples of Ae. aegypti from Deer Park, Harris County, Texas. Wolbachia detected in Ae. aegypti from both New Mexico and Florida was the wAlbB strain of Wolbachia pipientis. A Wolbachia‐positive colony of Ae. aegypti was established from pupae collected in Las Cruces, New Mexico, in 2018. The infected females of this strain transmitted Wolbachia to their progeny when crossed with males of Rockefeller strain of Ae. aegypti, which does not carry Wolbachia. In contrast, none of the progeny of Las Cruces males mated to Rockefeller females were infected with Wolbachia.     

Population genetics of Anopheles koliensis through Papua New Guinea: New cryptic species and landscape topography effects on genetic connectivity

New Guinea is a topographically and biogeographically complex region that supports unique endemic fauna. Studies describing the population connectivity of species through this region are scarce. We present a population and landscape genetic study on the endemic malaria‐transmitting mosquito, Anopheles koliensis (Owen). Using mitochondrial and nuclear sequence data, as well as microsatellites, we show the evidence of geographically discrete population structure within Papua New Guinea (PNG). We also confirm the existence of three rDNA ITS2 genotypes within this mosquito and assess reproductive isolation between individuals carrying different genotypes. Microsatellites reveal the clearest population structure and show four clear population units. Microsatellite markers also reveal probable reproductive isolation between sympatric populations in northern PNG with different ITS2 genotypes, suggesting that these populations may represent distinct cryptic species. Excluding individuals belonging to the newly identified putative cryptic species (ITS2 genotype 3), we modeled the genetic differences between A. koliensis populations through PNG as a function of terrain and find that dispersal is most likely along routes with low topographic relief. Overall, these results show that A. koliensis is made up of geographically and genetically discrete populations in Papua New Guinea with landscape topography being important in restricting dispersal.     

Sunshine versus gold: The effect of population age on genetic structure of an invasive mosquito

The genetic diversity and structure of invasive species are affected by the time since invasion, but it is not well understood how. We compare likely the oldest populations of Aedes aegypti in continental North America with some of the newest to illuminate the range of genetic diversity and structure that can be found within the invasive range of this important disease vector. Aedes aegypti populations in Florida have probably persisted since the 1600‐1700s, while populations in southern California derive from new invasions that occurred in the last 10 years. For this comparison, we genotyped 1,193 individuals from 28 sites at 12 highly variable microsatellites and a subset of these individuals at 23,961 single nucleotide polymorphisms (SNPs). This is the largest sample analyzed for genetic structure for either region, and it doubles the number of southern California populations previously analyzed. As predicted, the older populations (Florida) showed fewer indicators of recent founder effect and bottlenecks; in particular, these populations have dramatically higher genetic diversity and lower genetic structure. Geographic distance and driving distance were not good predictors of genetic distance in either region, especially southern California. Additionally, southern California had higher levels of genetic differentiation than any comparably sized documented region throughout the worldwide distribution of the species. Although population age and demographic history are likely driving these differences, differences in climate and transportation practices could also play a role.     

Genetic evidence for the origin of Aedes aegypti,the yellow fever mosquito,in the southwestern Indian Ocean

Aedes aegypti is among the best‐studied mosquitoes due to its critical role as a vector of human pathogens and ease of laboratory rearing. Until now, this species was thought to have originated in continental Africa, and subsequently colonized much of the world following the establishment of global trade routes. However, populations of this mosquito on the islands in the southwestern Indian Ocean (SWIO), where the species occurs with its nearest relatives referred to as the Aegypti Group, have received little study. We re‐evaluated the evolutionary history of Ae. aegypti and these relatives, using three data sets: nucleotide sequence data, 18,489 SNPs and 12 microsatellites. We found that: (a) the Aegypti Group diverged 16 MYA (95% HPD: 7–28 MYA) from its nearest African/Asian ancestor; (b) SWIO populations of Ae. aegypti are basal to continental African populations; (c) after diverging 7 MYA (95% HPD: 4–15 MYA) from its nearest formally described relative (Ae. mascarensis), Ae. aegypti moved to continental Africa less than 85,000 years ago, where it recently (<1,000 years ago) split into two recognized subspecies Ae. aegypti formosus and a human commensal, Ae. aegypti aegypti; (d) the Madagascar samples form a clade more distant from all other Ae. aegypti than the named species Ae. mascarensis, implying that Madagascar may harbour a new cryptic species; and (e) there is evidence of introgression between Ae. mascarensis and Ae. aegypti on Réunion, and between the two subspecies elsewhere in the SWIO, a likely consequence of recent introductions of domestic Ae. aegypti aegypti from Asia.     

At least two indigenous species of the Bemisia tabaci complex are present in Brazil

Bemisia tabaci is one of the most important global agricultural insect pests, being a vector of emerging plant viruses such as begomoviruses and criniviruses that cause serious problems in many countries. Although knowledge of the genetic diversity of B. tabaci populations is important for controlling this pest and understanding viral epidemics, limited information is available on this pest in Brazil. A survey was conducted in different locations of São Paulo and Mato Grosso states, and the phylogenetic relationships of B. tabaci individuals from 43 populations sampled from different hosts were analysed based on partial mitochondrial cytochrome oxidase 1 gene (mtCOI) sequences. According to the recently proposed classification of the B. tabaci complex, which employs the 3.5% mtCOI sequence divergence threshold for species demarcation, most of the specimens collected were found to belong to the Middle East‐Asia Minor 1 species, which includes the invasive populations of the commonly known B biotype, within the Africa/Middle East/Asia Minor high‐level group. Three specimens collected from Solanun gilo and Ipomoea sp. were grouped together and could be classified in the New World species that includes the commonly known A biotype. However, six specimens collected from Euphorbia heterophylla, Xanthium cavanillesii and Glycine maxima could not be classified into any of the 28 previously proposed species, although according to the 11% mtCOI sequence divergence threshold, they belong to the New World high‐level group. These specimens were classified into a new recently proposed species named New World 2 that includes populations from Argentina. Middle East‐Asia Minor 1, New World and New World 2 were differentiated by RFLP analysis of the mtCOI gene using TaqI enzyme. Taq I analysis in silico also differentiates these from Mediterranean species, thus making this method a convenient tool to determine population dynamics, especially critical for monitoring the presence of this exotic pest in Brazil.     

A layover in Europe: Reconstructing the invasion route of asexual lineages of a New Zealand snail to North America

Non‐native invasive species are threatening ecosystems and biodiversity worldwide. High genetic variation is thought to be a critical factor for invasion success. Accordingly, the global invasion of a few clonal lineages of the gastropod Potamopyrgus antipodarum is thus both puzzling and has the potential to help illuminate why some invasions succeed while others fail. Here, we used SNP markers and a geographically broad sampling scheme (N = 1617) including native New Zealand populations and invasive North American and European populations to provide the first widescale population genetic assessment of the relationships between and among native and invasive P. antipodarum. We used a combination of traditional and Bayesian molecular analyses to demonstrate that New Zealand populations harbour very high diversity relative to the invasive populations and are the source of the two main European genetic lineages. One of these two European lineages was in turn the source of at least one of the two main North American genetic clusters of invasive P. antipodarum, located in Lake Ontario. The other widespread North American group had a more complex origin that included the other European lineage and two New Zealand clusters. Altogether, our analyses suggest that just a small handful of clonal lineages of P. antipodarum were responsible for invasion across continents. Our findings provide critical information for prevention of additional invasions and control of existing invasive populations and are of broader relevance towards understanding the establishment and evolution of asexual populations and the forces driving biological invasion.     

Wolbachia infection alters the relative abundance of resident bacteria in adult Aedes aegypti mosquitoes,but not larvae

Insect–symbiont interactions are known to play key roles in host functions and fitness. The common insect endosymbiont Wolbachia can reduce the ability of several human pathogens, including arboviruses and the malaria parasite, to replicate in insect hosts. Wolbachia does not naturally infect Aedes aegypti, the primary vector of dengue virus, but transinfected Ae. aegypti have antidengue virus properties and are currently being trialled as a dengue biocontrol strategy. Here, we assess the impact of Wolbachia infection of Ae. aegypti on the microbiome of wild mosquito populations (adults and larvae) collected from release sites in Cairns, Australia, by profiling the 16S rRNA gene using next‐generation sequencing. Our data indicate that Wolbachia reduces the relative abundance of a large proportion of bacterial taxa in Ae. aegypti adults, that is in accordance with the known pathogen‐blocking effects of Wolbachia on a variety of bacteria and viruses. In adults, several of the most abundant bacterial genera were found to undergo significant shifts in relative abundance. However, the genera showing the greatest changes in relative abundance in Wolbachia‐infected adults represented a low proportion of the total microbiome. In addition, there was little effect of Wolbachia infection on the relative abundance of bacterial taxa in larvae, or on species diversity (accounting for species richness and evenness together) detected in adults or larvae. These results offer insight into the effects of Wolbachia on the Ae. aegypti microbiome in a native setting, an important consideration for field releases of Wolbachia into the population.     

New insights into the evolutionary history of Plasmodium falciparum from mitochondrial genome sequence analyses of Indian isolates

Estimating genetic diversity and inferring the evolutionary history of Plasmodium falciparum could be helpful in understanding origin and spread of virulent and drug‐resistant forms of the malaria pathogen and therefore contribute to malaria control programme. Genetic diversity of the whole mitochondrial (mt) genome of P. falciparum sampled across the major distribution ranges had been reported, but no Indian P. falciparum isolate had been analysed so far, even though India is highly endemic to P. falciparum malaria. We have sequenced the whole mt genome of 44 Indian field isolates and utilized published data set of 96 genome sequences to present global genetic diversity and to revisit the evolutionary history of P. falciparum. Indian P. falciparum presents high genetic diversity with several characteristics of ancestral populations and shares many of the genetic features with African and to some extent Papua New Guinean (PNG) isolates. Similar to African isolates, Indian P. falciparum populations have maintained high effective population size and undergone rapid expansion in the past with oldest time to the most recent common ancestor (TMRCA). Interestingly, one of the four single nucleotide polymorphisms (SNPs) that differentiates P. falciparum from P. falciparum‐like isolates (infecting non‐human primates in Africa) was found to be segregating in five Indian P. falciparum isolates. This SNP was in tight linkage with other two novel SNPs that were found exclusively in these five Indian isolates. The results on the mt genome sequence analyses of Indian isolates on the whole add to the current understanding on the evolutionary history of P. falciparum.     

Population genomics of the Asian tiger mosquito,Aedes albopictus: insights into the recent worldwide invasion

Aedes albopictus, the “Asian tiger mosquito,” is an aggressive biting mosquito native to Asia that has colonized all continents except Antarctica during the last ~30–40 years. The species is of great public health concern as it can transmit at least 26 arboviruses, including dengue, chikungunya, and Zika viruses. In this study, using double‐digest Restriction site‐Associated DNA (ddRAD) sequencing, we developed a panel of ~58,000 single nucleotide polymorphisms (SNPs) based on 20 worldwide Ae. albopictus populations representing both the invasive and the native range. We used this genomic‐based approach to study the genetic structure and the differentiation of Ae. albopictus populations and to understand origin(s) and dynamics of the recent invasions. Our analyses indicated the existence of two major genetically differentiated population clusters, each one including both native and invasive populations. The detection of additional genetic structure within each major cluster supports that these SNPs can detect differentiation at a global and local scale, while the similar levels of genomic diversity between native and invasive range populations support the scenario of multiple invasions or colonization by a large number of propagules. Finally, our results revealed the possible source(s) of the recent invasion in Americas, Europe, and Africa, a finding with important implications for vector‐control strategies.     

Isolation and characterization of microsatellite loci in Bemisia tabaci

Bemisia tabaci (Hemiptera: Aleyrodidae) is a haplo‐diploid species with a global distribution demonstrating strong geographical structure with eight recognizable genetic groups. Fifteen microsatellite loci (335 alleles, 6–44 alleles per locus) were derived from four of the eight groups and were then screened across 33 populations. These loci clearly differentiate the populations. The microsatellites amplified best in individuals from genetic groups rep‐resenting the Mediterranean, Middle East, Asia (three groups) and Australasia/Oceania and amplified less well with populations from sub‐Saharan Africa and the New World. This differential amplification pattern is a direct result of the relatedness to the microsatellite source material.     

Electrophysiological and behavioural response of Aedes albopictus to n‐heinecosane,an ovipositional pheromone of Aedes aegypti

Aedes aegypti (L.) and Aedes albopictus (Skuse) (Diptera: Culicidae) are highly anthropophilic mosquito species and potential vectors of dengue and yellow fever. The location of suitable sites for oviposition requires a set of visual, tactile, and olfactory cues that influence females before they lay their eggs. In this study, the effect of n‐heneicosane, a recognized oviposition pheromone of Ae. aegypti, on the olfactory receptors of the antennae of Ae. aegypti and Ae. albopictus was studied using electroantennographic detection coupled to gas chromatography (GC‐EAD). A significant electroantennographic response to n‐heneicosane in adult females of both mosquito species was observed. In addition, gravid Ae. albopictus females laid more eggs in substrate treated with n‐heneicosane at 0.1, 1, or 10 p.p.m. than in the control, denoting oviposition attractancy. Conversely, at 30, 50, 100, and 200 p.p.m., more eggs were laid in the control substrate, indicating oviposition repellency. Analysis of the larval cuticle by GC and mass spectrometry confirmed the presence of n‐heneicosane in the cuticles of Ae. albopictus larvae. The species‐specific role of n‐heneicosane as an oviposition pheromone in Ae. aegypti and its significance as a behaviour modifier of Ae. albopictus in breeding sites is discussed.     

Rapid evolution of Medicago polymorpha during invasion shifts interactions with the soybean looper

The Enemy Release Hypothesis posits that invasion of novel habitats can be facilitated by the absence of coevolved herbivores. However, a new environment and interactions with unfamiliar herbivores may impose selection on invading plants for traits that reduce their attractiveness to herbivores or for enhanced defenses compared to native host plants, leading to a pattern similar to enemy release but driven by evolutionary change rather than ecological differences. The Shifting Defense Hypothesis posits that plants in novel habitats will shift from specialized defense mechanisms to defense mechanisms effective against generalist herbivores in the new range. We tested these ideas by comparing herbivore preference and performance of native (Eurasia)‐ and invasive (New World)‐range Medicago polymorpha, using a generalist herbivore, the soybean looper, that co‐occurs with M. polymorpha in its New World invaded range. We found that soybean loopers varied in preference and performance depending on host genotype and that overall the herbivore preferred to consume plant genotypes from naïve populations from Eurasia. This potentially suggests that range expansion of M. polymorpha into the New World has led to rapid evolution of a variety of traits that have helped multiple populations become established, including those that may allow invasive populations to resist herbivory. Thus, enemy release in a novel range can occur through rapid evolution by the plant during invasion, as predicted by the Shifting Defense Hypothesis, rather than via historical divergence.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.