Sequence capture and next‐generation sequencing of ultraconserved elements in a large‐genome salamander

Amidst the rapid advancement in next‐generation sequencing (NGS) technology over the last few years, salamanders have been left behind. Salamanders have enormous genomes—up to 40 times the size of the human genome—and this poses challenges to generating NGS data sets of quality and quantity similar to those of other vertebrates. However, optimization of laboratory protocols is time‐consuming and often cost prohibitive, and continued omission of salamanders from novel phylogeographic research is detrimental to species facing decline. Here, we use a salamander endemic to the southeastern United States, Plethodon serratus, to test the utility of an established protocol for sequence capture of ultraconserved elements (UCEs) in resolving intraspecific phylogeographic relationships and delimiting cryptic species. Without modifying the standard laboratory protocol, we generated a data set consisting of over 600 million reads for 85 P. serratus samples. Species delimitation analyses support recognition of seven species within P. serratus sensu lato, and all phylogenetic relationships among the seven species are fully resolved under a coalescent model. Results also corroborate previous data suggesting nonmonophyly of the Ouachita and Louisiana regions. Our results demonstrate that established UCE protocols can successfully be used in phylogeographic studies of salamander species, providing a powerful tool for future research on evolutionary history of amphibians and other organisms with large genomes.     

Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

An empirical comparison of character‐based and coalescent‐based approaches to species delimitation in a young avian complex

The process of discovering species is a fundamental responsibility of systematics. Recently, there has been a growing interest in coalescent‐based methods of species delimitation aimed at objectively identifying species early in the divergence process. However, few empirical studies have compared these new methods with character‐based approaches for discovering species. In this study, we applied both a character‐based and a coalescent‐based approaches to delimit species in a closely related avian complex, the light‐vented/Taiwan bulbul (Pycnonotus sinensis/Pycnonotus taivanus). Population aggregation analyses of plumage, mitochondrial and 13 nuclear intron character data sets produced conflicting species hypotheses with plumage data suggesting three species, mitochondrial data suggesting two species, and nuclear intron data suggesting one species. Such conflict is expected among recently diverged species, and by integrating all sources of data, we delimited three species verified with independently congruent character evidence as well as a more weakly supported fourth species identified by a single character. Attempts to validate species hypothesis using Bayesian Phylogenetics and Phylogeography (BPP), a coalescent‐based method of species delimitation, revealed several issues that can seemingly affect statistical support for species recognition. We found that θ priors had a dramatic impact on speciation probabilities, with lower values consistently favouring splitting and higher values consistently favouring lumping. More resolved guide trees also resulted in overall higher speciation probabilities. Finally, we found suggestive evidence that BPP is sensitive to the divergent effects of nonrandom mating caused by intraspecific processes such as isolation‐with‐distance, and therefore, BPP may not be a conservative method for delimiting independently evolving population lineages. Based on these concerns, we questioned the reliability of BPP results and based our conclusions about species limits exclusively on character data.     

Divergence,gene flow,and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison‐dart frogs

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between ‘leapfroging’ taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome‐averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.     

Gut passage effect of the introduced red‐whiskered bulbul (Pycnonotus jocosus) on germination of invasive plant species in Mauritius

In Mauritius, many of the worst invasive plant species have fleshy fruits and rely on animals for dispersal. The introduced red‐whiskered bulbul (Pycnonotus jocosus) feeds on many fleshy‐fruited species, and often moves from invaded and degraded habitats into higher quality native forests, thus potentially acting as a mediator of continued plant invasion into these areas. Furthermore, gut passage may influence seed germination. To investigate this, we fed fleshy fruits of two invasive plant species, Ligustrum robustum and Clidemia hirta, to red‐whiskered bulbuls. Gut passage times of seeds were recorded. Gut‐passed seeds were sown and their germination rate and germination success compared with that of hand‐cleaned seeds, as well as that of seeds in whole fruits. Gut passage and hand‐cleaning had significant positive effects on germination of both species. Gut‐passed seeds of both C. hirta and L. robustum germinated faster than hand‐cleaned seeds. However, for L. robustum, this was only true when compared with hand‐cleaned seeds with intact endocarp; when compared with hand‐cleaned seeds without endocarp, there was no difference. For overall germination success, there was a positive effect of gut passage for C. hirta, but not for L. robustum. For both C. hirta and L. robustum, no seeds in intact fruits geminated, suggesting that removal of pulp is essential for germination. Our results suggest that, first, the initial invasion of native forests in Mauritius may not have happened so rapidly without efficient avian seed dispersers like the red‐whiskered bulbul. Second, the bulbul is likely to be a major factor in the continued re‐invasion of C. hirta and L. robustum into weeded and restored conservation management areas.     

Evolutionary history of rat‐borne Bartonella: the importance of commensal rats in the dissemination of bacterial infections globally

Emerging pathogens that originate from invasive species have caused numerous significant epidemics. Some bacteria of genus Bartonella are rodent‐borne pathogens that can cause disease in humans and animals alike. We analyzed gltA sequences of 191 strains of rat‐associated bartonellae from 29 rodent species from 17 countries to test the hypotheses that this bacterial complex evolved and diversified in Southeast Asia before being disseminated by commensal rats Rattus rattus (black rat) and Rattus norvegicus (Norway rat) to other parts of the globe. The analysis suggests that there have been numerous dispersal events within Asia and introductions from Asia to other regions, with six major clades containing Southeast Asian isolates that appear to have been dispersed globally. Phylogeographic analyses support the hypotheses that these bacteria originated in Southeast Asia and commensal rodents (R. rattus and R. norvegicus) play key roles in the evolution and dissemination of this Bartonella complex throughout the world.     

Life‐history predicts past and present population connectivity in two sympatric sea stars

Life‐history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic‐developing sea star, Parvulastra exigua, has lower levels of within‐ and among‐population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic‐developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic‐developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post‐Pleistocene introgression.     

Invasion facilitates hybridization with introgression in the Rattus rattus species complex

Biological invasions result in novel species interactions, which can have significant evolutionary impacts on both native and invading taxa. One evolutionary concern with invasions is hybridization among lineages that were previously isolated, but make secondary contact in their invaded range(s). Black rats, consisting of several morphologically very similar but genetically distinct taxa that collectively have invaded six continents, are arguably the most successful mammalian invaders on the planet. We used mitochondrial cytochrome b sequences, two nuclear gene sequences (Atp5a1 and DHFR) and nine microsatellite loci to examine the distribution of three invasive black rat lineages (Rattus tanezumi, Rattus rattus I and R. rattus IV) in the United States and Asia and to determine the extent of hybridization among these taxa. Our analyses revealed two mitochondrial lineages that have spread to multiple continents, including a previously undiscovered population of R. tanezumi in the south‐eastern United States, whereas the third lineage (R. rattus IV) appears to be confined to Southeast Asia. Analyses of nuclear DNA (both sequences and microsatellites) suggested significant hybridization is occurring among R. tanezumi and R. rattus I in the United States and also suggest hybridization between R. tanezumi and R. rattus IV in Asia, although further sampling of the latter species pair in Asia is required. Furthermore, microsatellite analyses suggest unidirectional introgression from both R. rattus I and R. rattus IV into R. tanezumi. Within the United States, introgression appears to be occurring to such a pronounced extent that we were unable to detect any nuclear genetic signal for R. tanezumi, and a similar pattern was detected in Asia.     

Herbarium specimens reveal a historical shift in phylogeographic structure of common ragweed during native range disturbance

Invasive plants provide ample opportunity to study evolutionary shifts that occur after introduction to novel environments. However, although genetic characters pre‐dating introduction can be important determinants of later success, large‐scale investigations of historical genetic structure have not been feasible. Common ragweed (Ambrosia artemisiifolia L.) is an invasive weed native to North America that is known for its allergenic pollen. Palynological records from sediment cores indicate that this species was uncommon before European colonization of North America, and ragweed populations expanded rapidly as settlers deforested the landscape on a massive scale, later becoming an aggressive invasive with populations established globally. Towards a direct comparison of genetic structure now and during intense anthropogenic disturbance of the late 19th century, we sampled 45 natural populations of common ragweed across its native range as well as historical herbarium specimens collected up to 140 years ago. Bayesian clustering analyses of 453 modern and 473 historical samples genotyped at three chloroplast spacer regions and six nuclear microsatellite loci reveal that historical ragweed's spatial genetic structure mirrors both the palaeo‐record of Ambrosia pollen deposition and the historical pattern of agricultural density across the landscape. Furthermore, for unknown reasons, this spatial genetic pattern has changed substantially in the intervening years. Following on previous work relating morphology and genetic expression between plants collected from eastern North America and Western Europe, we speculate that the cluster associated with humans’ rapid transformation of the landscape is a likely source of these aggressive invasive populations.     

Micro‐evolutionary diversification among Indian Ocean parrots: temporal and spatial changes in phylogenetic diversity as a consequence of extinction and invasion

Almost 90% of global bird extinctions have occurred on islands. The loss of endemic species from island systems can dramatically alter evolutionary trajectories of insular species biodiversity, resulting in a loss of evolutionary diversity important for species adaptation to changing environments. The Western Indian Ocean islands have been the scene of evolution for a large number of endemic parrots. Since their discovery in the 16th century, many of these parrots have become extinct or have declined in numbers. Alongside the extinction of species, a number of the Indian Ocean islands have experienced colonization by highly invasive parrots, such as the Ring‐necked Parakeet Psittacula krameri. Such extinctions and invasions can, on an evolutionary timescale, drive changes in species composition, genetic diversity and turnover in phylogenetic diversity, all of which can have important impacts on species potential for adaptation to changing environmental and climatic conditions. Using mtDNA cytochrome b data, we resolve the taxonomic placement of three extinct Indian Ocean parrots: the Rodrigues Psittacula exsul, Seychelles Psittacula wardi and Reunion Parakeets Psittacula eques. This case study quantifies how the extinction of these species has resulted in lost historical endemic phylogenetic diversity and reduced levels of species richness, and illustrates how it is being replaced by non‐endemic invasive forms such as the Ring‐necked Parakeet. Finally, we use our phylogenetic framework to identify and recommend a number of phylogenetically appropriate ecological replacements for the extinct parrots. Such replacements may be introduced once invasive forms have been cleared, to rejuvenate ecosystem function and restore lost phylogenetic diversity.     

Genetic diversity and haplotype distribution of Pachymeniopsis gargiuli sp. nov. and P. lanceolata (Halymeniales,Rhodophyta) in Korea,with notes on their non‐native distributions

The red alga Pachymeniopsis lanceolata, formerly known as Grateloupia lanceolata, is a component of the native algal flora of northeast Asia and has been introduced to European and North American waters. It has been confused with a cryptic species collected from Korea and Italy. Our analyses of rbcL, cox3 and ITS from P. lanceolata and this cryptic species has revealed two distinct entities, forming a clade, which were clearly separated from its congeners and positioned with other Asian species. Here, we describe the cryptic species as P. gargiuli sp. nov., a species that differs from others by molecular sequence and subtle anatomical characters. We hypothesize that P. gargiuli may have been recently dispersed by anthropogenic vectors, possibly at or near the same time as was P. lanceolata. Our cox3 data set revealed that one haplotype of P. gargiuli, shared between Korea and Italy, and two haplotypes of P. lanceolata, commonly occurring in Korea and USA, are invasive haplotypes. This is the first report of the utility of the mitochondrial coding cox3 sequences in red algae.     

The advantages of going large: genome‐wide SNPs clarify the complex population history and systematics of the threatened western pond turtle

As the field of phylogeography has matured, it has become clear that analyses of one or a few genes may reveal more about the history of those genes than the populations and species that are the targets of study. To alleviate these concerns, the discipline has moved towards larger analyses of more individuals and more genes, although little attention has been paid to the qualitative or quantitative gains that such increases in scale and scope may yield. Here, we increase the number of individuals and markers by an order of magnitude over previously published work to comprehensively assess the phylogeographical history of a well‐studied declining species, the western pond turtle (Emys marmorata). We present a new analysis of 89 independent nuclear SNP markers and one mitochondrial gene sequence scored for rangewide sampling of >900 individuals, and compare these to smaller‐scale, rangewide genetic and morphological analyses. Our enlarged SNP data fundamentally revise our understanding of evolutionary history for this lineage. Our results indicate that the gains from greatly increasing both the number of markers and individuals are substantial and worth the effort, particularly for species of high conservation concern such as the pond turtle, where accurate assessments of population history are a prerequisite for effective management.     

High‐throughput SNP‐genotyping analysis of the relationships among Ponto‐Caspian sturgeon species

Legally certified sturgeon fisheries require population protection and conservation methods, including DNA tests to identify the source of valuable sturgeon roe. However, the available genetic data are insufficient to distinguish between different sturgeon populations, and are even unable to distinguish between some species. We performed high‐throughput single‐nucleotide polymorphism (SNP)‐genotyping analysis on different populations of Russian (Acipenser gueldenstaedtii), Persian (A. persicus), and Siberian (A. baerii) sturgeon species from the Caspian Sea region (Volga and Ural Rivers), the Azov Sea, and two Siberian rivers. We found that Russian sturgeons from the Volga and Ural Rivers were essentially indistinguishable, but they differed from Russian sturgeons in the Azov Sea, and from Persian and Siberian sturgeons. We identified eight SNPs that were sufficient to distinguish these sturgeon populations with 80% confidence, and allowed the development of markers to distinguish sturgeon species. Finally, on the basis of our SNP data, we propose that the A. baerii‐like mitochondrial DNA found in some Russian sturgeons from the Caspian Sea arose via an introgression event during the Pleistocene glaciation.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Genetic diversity and demographic history of introduced sika deer on the Delmarva Peninsula

The introduction of non‐native species can have long‐term effects on native plant and animal communities. Introduced populations are occasionally not well understood and offer opportunities to evaluate changes in genetic structure through time and major population changes such as bottleneck and or founder events. Invasive species can often evolve rapidly in new and novel environments, which could be essential to their long‐term success. Sika deer are native to East Asia, and their introduction and establishment to the Delmarva Peninsula, USA, is poorly documented, but probably involved ≥1 founder and/or bottleneck events. We quantified neutral genetic diversity in the introduced population and compared genetic differentiation and diversity to the presumed source population from Yakushima Island, Japan, and a captive population of sika deer in Harrington, Delaware, USA. Based on the data from 10 microsatellite DNA loci, we observed reduced genetic variation attributable to founder events, support for historic hybridization events, and evidence that the population did originate from Yakushima Island stocks. Estimates of population structure through Bayesian clustering and demographic history derived from approximate Bayesian computation (ABC), were consistent with the hypothesized founder history of the introduced population in both timing and effective population size (approximately five effective breeding individuals, an estimated 36 generations ago). Our ABC results further supported a single introduction into the wild happening before sika deer spread throughout the Delmarva. We conclude that free‐ranging sika deer on Delmarva are descended from ca. five individuals introduced about 100 years ago from captive stocks of deer maintained in the United Kingdom. Free‐ranging sika deer on Delmarva have lost neutral diversity due to founder and bottleneck events, yet populations have expanded in recent decades and show no evidence of abnormalities associated with inbreeding. We suggest management practices including increasing harvest areas and specifically managing sika deer outside of Maryland.     

Influence of invasion history on rapid morphological divergence across island populations of an exotic bird

There is increasing evidence that exotic populations may rapidly differentiate from those in their native range and that differences also arise among populations within the exotic range. Using morphological and DNA‐based analyses, we document the extent of trait divergence among native North American and exotic Hawaiian populations of northern cardinal (Cardinalis cardinalis). Furthermore, using a combination of historical records and DNA‐based analyses, we evaluate the role of founder effects in producing observed trait differences. We measured and compared key morphological traits across northern cardinal populations in the native and exotic ranges to assess whether trait divergence across the Hawaiian Islands, where this species was introduced between 1929 and 1931, reflected observed variation across native phylogeographic clades in its native North America. We used and added to prior phylogenetic analyses based on a mitochondrial locus to identify the most likely native source clade(s) for the Hawaiian cardinal populations. We then used Approximate Bayesian Computation (ABC) to evaluate the role of founder effects in producing the observed differences in body size and bill morphology across native and exotic populations. We found cardinal populations on the Hawaiian Islands had morphological traits that diverged substantially across islands and overlapped the trait space of all measured native North American clades. The phylogeographic analysis identified the eastern North American clade (C. cardinalis cardinalis) as the most likely and sole native source for all the Hawaiian cardinal populations. The ABC analyses supported written accounts of the cardinal's introduction that indicate the original 300 cardinals shipped to Hawaii were simultaneously and evenly released across Hawaii, Kauai, and Oahu. Populations on each island likely experienced bottlenecks followed by expansion, with cardinals from the island of Hawaii eventually colonizing Maui unaided. Overall, our results suggest that founder effects had limited impact on morphological trait divergence of exotic cardinal populations in the Hawaiian archipelago, which instead reflect postintroduction events.     

An invasive tree facilitates the persistence of native rodents on an over‐grazed floodplain in tropical Australia

In an ecosystem under simultaneous threat from multiple alien species, one invader may buffer the impact of another. Our surveys on a remote floodplain in the Kimberley region of north western Australia show that invasive chinee apple trees (Ziziphus mauritiana) provide critical refuge habitat for native rodents (pale field rats, Rattus tunneyi). Feral horses (Equus caballus) have trampled most of the remaining floodplain, but are excluded from the area around each chinee apple tree by thorny foliage. Although chinee apple trees constituted <10% of trees along our transects, they represented >50% of trees that harboured rat burrows. The mean number of burrows under each chinee apple tree was twice as high as under most other tree species, and we trapped more than seven times as many rats under chinee apple trees as under other types of trees. The extensive burrow systems under chinee apple trees contained female as well as male rats, whereas we only captured males around the smaller burrow systems under other tree species. Our data suggest that this invasive tree plays a critical role in the persistence of pale field rat populations in this degraded ecosystem, and that managers should maintain these trees (despite their alien origins) at least until feral horses have been removed.     

Niche shifts and range expansions along cordilleras drove diversification in a high‐elevation endemic plant genus in the tropical Andes

The tropical Andes represent one of the world's biodiversity hot spots, but the evolutionary drivers generating their striking species diversity still remain poorly understood. In the treeless high‐elevation Andean environments, Pleistocene glacial oscillations and niche differentiation are frequently hypothesized diversification mechanisms; however, sufficiently densely sampled population genetic data supporting this are still lacking. Here, we reconstruct the evolutionary history of Loricaria (Asteraceae), a plant genus endemic to the Andean treeless alpine zone, based on comprehensive population‐level sampling of 289 individuals from 67 populations across the entire distribution ranges of its northern Andean species. Partly incongruent AFLP and plastid DNA markers reveal that the distinct genetic structure was shaped by a complex interplay of biogeography (spread along and across the cordilleras), history (Pleistocene glacial oscillations) and local ecological conditions. While plastid variation documents an early split or colonization of the northern Andes by at least two lineages, one of which further diversified, a major split in the AFLP data correlate with altitudinal ecological differentiation. This suggests that niche shifts may be important drivers of Andean diversification not only in forest–alpine transitions, but also within the treeless alpine zone itself. The patterns of genetic differentiation at the intraspecific level reject the hypothesized separation in spatially isolated cordilleras and instead suggest extensive gene flow among populations from distinct mountain chains. Our study highlights that leveraging highly variable markers against extensive population‐level sampling is a promising approach to address mechanisms of rapid species diversifications.     

Whole‐genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea

Anopheles mosquitoes are the vectors of several human diseases including malaria. In many malaria endemic areas, several species of Anopheles coexist, sometimes in the form of related sibling species that are morphologically indistinguishable. Determining the size and organization of Anopheles populations, and possible ongoing gene flow among them is important for malaria control and, in particular, for monitoring the spread of insecticide resistance alleles. However, these parameters have been difficult to evaluate in most Anopheles species due to the paucity of genetic data available. Here, we assess the extent of contemporary gene flow and historical variations in population size by sequencing and de novo assembling the genomes of wild‐caught mosquitoes from four species of the Anopheles punctulatus group of Papua New Guinea. Our analysis of more than 50 Mb of orthologous DNA sequences revealed no evidence of contemporary gene flow among these mosquitoes. In addition, investigation of the demography of two of the An. punctulatus species revealed distinct population histories. Overall, our analyses suggest that, despite their similarities in morphology, behaviour and ecology, contemporary sympatric populations of An. punctulatus are evolving independently.