Small populations of endangered species risk losing already eroded genetic diversity, important for adaptive potential, through the effects of genetic drift. The magnitude of drift can be mitigated by maximising the effective population size, as is the goal of genetic management strategies. Different mating systems, specifically those leading to reproductive skew, exacerbate genetic drift by distorting contributions. In the absence of an active management strategy, reproductive skew will have long-term effects on the genetic composition of a population, particularly where admixture is present. Here we examine the contrasting effects of conservation management strategies in two ex situ populations of the Critically Endangered eastern black rhino (Diceros bicornis michaeli), one managed as a semi-wild population in South Africa (SAx), and one managed under a mean-kinship breeding strategy in European zoos. We use molecular data to reconstruct pedigrees for both populations and validate the method using the zoo studbook. Using the reconstructed pedigree and studbook we show there is male sex-specific skew in both populations. However, the zoo’s mean-kinship breeding strategy effectively reduces reproductive skew in comparison to a semi-wild population with little genetic management. We also show that strong male reproductive skew in SAx has resulted in extensive admixture, which may require a re-evaluation of the population’s original intended role in the black rhino meta-population. With a high potential for admixture in many ex situ populations of endangered species, molecular and pedigree data remain vital tools for populations needing to balance drift and selection.

Understanding the demographic history and genetic make-up of colonizing species is critical for inferring population sources and colonization routes. This is of main interest for designing accurate control measures in areas newly colonized by vector species of economically important pathogens. The biting midge Culicoides imicola is a major vector of orbiviruses to livestock. Historically, the distribution of this species was limited to the Afrotropical region. Entomological surveys first revealed the presence of C. imicola in the south of the Mediterranean basin by the 1970s. Following recurrent reports of massive bluetongue outbreaks since the 1990s, the presence of the species was confirmed in northern areas. In this study, we addressed the chronology and processes of C. imicola colonization in the Mediterranean basin. We characterized the genetic structure of its populations across Mediterranean and African regions using both mitochondrial and nuclear markers, and combined phylogeographical analyses with population genetics and approximate Bayesian computation. We found a west/east genetic differentiation between populations, occurring both within Africa and within the Mediterranean basin. We demonstrated that three of these groups had experienced demographic expansions in the Pleistocene, probably because of climate changes during this period. Finally, we showed that C. imicola could have colonized the Mediterranean basin in the Late Pleistocene or Early Holocene through a single event of introduction; however, we cannot exclude the hypothesis involving two routes of colonization. Thus, the recent bluetongue outbreaks are not linked to C. imicola colonization event, but rather to biological changes in the vector or the virus.

Recent advances in sequencing allow population-genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction-site-associated DNA sequence (RAD-seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well-characterized single nucleotide polymorphism (SNP) data set from 21 three-spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single-outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population-genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population-demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population-genomic data set, making it especially valuable for nonmodel species.

Southeast Asia harbours abundant biodiversity, hypothesized to have been generated by Pliocene and Pleistocene climatic and environmental change. Vicariance between the island of Borneo, the remaining Indonesian archipelago and mainland Southeast Asia caused by elevated sea levels during interglacial periods has been proposed to lead to diversification in the littoral zone mosquito Anopheles (Cellia) sundaicus (Rodenwaldt) sensu lato. To test this biogeographical hypothesis, we inferred the population history and assessed gene flow of A. sundaicus s.l. sampled from 18 populations across its pan-Asian species range, using sequences from mitochondrial cytochrome c oxidase subunit 1 (CO1), the internal transcribed spacer 2 (ITS2) and the mannose phosphate isomerase (Mpi) gene. A hypothesis of ecological speciation for A. sundaicus involving divergent adaptation to brackish and freshwater larval habitats was also previously proposed, based on a deficiency of heterozygotes for Mpi allozyme alleles in sympatry. This hypothesis was not supported by Mpi sequence data, which exhibited no fixed differences between brackish and freshwater larval habitats. Mpi and CO1 supported the presence of up to eight genetically distinct population groupings. Counter to the hypothesis of three allopatric species, divergence was often no greater between Borneo, Sumatra/Java and the Southeast Asian mainland than it was between genetic groupings within these landmasses. An isolation-with-migration (IM) model indicates recurrent gene flow between the current major landmasses. Such gene flow would have been possible during glacial periods when the current landmasses merged, presenting opportunities for dispersal along expanding and contracting coastlines. Consequently, Pleistocene climatic variation has proved a homogenizing, rather than diversifying, force for A. sundaicus diversity.

Southeast Asia is one of the world’s richest regions in terms of biodiversity. An understanding of the distribution of diversity and the factors shaping it is lacking, yet essential for identifying conservation priorities for the region’s highly threatened biodiversity. Here we take a large scale comparative approach, combining data from nine forest associated Anopheles mosquito species and using statistical phylogeographic methods to disentangle the effects of environmental history, species specific ecology, and random coalescent effects. Spatially explicit modelling of Pleistocene demographic history supports a common influence of environmental events in shaping the genetic diversity of all species examined, despite differences in species' mtDNA gene trees. Populations were periodically restricted to allopatric northeastern and northwestern refugia, most likely due to Pleistocene forest fragmentation. Subsequent southwards post-glacial recolonisation is supported by a north-south gradient of decreasing genetic diversity. Repeated allopatric fragmentation and recolonisation has led to the formation of deeply divergent geographical lineages within four species and a suture zone where these intraspecific lineages meet along the Thai-Myanmar border. A common environmental influence for this divergence was further indicated by strong support for simultaneous divergence within the same four species, dating to approximately 900 kya. Differences in the geographical structuring of genetic diversity between species are likely the result of varying species’ biology. The findings have important implications for conservation planning; if the refugial regions and suture zone identified here are shared by other forest taxa, the unique and high levels of genetic diversity they house will make these areas conservation priorities.

Tropical forests have undergone repeated fragmentation and expansion during Pleistocene glacial and interglacial periods, respectively. The effects of this repeated forest fragmentation in driving vicariance in tropical taxa have been well studied. However, relatively little is known about how often this process results in allopatric speciation, since it may be inhibited by recurrent gene flow during repeated secondary contact, or to what extent Pleistocene-dated speciation results from ecological specialisation in the face of gene flow. Here, divergence times and gene flow between three closely-related mosquito species of the Anopheles dirus species complex endemic to the forests of Southeast Asia, are inferred using coalescent based Bayesian analysis. An Isolation with Migration model is applied to sequences of two mitochondrial and three nuclear genes, and 11 microsatellites. The divergence of An. scanloni has occurred despite unidirectional nuclear gene flow from this species into An. dirus. The inferred asymmetric gene flow may result from the unique evolutionary adaptation of An. scanloni to limestone karst habitat, and therefore the fitness advantage of this species over An. dirus in regions of sympatry. Mitochondrial introgression has led to the complete replacement of An. dirus haplotypes with those of An. baimaii through a recent (~62 kya) selective sweep. Speciation of An. baimaii and An. dirus is inferred to have involved allopatric divergence throughout much of the Pleistocene. Secondary contact and bidirectional gene flow has occurred only within the last 100,000 years, by which time the process of allopatric speciation seems to have been largely completed.