Phylogeography of the heavily poached African common pangolin (Pholidota,Manis tricuspis) reveals six cryptic lineages as traceable signatures of Pleistocene diversification

Knowledge on faunal diversification in African rainforests remains scarce. We used phylogeography to assess (i) the role of Pleistocene climatic oscillations in the diversification of the African common pangolin (Manis tricuspis) and (ii) the utility of our multilocus approach for taxonomic delineation and trade tracing of this heavily poached species. We sequenced 101 individuals for two mitochondrial DNA (mtDNA), two nuclear DNA and one Y‐borne gene fragments (totalizing 2602 bp). We used a time‐calibrated, Bayesian inference phylogenetic framework and conducted character‐based, genetic and phylogenetic delineation of species hypotheses within African common pangolins. We identified six geographic lineages partitioned into western Africa, Ghana, the Dahomey Gap, western central Africa, Gabon and central Africa, all diverging during the Middle to Late Pleistocene. MtDNA (cytochrome b + control region) was the sole locus to provide diagnostic characters for each of the six lineages. Tree‐based Bayesian delimitation methods using single‐ and multilocus approaches gave high support for ‘species’ level recognition of the six African common pangolin lineages. Although the diversification of African common pangolins occurred during Pleistocene cyclical glaciations, causative correlation with traditional rainforest refugia and riverine barriers in Africa was not straightforward. We conclude on the existence of six cryptic lineages within African common pangolins, which might be of major relevance for future conservation strategies. The high discriminative power of the mtDNA markers used in this study should allow an efficient molecular tracing of the regional origin of African common pangolin seizures.     

The emergence of a commercial trade in pangolins from Gabon

Recent seizures of illegally held wildlife indicate a mounting global trade in pangolins involving all eight species. Seizures of illegally traded African pangolins are increasing as wild populations of Asian species decline. We investigated trade in pangolins and law enforcement efforts in Gabon; a country likely to have intact wild populations of three of the four species of African pangolin. We compared village sales and trade chains between 2002‐2003 and 2014. Hunters reported pangolins to be the most frequently requested species in 2014, and the value of pangolins had increased at every point along their trade chain. In Libreville, giant pangolin prices increased 211% and arboreal pangolin prices 73% whilst inflation rose only 4.6% over the same period. We documented a low rate of interception of illegally traded pangolins despite increased law enforcement. Surveys of potential export routes detected exports across forest borders, in conjunction with ivory, but not through public transport routes. We conclude that whilst there is clear potential and likelihood that a wild pangolin export trade is emerging from Gabon, traditional bushmeat trade chains may not be the primary supply route. We recommend adjusting conservation policies and actions to impede further development of illegal trade within and from Gabon .     

Reproductive behavior of the captive Sunda pangolin (Manis javanica Desmarest, 1822)

Understanding reproductive behavior is important for the conservation of endangered species, but research on the reproductive behavior of Sunda pangolins (Manis javanica Desmarest, 1822) is still very scarce. In this study, we used focal animal sampling and all-occurrence recording by an infrared monitor to observe the behaviors of two Sunda pangolins during a 5-day breeding period at the Pangolin Research Base for Artificial Rescue and Conservation Breeding of South China Normal University (PRB-SCNU). The behavioral characteristics and breeding strategies were analyzed, and the results were discussed together with information on other captive Sunda pangolin pairs at PRB-SCNU. The results found that there was no obvious estrus behavior in the captive female, while the male could exhibit sexual excitement and courtship behavior after a brief introduction period. Repeated copulation continued over many days after the female accepted the courtship. The average duration of copulatory behavior was 248.9 ± 148.7 s (n = 25), and all copulation occurred between 20:00 and 08:00 hr in the natural day–night photoperiod. The mating position of Sunda pangolins was lateral–ventral and was classified as the ninth or eleventh pattern under both Dewsbury's and Dixson's classification systems. This study can provide scientific guidance for the captive breeding and management of Sunda pangolins and other pangolin species, which is of great significance for ex situ conservation tactics.     

Isolation and characterization of microsatellite markers in pangolins (Mammalia,Pholidota, Manis spp.)

Thirty‐four polymorphic dinucleotide microsatellite loci were developed in the Malayan pangolin Manis javanica. Of the 34 markers, 32 and 18 were also amplified, respectively, in the Chinese pangolin (Manis pentadactyla) and the African tree pangolin (Manis tricuspis). Analysis of 24 Malayan, 12 Chinese and 2 African tree pangolins showed high levels of variability (heterozygosity ranging from 0.321 to 0.708). These are the first available microsatellite markers in Pholidota and will be an invaluable tool for evolutionary and conservation genetic studies in pangolins.     

Nuclear and chloroplast DNA phylogeography reveals Pleistocene divergence and subsequent secondary contact of two genetic lineages of the tropical rainforest tree species Shorea leprosula (Dipterocarpaceae) in South‐East Asia

Tropical rainforests in South‐East Asia have been affected by climatic fluctuations during past glacial eras. To examine how the accompanying changes in land areas and temperature have affected the genetic properties of rainforest trees in the region, we investigated the phylogeographic patterns of a widespread dipterocarp species, Shorea leprosula. Two types of DNA markers were used: expressed sequence tag‐based simple sequence repeats and chloroplast DNA (cpDNA) sequence variations. Both sets of markers revealed clear genetic differentiation between populations in Borneo and those in the Malay Peninsula and Sumatra (Malay/Sumatra). However, in the south‐western part of Borneo, genetic admixture of the lineages was observed in the two marker types. Coalescent simulation based on cpDNA sequence variation suggested that the two lineages arose 0.28–0.09 million years before present and that following their divergence migration from Malay/Sumatra to Borneo strongly exceeded migration in the opposite direction. We conclude that the genetic structure of S. leprosula was largely formed during the middle Pleistocene and was subsequently modified by eastward migration across the subaerially exposed Sunda Shelf.     

Phylogeographic structure in mitochondrial DNA of a South-east Asian freshwater fish, Hemibagrus nemurus (Siluroidei; Bagridae) and Pleistocene sea-level changes on the Sunda shelf

We postulated that the biogeographical history of South-east Asia contributed to extensive admixture during Pleistocene low sea levels of genetic groups of an obligate freshwater fish (the river catfish, Hemibagrus nemurus) isolated during periods of high sea levels. During Pleistocene glacial maxima, the sea level was lower than at present and the islands of the Sunda shelf (Sumatra, Borneo and Java) and the Asian mainland were connected by lowlands traversed by rivers. Restriction fragment length polymorphisms in mitochondrial DNA were documented for 140 putative H. nemurus analysed from 13 sampling sites resulting in the definition of 35 haplotypes. The high level of haplotype differentiation (mean P × 100 = 2.22, SD = 1.33) indicates that the subdivision of the ancestral H. netnurus group was extensive and probably occurred early in the Pleistocene. The occurrence of some genetically divergent groups of the H. netnurus complex occurring in sympatry in widely separated locations supports the proposition that low sea levels aided the dispersion and mingling of genetic groups. Based on both genetic and morphological evidence, the main H. nemurus line gave rise to three regional groups: (1) a morphologically distinct ‘Indochinese’ group composed of two mtDNA clades overlapping in east peninsular Malaysia; (2) a ‘Sundaic’ group composed of various lineages of differing morphology and genetic identity; (3) a genetically distinct ‘Sarawak’ group in west Borneo, similar in morphology to the ‘Sundaic’ and ‘Indochinese’ groups, but including a small, golden colour morph as a distinct dade. The morphologically similar Sundaic forms from west Java, Sumatra and west Borneo show some degree of genetic divergence, but their phylogenetic relationships are poorly resolved. The most genetically and morphologically distinct Sundaic dade, assigned to H. hoevenii, colonized the Kapuas river (west Borneo), east Sumatra and south peninsular Malaysia. Contrary to our original hypothesis and present biogeographical theory, little exchange of genetic groups has apparently occurred between the mainland and the Sunda Islands during recent glaciations.     

2010-2020年中华穿山甲在中国的发现记录及保护现状

中国是中华穿山甲(Manis pentadactyla)历史分布区面积和野生种群数量最大的国家.中华穿山甲曾广泛见于我国长江以南各省,但20世纪中期以来,由于其甲片被作为贵重的中药材原料,加之地下野味市场的需求,大量非法捕猎使得我国野生穿山甲经历了剧烈的分布区缩减和种群下降.目前,中华穿山甲已被列为我国一级重点保护野生...     

Genome-wide genetic diversity,population structure and admixture analysis in African and Asian cattle breeds

Knowledge about genetic diversity and population structure is useful for designing effective strategies to improve the production, management and conservation of farm animal genetic resources. Here, we present a comprehensive genome-wide analysis of genetic diversity, population structure and admixture based on 244 animals sampled from 10 cattle populations in Asia and Africa and genotyped for 69 903 autosomal single-nucleotide polymorphisms (SNPs) mainly derived from the indicine breed. Principal component analysis, STRUCTURE and distance analysis from high-density SNP data clearly revealed that the largest genetic difference occurred between the two domestic lineages (taurine and indicine), whereas Ethiopian cattle populations represent a mosaic of the humped zebu and taurine. Estimation of the genetic influence of zebu and taurine revealed that Ethiopian cattle were characterized by considerable levels of introgression from South Asian zebu, whereas Bangladeshi populations shared very low taurine ancestry. The relationships among Ethiopian cattle populations reflect their history of origin and admixture rather than phenotype-based distinctions. The high within-individual genetic variability observed in Ethiopian cattle represents an untapped opportunity for adaptation to changing environments and for implementation of within-breed genetic improvement schemes. Our results provide a basis for future applications of genome-wide SNP data to exploit the unique genetic makeup of indigenous cattle breeds and to facilitate their improvement and conservation.     

海南尖峰岭中华穿山甲的分布与保护现状

中华穿山甲(Manis pentadactyla)是国家I级重点保护野生动物, 被IUCN红色名录列为极危(CR)物种, 也被列入CITES附录I。分布现状信息的匮乏是制约该物种保护规划制定与保护行动开展的关键因素。本文利用红外相机陷阱法和样线调查法, 于2020年8月至2021年11月在海南尖峰岭林区进行中华穿山甲海南亚种(M. p. pusilla)的监测调查, 分析评估了其分布与保护现状。调查发现: (1)尖峰岭林区7个公里网格内的8台红外相机拍摄到10次中华穿山甲出没的影像, 且在11个网格内见到35个觅食洞穴, 其中在红外相机损失率较低的南中区域记录到的中华穿山甲实体数量和洞穴数最多; (2)中华穿山甲主要分布于尖峰岭林区400-1,000 m海拔区域, 林区的国家公园一般控制区内仍然有3台红外相机拍摄到中华穿山甲及调查记录到18个觅食洞穴。结果表明: 海南尖峰岭林区仍然存在中华穿山甲种群, 目前可能被海拔等因素隔离为3个小种群; 人为干扰是影响该物种种群恢复的重要因素之一。因此, 本文提出如下建议: (1)管理部门要落实各项管护制度并加强巡护管理力度以把人为干扰程度降到最低; (2)管理部门在今后国家公园总体规划调整时可将有中华穿山甲等极危物种分布的一般控制区调整为核心保护区, 而在未调整规划前则要重点加强该区域的管理, 对建设项目布局等要尽量避开该区域; (3)尽快开展尖峰岭林区中华穿山甲的生境适宜性分析及廊道研究, 以使该物种得到更有效的保护和管理; (4)今后要定期开展其监测与保护研究, 明确尖峰岭林区中华穿山甲种群数量的动态变化; (5)在海南岛范围内开展中华穿山甲资源调查, 明确该物种在海南的种群数量及分布等情况。     

Comparative study of gut microbiota from captive and confiscated-rescued wild pangolins

Pangolins are among the most critically endangered animals due to widespread poaching and worldwide trafficking. Captive breeding is considered to be one way to protect them and increase the sizes of their populations. However, comparative studies of captive and wild pangolins in the context of gut microbiota are rare. Here, the gut microbiome of captive and confiscated-rescued wild pangolins is compared, and the effects of different periods of captivity and captivity with and without antibiotic treatment are considered. We show that different diets and periods of captivity, as well as the application of antibiotic therapy, can alter gut community composition and abundance in pangolins. Compared to wild pangolins, captive pangolins have an increased capacity for chitin and cellulose/hemicellulose degradation, fatty acid metabolism, and short-chain fatty acid synthesis, but a reduced ability to metabolize exogenous substances. In addition to increasing the ability of the gut microbiota to metabolize nutrients in captivity, captive breeding imposes some risks for survival by resulting in a greater abundance of antibiotic resistance genes and virulence factors in captive pangolins than in wild pangolins. Our study is important for the development of guidelines for pangolin conservation, including health assessment, disease prevention, and rehabilitation of wild pangolin populations.     

Historical processes and contemporary ocean currents drive genetic structure in the seagrass Thalassia hemprichii in the Indo‐Australian Archipelago

Understanding spatial patterns of gene flow and genetic structure is essential for the conservation of marine ecosystems. Contemporary ocean currents and historical isolation due to Pleistocene sea level fluctuations have been predicted to influence the genetic structure in marine populations. In the Indo‐Australian Archipelago (IAA), the world's hotspot of marine biodiversity, seagrasses are a vital component but population genetic information is very limited. Here, we reconstructed the phylogeography of the seagrass Thalassia hemprichii in the IAA based on single nucleotide polymorphisms (SNPs) and then characterized the genetic structure based on a panel of 16 microsatellite markers. We further examined the relative importance of historical isolation and contemporary ocean currents in driving the patterns of genetic structure. Results from SNPs revealed three population groups: eastern Indonesia, western Indonesia (Sunda Shelf) and Indian Ocean; while the microsatellites supported five population groups (eastern Indonesia, Sunda Shelf, Lesser Sunda, Western Australia and Indian Ocean). Both SNPs and microsatellites showed asymmetrical gene flow among population groups with a trend of southwestward migration from eastern Indonesia. Genetic diversity was generally higher in eastern Indonesia and decreased southwestward. The pattern of genetic structure and connectivity is attributed partly to the Pleistocene sea level fluctuations modified to a smaller level by contemporary ocean currents.     

Evolutionary history of the Paradoxurus palm civets – a new model for Asian biogeography

Aim Using molecular data and dental features, we investigated the genetic and morphological diversity among species of palm civets in the genus Paradoxurus, with a focus on the common palm civet, Paradoxurus hermaphroditus (Carnivora, Viverridae), in order to address biogeographic scenarios and provide recommendations for a taxonomic revision. Location Asia: Pakistan to the Lesser Sunda Islands. Methods We investigated the genetic diversity within Paradoxurus using two mitochondrial (cytochrome b, control region) and one nuclear (intron 7 of the β‐fibrinogen) markers. We used samples from 85 individuals of P. hermaphroditus (including 20 museum specimens) and one representative of each of the other species in the genus Paradoxurus: Paradoxurus jerdoni and Paradoxurus zeylonensis. DNA sequences were analysed using phylogenetic and haplotype network methods, and divergence dates were estimated for the clades retrieved. Furthermore, we examined dental characters from a large series of specimens and compared the morphological variation with the molecular data. Results Our phylogenetic analyses revealed that P. hermaphroditus is paraphyletic. We identified three major lineages distributed: (1) in the Indian subcontinent, south China, Hainan and in areas above 200 m in Indochina; (2) in Peninsular Malaysia, Java, Sumatra and in areas below 200 m in Indochina; and (3) in Borneo, the Philippines and the Mentawai archipelago. Our morphological observations were congruent with these three molecular lineages. Divergence date estimates inferred a Pliocene origin for Paradoxurus (2.8–5.7 Ma), with the three main clades diversifying from the mid–Early Pliocene to the end of the Pliocene. We suggest that the flooding of the Isthmus of Kra during the Pliocene was a major event shaping the diversification of Paradoxurus palm civets. We also hypothesize that the elevational segregation of the two lineages on the mainland could have resulted from the vegetational changes that were induced by Late Pliocene glacial episodes. Main conclusions The Isthmus of Kra is a major boundary between two major lineages of P. hermaphroditus. There is a need for a taxonomic revision for P. hermaphroditus, and we suggest that this species should be split into at least three species.     

Knowledge and Uses of African Pangolins as a Source of Traditional Medicine in Ghana

Traditional medicine has been practised in Ghana for centuries with the majority of Ghanaians still patronising the services of traditional healers. Throughout Africa a large number of people use pangolins as a source of traditional medicine, however, there is a dearth of information on the use of animals in folk medicine in Ghana, in particular the use of pangolins. The aim of this study was to determine the prevalent use of pangolins and the level of knowledge of pangolin use among traditional healers in Ghana for the treatment of human ailments. Data was gathered from 48 traditional healers using semi-structured interviews on the traditional medicinal use of pangolin body parts in the Kumasi metropolis of Ghana. The cultural importance index, relative frequency of citation, informant agreement ratio and use agreement values were calculated to ascertain the most culturally important pangolin body part as well as the level of knowledge dissemination among traditional healers with regards pangolin body parts. Our study revealed that 13 body parts of pangolins are used to treat various medicinal ailments. Pangolin scales and bones were the most prevalent prescribed body parts and indicated the highest cultural significance among traditional healing practices primarily for the treatment of spiritual protection, rheumatism, financial rituals and convulsions. Despite being classified under Schedule 1 of Ghana’s Wildlife Conservation Act of 1971 (LI 685), that prohibits anyone from hunting or being in possession of a pangolin, our results indicated that the use of pangolins for traditional medicinal purposes is widespread among traditional healers in Ghana. A study on the population status and ecology of the three species of African pangolins occurring in Ghana is urgently required in order to determine the impact this harvest for traditional medical purposes has on their respective populations as current levels appear to be unmonitored and unsustainable.     

Single‐nucleotide polymorphism discovery and panel characterization in the African forest elephant

The continuing decline in forest elephant (Loxodonta cyclotis) numbers due to poaching and habitat reduction is driving the search for new tools to inform management and conservation. For dense rainforest species, basic ecological data on populations and threats can be challenging and expensive to collect, impeding conservation action in the field. As such, genetic monitoring is being increasingly implemented to complement or replace more burdensome field techniques. Single‐nucleotide polymorphisms (SNPs) are particularly cost‐effective and informative markers that can be used for a range of practical applications, including population census, assessment of human impact on social and genetic structure, and investigation of the illegal wildlife trade. SNP resources for elephants are scarce, but next‐generation sequencing provides the opportunity for rapid, inexpensive generation of SNP markers in nonmodel species. Here, we sourced forest elephant DNA from 23 samples collected from 10 locations within Gabon, Central Africa, and applied double‐digest restriction‐site‐associated DNA (ddRAD) sequencing to discover 31,851 tags containing SNPs that were reduced to a set of 1,365 high‐quality candidate SNP markers. A subset of 115 candidate SNPs was then selected for assay design and validation using 56 additional samples. Genotyping resulted in a high conversion rate (93%) and a low per allele error rate (0.07%). This study provides the first panel of 107 validated SNP markers for forest elephants. This resource presents great potential for new genetic tools to produce reliable data and underpin a step‐change in conservation policies for this elusive species.     

Conservation genomics of anadromous Atlantic salmon across its North American range: outlier loci identify the same patterns of population structure as neutral loci

Anadromous Atlantic salmon (Salmo salar) is a species of major conservation and management concern in North America, where population abundance has been declining over the past 30 years. Effective conservation actions require the delineation of conservation units to appropriately reflect the spatial scale of intraspecific variation and local adaptation. Towards this goal, we used the most comprehensive genetic and genomic database for Atlantic salmon to date, covering the entire North American range of the species. The database included microsatellite data from 9142 individuals from 149 sampling locations and data from a medium‐density SNP array providing genotypes for >3000 SNPs for 50 sampling locations. We used neutral and putatively selected loci to integrate adaptive information in the definition of conservation units. Bayesian clustering with the microsatellite data set and with neutral SNPs identified regional groupings largely consistent with previously published regional assessments. The use of outlier SNPs did not result in major differences in the regional groupings, suggesting that neutral markers can reflect the geographic scale of local adaptation despite not being under selection. We also performed assignment tests to compare power obtained from microsatellites, neutral SNPs and outlier SNPs. Using SNP data substantially improved power compared to microsatellites, and an assignment success of 97% to the population of origin and of 100% to the region of origin was achieved when all SNP loci were used. Using outlier SNPs only resulted in minor improvements to assignment success to the population of origin but improved regional assignment. We discuss the implications of these new genetic resources for the conservation and management of Atlantic salmon in North America.     

Construction, characterization and FISH mapping of a bacterial artificial chromosome library of Chinese pangolin (Manis pentadactyla)

Chinese pangolins as a representative species in the order Pholidota have highly specified morphological characters and occupy an important place in the mammalian phylogenetic tree. To obtain genomic data for this species, we have constructed a bacterial artificial chromosome (BAC) library of Chinese pangolin. The library contains 208,272 clones with an average insert size of 122.1 kb and represents approximately eight times the Chinese pangolin haploid genome (if we assume that the Chinese pangolins have a genome size similar to human). One hundred and twenty randomly-selected BAC clones were mapped onto Chinese pangolin chromosomes by fluorescence in situ hybridization (FISH), showing a largely unbiased chromosomal distribution. Several clones containing repetitive DNA and ribosomal DNA genes were also found. The BAC library and FISH mapped BAC clones are useful resources for comparative genomics and cytogenetics of mammals and in particular, the ongoing genome sequencing project of Chinese pangolins.     

The first authenticated record of the pangolin tick <Emphasis Type="Italic">Amblyomma javanense</Emphasis> (Acari: Ixodidae) in Singapore,with notes on its biology and conservation

The first authenticated record of the co-endangered pangolin tick (Amblyomma javanense) is reported from the Republic of Singapore, based on specimens of males and females collected from the Sunda pangolin (Manis javanica). The biology of the species is reviewed and steps for its conservation are outlined to ensure its continued survival across its range.     

Wallace's Line and plant distributions: two or three phytogeographical areas and where to group Java?

Wallace's Line or its variants divide the Malay Archipelago or Malesia into a western and eastern area, but is this suitable for plant distributions? Indeed, all boundaries satisfactorily divide Malesia into two parts, stopping far more species east or west of a line than disperse over the boundary. However, phenetic analyses (principal components analysis, nonmetric multidimensional scaling analysis and the unweighted pair group method with arithmetic mean) of 7340 species distributions revealed a stronger partitioning of Malesia into three instead of two regions: the western Sunda Shelf minus Java (Malay Peninsula, Sumatra, Borneo), central Wallacea (Philippines, Sulawesi, Lesser Sunda Islands, Moluccas, with Java), and the eastern Sahul Shelf (New Guinea). Java always appears to be part of Wallacea, probably because of its mainly dry monsoon climate. The three phytogeographic areas equal the present climatic division of Malesia. An everwet climate exists on the Sunda and Sahul Shelves, whereas most of Wallacea has a yearly dry monsoon. During glacial maxima, the Sunda and Sahul Shelves became land areas connected with Asia and Australia, respectively, whereas sea barriers remained within Wallacea. Consequently, the flora of the two shelves is more homogeneous than the Wallacean flora. Wallacea is a distinct area because it comprises many endemic, drought tolerant floristic elements. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 531–545.     

Neanderthal and Denisova genetic affinities with contemporary humans: Introgression versus common ancestral polymorphisms

Analyses of the genetic relationships among modern humans, Neanderthals and Denisovans have suggested that 1–4% of the non-Sub-Saharan African gene pool may be Neanderthal derived, while 6–8% of the Melanesian gene pool may be the product of admixture between the Denisovans and the direct ancestors of Melanesians. In the present study, we analyzed single nucleotide polymorphism (SNP) diversity among a worldwide collection of contemporary human populations with respect to the genetic constitution of these two archaic hominins and Pan troglodytes (chimpanzee). We partitioned SNPs into subsets, including those that are derived in both archaic lineages, those that are ancestral in both archaic lineages and those that are only derived in one archaic lineage. By doing this, we have conducted separate examinations of subsets of mutations with higher probabilities of divergent phylogenetic origins. While previous investigations have excluded SNPs from common ancestors in principal component analyses, we included common ancestral SNPs in our analyses to visualize the relative placement of the Neanderthal and Denisova among human populations. To assess the genetic similarities among the various hominin lineages, we performed genetic structure analyses to provide a comparison of genetic patterns found within contemporary human genomes that may have archaic or common ancestral roots. Our results indicate that 3.6% of the Neanderthal genome is shared with roughly 65.4% of the average European gene pool, which clinally diminishes with distance from Europe. Our results suggest that Neanderthal genetic associations with contemporary non-Sub-Saharan African populations, as well as the genetic affinities observed between Denisovans and Melanesians most likely result from the retention of ancient mutations in these populations.     

Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing

The KwaZulu‐Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu‐Natal Province, South Africa. In this study, we developed a single‐nucleotide polymorphism (SNP) dataset from double‐digest restriction site‐associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focusing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high‐quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well‐supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north–south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.