Landscape genetics indicate recently increased habitat fragmentation in African forest‐associated chafers

Today, indigenous forests cover less than 0.6% of South Africa's land surface and are highly fragmented. Most forest relicts are very small and typically occur in fire‐protected gorges along the eastern Great Escarpment. Yet, they hold a unique and valuable fauna with high endemism and ancient phylogenetic lineages, fostered by long‐term climatic stability and complex microclimates. Despite numerous studies on southern African vegetation cover, the current state of knowledge about the natural extension of indigenous forests is rather fragmentary. We use an integrated approach of population‐level phylogeography and climatic niche modeling of forest‐associated chafer species to assess connectivity and extent of forest habitats since the last glacial maximum. Current and past species distribution models ascertained potential fluctuations of forest distribution and supported a much wider potential current extension of forests based on climatic data. Considerable genetic admixture of mitochondrial and nuclear DNA among many populations and an increase in mean population mutation rate in Extended Bayesian Skyline Plots of all species indicated more extended or better connected forests in the recent past (<5 kya). Genetic isolation of certain populations, as revealed by population differentiation statistics (), as well as landscape connectivity statistics and habitat succession scenarios suggests considerable loss of habitat connectivity. As major anthropogenic influence is likely, conservational actions need to be considered.     

Scope for genetic rescue of an endangered subspecies though re‐establishing natural gene flow with another subspecies

Genetic diversity is positively linked to the viability and evolutionary potential of species but is often compromised in threatened taxa. Genetic rescue by gene flow from a more diverse or differentiated source population of the same species can be an effective strategy for alleviating inbreeding depression and boosting evolutionary potential. The helmeted honeyeater Lichenostomus melanops cassidix is a critically endangered subspecies of the common yellow‐tufted honeyeater. Cassidix has declined to a single wild population of ~130 birds, despite being subject to intensive population management over recent decades. We assessed changes in microsatellite diversity in cassidix over the last four decades and used population viability analysis to explore whether genetic rescue through hybridization with the neighbouring Lichenostomus melanops gippslandicus subspecies constitutes a viable conservation strategy. The contemporary cassidix population is characterized by low genetic diversity and effective population size (N_e < 50), suggesting it is vulnerable to inbreeding depression and will have limited capacity to evolve to changing environments. We find that gene flow from gippslandicus to cassidix has declined substantially relative to pre‐1990 levels and argue that natural levels of gene flow between the two subspecies should be restored. Allowing gene flow (~4 migrants per generation) from gippslandicus into cassidix (i.e. genetic rescue), in combination with continued annual release of captive‐bred cassidix (i.e. demographic rescue), should lead to positive demographic and genetic outcomes. Although we consider the risk of outbreeding depression to be low, we recommend that genetic rescue be managed within the context of the captive breeding programme, with monitoring of outcomes.     

Genomic signatures of evolution in Nautilus—An endangered living fossil

Living fossils are survivors of previously more diverse lineages that originated millions of years ago and persisted with little morphological change. Therefore, living fossils are model organisms to study both long‐term and ongoing adaptation and speciation processes. However, many aspects of living fossil evolution and their persistence in the modern world remain unclear. Here, we investigate three major aspects of the evolutionary history of living fossils: cryptic speciation, population genetics and effective population sizes, using members of the genera Nautilus and Allonautilus as classic examples of true living fossils. For this, we analysed genomewide ddRAD‐Seq data for all six currently recognized nautiloid species throughout their distribution range. Our analyses identified three major allopatric Nautilus clades: a South Pacific clade, subdivided into three subclades with no signs of admixture between them; a Coral Sea clade, consisting of two genetically distinct populations with significant admixture; and a widespread Indo‐Pacific clade, devoid of significant genetic substructure. Within these major clades, we detected five Nautilus groups, which likely correspond to five distinct species. With the exception of Nautilus macromphalus, all previously described species are at odds with genomewide data, testifying to the prevalence of cryptic species among living fossils. Detailed F_ST analyses further revealed significant genome‐wide and locus‐specific signatures of selection between species and differentiated populations, which is demonstrated here for the first time in a living fossil. Finally, approximate Bayesian computation (ABC) simulations suggest large effective population sizes, which may explain the low levels of population differentiation commonly observed in living fossils.     

Large‐scale adaptive divergence in Boechera fecunda,an endangered wild relative of Arabidopsis

Many biological species are threatened with extinction because of a number of factors such as climate change and habitat loss, and their preservation depends on an accurate understanding of the extent of their genetic variability within and among populations. In this study, we assessed the genetic divergence of five quantitative traits in 10 populations of an endangered cruciferous species, Boechera fecunda, found in only several populations in each of two geographic regions (WEST and EAST) in southwestern Montana. We analyzed variation in quantitative traits, neutral molecular markers, and environmental factors and provided evidence that despite the restricted geographical distribution of this species, it exhibits a high level of genetic variation and regional adaptation. Conservation efforts therefore should be directed to the preservation of populations in each of these two regions without attempting transplantation between regions. Heritabilities and genetic coefficients of variation estimated from nested ANOVAs were generally high for leaf and rosette traits, although lower (and not significantly different from 0) for water‐use efficiency. Measures of quantitative genetic differentiation, Q_ST, were calculated for each trait from each pair of populations. For three of the five traits, these values were significantly higher between regions compared with those within regions (after adjustment for neutral genetic variation, F_ST). This suggested that natural selection has played an important role in producing regional divergence in this species. Our analysis also revealed that the B. fecunda populations appear to be locally adapted due, at least in part, to differences in environmental conditions in the EAST and WEST regions.     

Low‐coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp

Population divergence and gene flow are key processes in evolution and ecology. Model‐based analysis of genome‐wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome‐level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.     

Landscape genomic approach to detect selection signatures in locally adapted Brazilian swine genetic groups

Samples of 191 animals from 18 different Brazilian locally adapted swine genetic groups were genotyped using Illumina Porcine SNP60 BeadChip in order to identify selection signatures related to the monthly variation of Brazilian environmental variables. Using BayeScan software, 71 SNP markers were identified as F_ST outliers and 60 genotypes (58 markers) were found by Samβada software in 371 logistic models correlated with 112 environmental variables. Five markers were identified in both methods, with a Kappa value of 0.073 (95% CI: 0.011–0.134). The frequency of these markers indicated a clear north–south country division that reflects Brazilian environmental differences in temperature, solar radiation, and precipitation. Global spatial territory correlation for environmental variables corroborates this finding (average Moran's I = 0.89, range from 0.55 to 0.97). The distribution of alleles over the territory was not strongly correlated with the breed/genetic groups. These results are congruent with previous mtDNA studies and should be used to direct germplasm collection for the National gene bank.     

Using genomic information for management planning of an endangered perennial,Viola uliginosa

Species occupying habitats subjected to frequent natural and/or anthropogenic changes are a challenge for conservation management. We studied one such species, Viola uliginosa, an endangered perennial wetland species typically inhabiting sporadically flooded meadows alongside rivers/lakes. In order to estimate genomic diversity, population structure, and history, we sampled five sites in Finland, three in Estonia, and one each in Slovenia, Belarus, and Poland using genomic SNP data with double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq). We found monophyletic populations, high levels of inbreeding (mean population F_SNP = 0.407–0.945), low effective population sizes (N_e = 0.8–50.9), indications of past demographic expansion, and rare long‐distance dispersal. Our results are important in implementing conservation strategies for V. uliginosa, which should include founding of seed banks, ex situ cultivations, and reintroductions with individuals of proper origin, combined with continuous population monitoring and habitat management.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

Genetic depletion at adaptive but not neutral loci in an endangered bird species

Many endangered species suffer from the loss of genetic diversity, but some populations may be able to thrive even if genetically depleted. To investigate the underlying genetic processes of population bottlenecks, we apply an innovative approach for assessing genetic diversity in the last known population of the endangered Pale‐headed Brushfinch (Atlapetes pallidiceps) in Ecuador. First, we measure genetic diversity at eleven neutral microsatellite loci and adaptive SNP variation in five Toll‐like receptor (TLR) immune system genes. Bottleneck tests confirm genetic drift as the main force shaping genetic diversity in this species and indicate a 99 % reduction in population size dating back several hundred years. Second, we compare contemporary microsatellite diversity with historic museum samples of A. pallidiceps, finding no change in genetic diversity. Third, we compare genetic diversity in the Pale‐headed Brushfinch with two co‐occurring‐related brushfinch species (Atlapetes latinuchus, Buarremon torquatus), finding a reduction of up to 91% diversity in the immune system genes but not in microsatellites. High TLR diversity is linked to decreased survival probabilities in A. pallidiceps. Low TLR diversity is thus probably an adaptation to the specific selection regime within its currently very restricted distribution (approximately 200 ha), but could severely restrict the adaptive potential of the species in the long run. Our study illustrates the importance of investigating both neutral and adaptive markers to assess the effect of population bottlenecks and for recommending specific management plans in endangered species.     

Population genomic analysis uncovers African and European admixture in Drosophila melanogaster populations from the south‐eastern United States and Caribbean Islands

Drosophila melanogaster is postulated to have colonized North America in the past several 100 years in two waves. Flies from Europe colonized the east coast United States while flies from Africa inhabited the Caribbean, which if true, make the south‐east US and Caribbean Islands a secondary contact zone for African and European D. melanogaster. This scenario has been proposed based on phenotypes and limited genetic data. In our study, we have sequenced individual whole genomes of flies from populations in the south‐east US and Caribbean Islands and examined these populations in conjunction with population sequences from the west coast US, Africa, and Europe. We find that west coast US populations are closely related to the European population, likely reflecting a rapid westward expansion upon first settlements into North America. We also find genomic evidence of African and European admixture in south‐east US and Caribbean populations, with a clinal pattern of decreasing proportions of African ancestry with higher latitude. Our genomic analysis of D. melanogaster populations from the south‐east US and Caribbean Islands provides more evidence for the Caribbean Islands as the source of previously reported novel African alleles found in other east coast US populations. We also find the border between the south‐east US and the Caribbean island to be the admixture hot zone where distinctly African‐like Caribbean flies become genomically more similar to European‐like south‐east US flies. Our findings have important implications for previous studies examining the generation of east coast US clines via selection.     

Exploring a Pool‐seq‐only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool‐seq data to generate a de novo genome assembly for mining exons, upon which Pool‐seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual‐based single nucleotide polymorphisms [SNPs]) and from mapping the Pool‐seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F_ST) between the two introduced populations exceeds that of the naturally sympatric populations (F_ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high‐quality reference assembly from a divergent species. We conclude that the Pool‐seq‐only approach can be suitable for detecting and quantifying genome‐wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.     

Genomic signatures of rapid adaptive evolution in the bluespotted cornetfish,a Mediterranean Lessepsian invader

Biological invasions are increasingly creating ecological and economical problems both on land and in aquatic environments. For over a century, the Mediterranean Sea has steadily been invaded by Indian Ocean/Red Sea species (called Lessepsian invaders) via the Suez Canal, with a current estimate of ~450 species. The bluespotted cornetfish, Fistularia commersonii, considered a ‘Lessepsian sprinter’, entered the Mediterranean in 2000 and by 2007 had spread through the entire basin from Israel to Spain. The situation is unique and interesting both because of its unprecedented rapidity and by the fact that it took this species c. 130 years to immigrate into the Mediterranean. Using genome scans, with restriction site‐associated DNA (RAD) sequencing, we evaluated neutral and selected genomic regions for Mediterranean vs. Red Sea cornetfish individuals. We found that few fixed neutral changes were detectable among populations. However, almost half of the genes associated with the 47 outlier loci (potentially under selection) were related to disease resistance and osmoregulation. Due to the short time elapsed from the beginning of the invasion to our sampling, we interpret these changes as signatures of rapid adaptation that may be explained by several mechanisms including preadaptation and strong local selection. Such genomic regions are therefore good candidates to further study their role in invasion success.     

Roads to isolation: Similar genomic history patterns in two species of freshwater crabs with contrasting environmental tolerances and range sizes

Freshwater species often show high levels of endemism and risk of extinction owing to their limited dispersal abilities. This is exemplified by the stenotopic freshwater crab, Johora singaporensis which is one of the world's 100 most threatened species, and currently inhabits less than 0.01 km² of five low order hill streams within the highly urbanized island city‐state of Singapore. We compared populations of J. singaporensis with that of the non‐threatened, widespread, abundant, and eurytopic freshwater crab, Parathelphusa maculata, and found surprisingly high congruence between their population genomic histories. Based on 2,617 and 2,470 genome‐wide SNPs mined via the double‐digest restriction‐associated DNA sequencing method for ~90 individuals of J. singaporensis and P. maculata, respectively, the populations are strongly isolated (F_ST = 0.146–0.371), have low genetic diversity for both species (also for COI), and show signatures of recent genetic bottlenecks. The most genetically isolated populations for both species are separated from other populations by one of the oldest roads in Singapore. These results suggest that anthropogenic developments may have impacted stream‐dependent species in a uniform manner, regardless of ubiquity, habitat preference, or dispersal modes of the species. While signs of inbreeding were not detected for the critically endangered species, the genetic distinctiveness and low diversity of the populations call for genetic rescue and connecting corridors between the remaining fragments of the natural habitat.     

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.     

Gene flow and population structure of a solitary top carnivore in a human‐dominated landscape

While African leopard populations are considered to be continuous as demonstrated by their high genetic variation, the southernmost leopard population exists in the Eastern and Western Cape, South Africa, where anthropogenic activities may be affecting this population's structure. Little is known about the elusive, last free‐roaming top predator in the region and this study is the first to report on leopard population structuring using nuclear DNA. By analyzing 14 microsatellite markers from 40 leopard tissue samples, we aimed to understand the populations' structure, genetic distance, and gene flow (Nm). Our results, based on spatially explicit analysis with Bayesian methods, indicate that leopards in the region exist in a fragmented population structure with lower than expected genetic diversity. Three population groups were identified, between which low to moderate levels of gene flow were observed (Nm 0.5 to 3.6). One subpopulation exhibited low genetic differentiation, suggesting a continuous population structure, while the remaining two appear to be less connected, with low emigration and immigration between these populations. Therefore, genetic barriers are present between the subpopulations, and while leopards in the study region may function as a metapopulation, anthropogenic activities threaten to decrease habitat and movement further. Our results indicate that the leopard population may become isolated within a few generations and suggest that management actions should aim to increase habitat connectivity and reduce human–carnivore conflict. Understanding genetic diversity and connectivity of populations has important conservation implications that can highlight management of priority populations to reverse the effects of human‐caused extinctions.     

Population genomic analyses from low‐coverage RAD‐Seq data: a case study on the non‐model cucurbit bottle gourd

Restriction site‐associated DNA sequencing (RAD‐Seq), a next‐generation sequencing‐based genome ‘complexity reduction’ protocol, has been useful in population genomics in species with a reference genome. However, the application of this protocol to natural populations of genomically underinvestigated species, particularly under low‐to‐medium sequencing depth, has not been well justified. In this study, a Bayesian method was developed for calling genotypes from an F₂ population of bottle gourd [Lagenaria siceraria (Mol.) Standl.] to construct a high‐density genetic map. Low‐depth genome shotgun sequencing allowed the assembly of scaffolds/contigs comprising approximately 50% of the estimated genome, of which 922 were anchored for identifying syntenic regions between species. RAD‐Seq genotyping of a natural population comprising 80 accessions identified 3226 single nuclear polymorphisms (SNPs), based on which two sub‐gene pools were suggested for association with fruit shape. The two sub‐gene pools were moderately differentiated, as reflected by the Hudson's F_ST value of 0.14, and they represent regions on LG7 with strikingly elevated F_ST values. Seven‐fold reduction in heterozygosity and two times increase in LD (r²) were observed in the same region for the round‐fruited sub‐gene pool. Outlier test suggested the locus LX3405 on LG7 to be a candidate site under selection. Comparative genomic analysis revealed that the cucumber genome region syntenic to the high F_ST island on LG7 harbors an ortholog of the tomato fruit shape gene OVATE. Our results point to a bright future of applying RAD‐Seq to population genomic studies for non‐model species even under low‐to‐medium sequencing efforts. The genomic resources provide valuable information for cucurbit genome research.     

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

The adaptive genomic landscape of beak morphology in Darwin's finches

Beak shape in Darwin's ground finches (Geospiza) is emblematic of natural selection and adaptive radiation, yet our understanding of the genetic basis of beak shape variation, and thus the genetic target of natural selection, is still evolving. Here we reveal the genomic architecture of beak shape variation using genomewide comparisons of four closely related and hybridizing species across 13 islands subject to parallel natural selection. Pairwise contrasts among species were used to identify a large number of genomic loci that are consistently related to species differences across a complex landscape. These loci are associated with hundreds of genes that have enriched GO categories significantly associated with development. One genomic region of particular interest is a section of Chromosome 1A with many candidate genes and increased linkage. The distinct, pointed beak shape of the cactus finch is linked to an excess of intermediate frequency alleles and increased heterozygosity in significant SNPs, but not across the rest of the genome. Alleles associated with pointier beaks among species were associated with pointier‐beaked populations within each species, thus establishing a common basis for natural selection, species divergence and adaptive radiation. The adaptive genomic landscape for Darwin's finches mirrors theoretical expectations based on morphological variation. The implication that a large number of genes are actively maintained to facilitate beak variation across parallel populations with documented interspecies admixture challenges our understanding of evolutionary processes in the wild.