GENETIC EVIDENCE FOR A PLEISTOCENE POPULATION EXPLOSION

Phylogeographic Diversity Analysis of Bipolaris sorokiniana (Sacc.) Shoemaker Causing Spot Blotch Disease in Wheat and Barley

Bipolaris sorokiniana is a fungal pathogen that infects wheat, barley, and other crops, causing spot blotch disease. The disease is most common in humid, warm, wheat-growing regions, with South Asia's Eastern Gangetic Plains serving as a hotspot. There is very little information known about its genetic variability, demography, and divergence period. The current work is the first to study the phylogeographic patterns of B. sorokiniana isolates obtained from various wheat and barley-growing regions throughout the world, with the goal of elucidating the demographic history and estimating divergence times. In this study, 162 ITS sequences, 18 GAPDH sequences, and 74 TEF-1α sequences from B. sorokiniana obtained from the GenBank, including 21 ITS sequences produced in this study, were used to analyse the phylogeographic pattern of distribution and evolution of B. sorokiniana infecting wheat and barley. The degrees of differentiation among B. sorokiniana sequences from eighteen countries imply the presence of a broad and geographically undifferentiated global population. The study provided forty haplotypes. The H_1 haplotype was identified to be the ancestral haplotype, followed by H_29 and H_27, with H_1 occupying a central position in the median-joining network and being shared by several populations from different continents. The phylogeographic patterns of species based on multi-gene analysis, as well as the predominance of a single haplotype, suggested that human-mediated dispersal may have played a significant role in shaping this pathogen's population. According to divergence time analysis, haplogroups began at the Plio/Pleistocene boundary.

Population genetic structure and evolutionary genetics of Anopheles sinensis based on knockdown resistance (kdr) mutations and mtDNA-COII gene in China-Laos, Thailand-Laos, and Cambodia-Laos borders

BACKGROUND

Vector control is still a pivotal method for preventing malaria, and its potency is weakened by the increasing resistance of vectors to chemical insecticides. As the most abundant and vital malaria vector in Southeast Asia, the chemical insecticide resistance status in Anopheles sinensis remains elusive in Laos, which makes it imperative to evaluate the true nature of chemical insecticide resistance-associated genetic mutations in An. sinensis in Laos.

METHODS

Adult An. sinensis were collected from three border regions in Laos. DNA was extracted from individual mosquitoes. PCR amplification and DNA sequencing of a fragment containing codon 1014 of the voltage-gated sodium channel (vgsc) gene were completed to study the kdr allele frequency distribution, kdr intron polymorphism, population genetic diversity, and the evolutionary status of the kdr codon. The mitochondrial cytochrome c oxidase subunit II gene (COII) was amplified and sequenced to examine population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

Nine wild kdr haplotypes of the vgsc gene were detected in this study, and eight of them, namely 1014L1, 1014L2, 1014L4, 1014L7, 1014L9, 1014L10, 1014L11, and 1014L21, were discovered in the China-Laos border (northern Laos), while 1014L3 was only detected in the Thailand-Laos border (northwestern Laos) and Cambodia-Laos border (southern Laos). The newly identified haplotype, 1014L21, was uniquely distributed in the China-Laos border and was not identified in other countries. Based on sequence analysis of the mitochondrial COII genes, significant genetic differentiation and limited gene flow were detected between the China-Laos and Cambodia-Laos An. sinensis populations, which suggested that those two regions were genetically isolated. The distinct distribution of the kdr haplotype frequencies is probably the result of geographical isolation in mosquito populations.

CONCLUSIONS

Lack of kdr mutations in the vgsc gene was probably due to genetic isolation and the absence of intense selection pressure in the three border regions of Laos. This study reveals that pyrethroid-based chemical insecticides are still appropriate for battling An. sinensis in parts of Laos, and routine monitoring of chemical insecticide resistance should be continuously implemented and focused on more restricted areas as part of chemical insecticide resistance management.

Spatial Variation of Acanthophlebia cruentata (Ephemeroptera), a Mayfly Endemic to Te Ika-a-Māui—North Island of Aotearoa, New Zealand

Simple Summary

Patterns of variation within a widespread species can provide evidence of population history. Adjacent stable populations with gene flow show clinal genetic divergence forming a pattern of isolation by distance. Populations that have grown due to an increase in potential habitat such as forest cover since the last glacial maximum will have low genetic variation showing patterns of range expansion. The mayfly Acanthophlebia cruentata of Aotearoa, New Zealand, is widespread in streams of North Island but absent from the cooler South Island. Mayfly nymphs are restricted to streams but adults fly, facilitating gene flow among catchments. We detected higher genetic diversity at lower latitudes of this mayfly’s range compared to most of its distribution, concordant with predictions of limited forest cover in New Zealand during Pleistocene glacial periods. A signature of recent range expansion was observed in the higher latitudes. Despite initial observation suggesting mayfly size correlated with latitude, we found sex, elevation and sampling date were significant predictors of size, and some size variation is also explained by three regional groups based on haplotype distribution.

Abstract

The mayfly Acanthophlebia cruentata of Aotearoa, New Zealand, is widespread in Te Ika-a-Māui North Island streams, but has never been collected from South Island despite land connection during the last glacial maximum. Population structure of this mayfly might reflect re-colonisation after volcanic eruptions in North Island c1800 years ago, climate cycling or conceal older, cryptic diversity. We collected population samples from 33 locations to estimate levels of population genetic diversity and to document phenotypic variation. Relatively low intraspecific haplotype divergence was recorded among mitochondrial cytb sequences from 492 individuals, but these resolved three geographic-haplotype regions (north, west, east). We detected a signature of isolation by distance at low latitudes (north) but evidence of recent population growth in the west and east. We did not detect an effect of volcanic eruptions but infer range expansion into higher latitudes from a common ancestor during the last glacial period. As judged from wing length, both sexes of adult mayflies were larger at higher elevation and we found that haplotype region was also a significant predictor of Acanthophlebia cruentata size. This suggests that our mitochondrial marker is concordant with nuclear genetic differences that might be explained by founder effect during range expansion.

Insights on the historical biogeography of Philippine domestic pigs and its relationship with continental domestic pigs and wild boars

The Philippine archipelago was believed to have never been connected to the Asian continent even during the severe Quaternary sea-level drops. As a result, the history of domestic pig (Sus scrofa) dispersal in the Philippines remains controversial and must have some anthropogenic origin associated with human migration events. In this study, the context of origin, dispersal, and the level of genetic introgression in Philippine domestic pigs were deduced using mitochondrial DNA D-loop analysis altogether with domestic pigs and wild boar corresponding to their geographic origin. The results revealed considerable genetic diversity (0.900±0.016) and widespread Asian pig-ancestry (94.60%) in the phylogenetic analysis, with admixed European pig-origin (5.10%) harboring various fractions of ancestry from Berkshire and Landrace. The close genetic connection between the continental wild boars and domestic pigs present in the Philippine domestic pigs corroborates our hypothesis of a genetic signal that may be associated with the recently reported multiple waves of human migrations to the Philippines. The Haplogroup D7, reported to occur only in Indo-Burma Biodiversity Hotspots, included a high frequency of Philippine domestic pig haplotypes (54.08%), which poses an interesting challenge because its distribution is not consistent with the hypothesized migration route of Neolithic Austronesian-speaking populations. We detected the first Pacific Clade signature and ubiquitously distributed D2 haplotypes (Asian major) on several Philippine islands. The analyses of mismatch distribution and neutrality test were consistent with the Bayesian skyline plot which showed a long stationary period of effective population size. The population decline was consistent with the pronounced population bottleneck in Asian and European pigs during the interglacial periods of the Pleistocene. The results of this study will support the conservation strategies and improvements of economically important genetic resources in the Philippines.

Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border

BACKGROUND

Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area.

METHODS

1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F_ST = 0.97824, 0.97343, P < 0.05), as well as lower gene flow (Nm = 0.01112, 0.01365) in the An. nivipes populations. Phylogenetic analysis and SAMOVA revealed a gene barrier restricting gene flow among three An. nivipes populations. Mantel test suggested a significant correlation between geography and gene distance in all An. nivipes populations (Z = 44,983.1865, r = 0.5575, P = 0.0070). Neutrality test and Mismatch distribution revealed a recent population expansion of An. nivipes in the Cambodia-Laos population.

CONCLUSIONS

Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. Based on sequences analysis of COII, a genetic barrier between Cambodia-Laos and two Indian populations was found, and a recent population expanding or selecting of An. nivipes occurred in the Cambodia-Laos population, suggesting that COII might be an effective marker for describing the molecular phylogeny of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.

Molecular phylogeography and species distribution modelling evidence of 'oceanic' adaptation for Actinidia eriantha with a refugium along the oceanic-continental gradient in a biodiversity hotspot

BACKGROUND

Refugia is considered to be critical for maintaining biodiversity; while discerning the type and pattern of refugia is pivotal for our understanding of evolutionary processes in the context of conservation. Interglacial and glacial refugia have been studied throughout subtropical China. However, studies on refugia along the oceanic-continental gradient have largely been ignored. We used a liana Actinidia eriantha, which occurs across the eastern moist evergreen broad-leaved forests of subtropical China, as a case study to test hypotheses of refugia along the oceanic-continental gradient and 'oceanic' adaptation.

RESULTS

The phylogeographic pattern of A. eriantha was explored using a combination of three cpDNA markers and 38 nuclear microsatellite loci, Species distribution modelling and dispersal corridors analysis. Our data showed intermediate levels of genetic diversity [haplotype diversity (h_T) = 0.498; unbiased expected heterozygosity (UH_E) = 0.510] both at the species and population level. Microsatellite loci revealed five clusters largely corresponding to geographic regions. Coalescent time of cpDNA lineages was dated to the middle Pliocene (ca. 4.03 Ma). Both geographic distance and climate difference have important roles for intraspecific divergence of the species. The Zhejiang-Fujian Hilly Region was demonstrated to be a refugium along the oceanic-continental gradient of the species and fit the 'refugia in refugia' pattern. Species distribution modelling analysis indicated that Precipitation of Coldest Quarter (importance of 44%), Temperature Seasonality (29%) and Mean Temperature of Wettest Quarter (25%) contributed the most to model development. By checking the isolines in the three climate layers, we found that A. eriantha prefer higher precipitation during the coldest quarter, lower seasonal temperature difference and lower mean temperature during the wettest quarter, which correspond to 'oceanic' adaptation. Actinidia eriantha expanded to its western distribution range along the dispersal corridor repeatedly during the glacial periods.

CONCLUSIONS

Overall, our results provide integrated evidence demonstrating that the Zhejiang-Fujian Hilly Region is a refugium along the oceanic-continental gradient of Actinidia eriantha in subtropical China and that speciation is attributed to 'oceanic' adaptation. This study gives a deeper understanding of the refugia in subtropical China and will contribute to the conservation and utilization of kiwifruit wild resources in the context of climate change.

Comparative phylogeography study reveals introgression and incomplete lineage sorting during rapid diversification of Rhodiola

BACKGROUND AND AIMS

Many plant taxa in the Qinghai-Tibetan Plateau (QTP) and the Hengduan Mountains (HM) radiated rapidly during the Quaternary but with frequent secondary contact between diverging populations. Incomplete lineage sorting and introgressive hybridization might be involved during the rapid radiation, but their effects on phylogeography have not been fully determined.

METHODS

We investigated the chloroplast DNA (cpDNA)/internal transcribed spacer (ITS) sequence variations of 611 samples of Rhodiola bupleuroides, R. discolor, R. fastigiata and R. chrysanthemifolia from the QTP and HM to compare the phylogeographic patterns between the four species with different evolutionary histories, geographic ranges and reproductive modes.

KEY RESULTS

The divergence times of these species were consistent with the last peak of in situ speciation in the HM. While closely related species exhibited different phylogeographic patterns, they shared several ribotypes and haplotypes in sympatric populations, suggesting introgressive hybridization. A significant phylogenetic discordance between ribotypes and haplotypes was detected in three species, implying incomplete lineage sorting. Rhodiola discolor houses an extraordinary richness of cpDNA haplotypes, and this finding may be attributed to adaptive radiation.

CONCLUSION

In addition to geographic isolation and climate oscillations during the Quaternary, both introgressive hybridization and incomplete lineage sorting play important roles in species that experienced rapid diversification in the QTP and HM.

Worldwide diversity, association potential, and natural selection in the superimposed taste genes, CD36 and GNAT3

CD36 and GNAT3 mediate taste responses, with CD36 acting as a lipid detector and GNAT3 acting as the α subunit of gustducin, a G protein governing sweet, savory, and bitter transduction. Strikingly, the genes encoding CD36 and GNAT3 are genomically superimposed, with CD36 completely encompassing GNAT3. To characterize genetic variation across the CD36-GNAT3 region, its implications for phenotypic diversity, and its recent evolution, we analyzed from ~2,500 worldwide subjects sequenced by the 1000 Genomes Project (1000GP). CD36-GNAT3 harbored extensive diversity including 8,688 single-nucleotide polymorphisms (SNPs), 414 indels, and other complex variants. Sliding window analyses revealed that nucleotide diversity and population differentiation across CD36-GNAT3 were consistent with genome-wide trends in the 1000GP (π = 0.10%, P = 0.64; FST = 9.0%, P = 0.57). In addition, functional predictions using SIFT and PolyPhen-2 identified 60 variants likely to alter protein function, and they were in weak linkage disequilibrium (r2 < 0.17), suggesting their effects are largely independent. However, the frequencies of predicted functional variants were low (P¯ = 0.0013), indicating their contributions to phenotypic variance on population scales are limited. Tests using Tajima's D statistic revealed that pressures from natural selection have been relaxed across most of CD36-GNAT3 during its recent history (0.39 < P < 0.67). However, CD36 exons showed signs of local adaptation consistent with prior reports (P < 0.035). Thus, CD36 and GNAT3 harbor numerous variants predicted to affect taste sensitivity, but most are rare and phenotypic variance on a population level is likely mediated by a small number of sites.

Plio-Pleistocene climatic change drives allopatric speciation and population divergence within the Scrophularia incisa complex (Scrophulariaceae) of desert and steppe subshrubs in Northwest China

Numerous temperate plants and animals on the Qinghai-Tibet Plateau (QTP) are hypothesized to have differentiated due to vicariant allopatric speciation associated with the geologic uplifts. However, this hypothesis has rarely been tested through a phylogeographic study of relative species in a broader geographic context, including the QTP, Tianshan Mountains, Mongolian Plateau, and surrounding regions. To understand the speciation and diversification process of plants across this wide area, phylogeographic analysis were examined from Scrophularia incisa and two other closely relative species comprising S. kiriloviana and S. dentata. Thirty-two populations of the three close relatives were genotyped using chloroplast DNA fragments and nuclear microsatellite loci to assess population structure and diversity, supplemented by phylogenetic dating, ancestral area reconstructions and species distribution modelings, as well as niche identity tests. Our chloroplast DNA (cpDNA) phylogeny showed that this monophyletic group of desert and steppe semi-shrub is derived from a Middle Pliocene ancestor of the Central Asia. Lineages in Central Asia vs. China diverged through climate/tectonic-induced vicariance during Middle Pliocene. Genetic and ENM data in conjunction with niche differentiation analyses support that the divergence of S. incisa, S. dentata and S. kiriloviana in China lineage proceeded through allopatric speciation, might triggered by early Pleistocene climate change of increase of aridification and enlargement of deserts, while subsequent climate-induced cycles of range contractions/expansions enhanced the geographical isolation and habit fragmentation of these taxa. These findings highlight the importance of the Plio-Pleistocene climate change in shaping genetic diversity and driving speciation in temperate steppes and deserts of Northwestern China.

Genetic diversification of the Kanehira bitterling Acheilognathus rhombeus inferred from mitochondrial DNA, with comments on the phylogenetic relationship with its sister species Acheilognathus barbatulus

The Kanehira bitterling, Acheilognathus rhombeus, is a freshwater fish, discontinuously distributed in western Japan and the Korean Peninsula. Unusually among bitterling it is an autumn-spawning species and shows developmental diapause. Consequently, the characterization of its evolutionary history is significant not only in the context of the fish assemblage of East Asia, but also for understanding life-history evolution. This study aimed to investigate the phylogeography of A. rhombeus and its sister species Acheilognathus barbatulus, distributed in China, using a mitochondrial analysis of the ND1 gene from 311 samples collected from 50 localities in Japan and continental Asia. Phylogenetic analysis revealed that A. barbatulus is included in A. rhombeus and genetically closer to Japanese A. rhombeus than to Korean A. rhombeus. Divergence of Korean A. rhombeus and A. barbatulus from Japanese A. rhombeus was estimated to be from the late Pliocene (3.44 Mya) and the early Pleistocene (1.98 Mya), respectively. Each event closely coincided with the time of the Japan Sea opening. Japanese A. rhombeus comprised seven lineages: three in Honshu and four in Kyushu. One lineage in central Kyushu was genetically closer to the Honshu lineages than to other lineages in northern Kyushu. Divergence of Japanese lineages was estimated to be from the early to middle Pleistocene (0.55-0.93 Mya), during a period of geological and paleoclimatic change, including volcanic activity. Population expansion in the late Pleistocene (<0.10 Ma) was suggested in many of the lineages, which accords with other freshwater fishes. Biogeographically the ancestral A. rhombeus/A. barbatulus was likely to have repeatedly colonized Japan from the continent through land bridges in the late Pliocene and the early Pleistocene. However, the close genetic relationship between Japanese A. rhombeus and A. barbatulus suggests another possibility, with the second colonization occurring in reverse, from Japan to China. The small genetic distance between them indicates that the colonization occurred later than colonization events of other freshwater fishes, including other bitterling species.

Assessing signals of selection and historical demography to develop conservation strategies in the Chilean emblematic Araucaria araucana

Loss of genetic diversity reduces the ability of species to evolve and respond to environmental change. Araucaria araucana is an emblematic conifer species from southern South America, with important ethnic value for the Mapuche people (Pehuenche); the Chilean Government has catalogued its conservation status as vulnerable. Climatic fluctuations were potentially a major impact in the genetic variation within many tree species. In this context, the restricted geographic distribution of A. araucana in Chile appears to be a consequence of the Last Glacial Maximum (LGM). During the past two centuries, strong human intervention has also affected the geographical distribution and population sizes of A. araucana. Reduction of population size may cause loss of genetic diversity, which could affect frequency of adaptive loci. The aims of this study were to know the existence of potential loci under selection and populations with genetic, demographic disequilibrium in the Chilean distribution of A. araucana. Based on 268 polymorphic AFLP loci, we have investigated potential loci under selection and genetic, demographic disequilibrium within seven Chilean populations of Araucaria araucana. Correlation of 41 outlier loci with the environmental variables of precipitation and temperature reveals signatures of selection, whereas 227 neutral loci provide estimates of demographic equilibrium and genetic population structure. Three populations are recommended as priorities for conservation.

Molecular phylogeny of the Anopheles hyrcanus group (Diptera: Culicidae) based on rDNA-ITS2 and mtDNA-COII

BACKGROUND

The Anopheles hyrcanus group, which includes 25 species, is widely distributed in the Oriental and Palaearctic regions. Given the difficulty in identifying cryptic or sibling species based on their morphological characteristics, molecular identification is regarded as an important complementary approach to traditional morphological taxonomy. The aim of this study was to reconstruct the phylogeny of the Hyrcanus group using DNA barcoding markers in order to determine the phylogenetic correlations of closely related taxa and to compare these markers in terms of identification efficiency and genetic divergence among species.

METHODS

Based on data extracted from the GenBank database and data from the present study, we used 399 rDNA-ITS2 sequences of 19 species and 392 mtDNA-COII sequences of 14 species to reconstruct the molecular phylogeny of the Hyrcanus group across its worldwide range. We also compared the performance of rDNA-ITS2 against that of mtDNA-COII to assess the genetic divergence of closely related species within the Hyrcanus group.

RESULTS

Average interspecific divergence for the rDNA-ITS2 sequence (0.376) was 125-fold higher than the average intraspecies divergence (0.003), and average interspecific divergence for the mtDNA-COII sequence (0.055) was eightfold higher than the average intraspecies divergence (0.007). The barcoding gap ranged from 0.015 to 0.073 for rDNA-ITS2, and from 0.017 to 0.025 for mtDNA-COII. Two sets of closely related species, namely, Anophels lesteri and An. paraliae, and An. sinensis, An. belenrae and An. kleini, were resolved by rDNA-ITS2. In contrast, the relationship of An. sinensis/An. belenrae/An. kleini was poorly defined in the COII tree. The neutrality test and mismatch distribution revealed that An. peditaeniatus, An. hyrcanus, An. sinensis and An. lesteri were likely to undergo hitchhiking or population expansion in accordance with both markers. In addition, the population of an important vivax malaria vector, An. sinensis, has experienced an expansion after a bottleneck in northern and southern Laos.

CONCLUSIONS

The topology of the Hyrcanus group rDNA-ITS2 and mtDNA-COII trees conformed to the morphology-based taxonomy for species classification rather than for that for subgroup division. rDNA-ITS2 is considered to be a more reliable diagnostic tool than mtDNA-COII in terms of investigating the phylogenetic correlation between closely related mosquito species in the Hyrcanus group. Moreover, the population expansion of an important vivax malaria vector, An. sinensis, has underlined a potential risk of malaria transmission in northern and southern Laos. This study contributes to the molecular identification of the Anopheles hyrcanus group in vector surveillance.

Genome-wide assessment elucidates connectivity and the evolutionary history of the highly dispersive marine invertebrate Littoraria flava (Littorinidae: Gastropoda)

An important goal of marine population genetics is to understand how spatial connectivity patterns are influenced by historical and evolutionary factors. In this study, we evaluate the demographic history and population structure of Littoraria flava, a highly dispersive marine gastropod in the Brazilian intertidal zone. To test the hypotheses that the species has (1) historically high levels of gene flow on a macrogeographical spatial scale and (2) a distribution in rocky shores that consists of subpopulations, we collected specimens along the Brazilian coastline and combined different sets of genetic markers (mitochondrial DNA, ITS-2 and single nucleotide polymorphisms) with niche-based modelling to predict its palaeodistribution. Low genetic structure was observed, as well as high gene flow over long distances. The demographic analyses suggest that L. flava has had periods of population bottlenecks followed by expansion. According to both palaeodistribution and coalescent simulations, these expansion events occurred during the Pleistocene interglacial cycles (21 kya) and the associated climatic changes were the probable drivers of the distribution of the species. This is the first phylogeographical study of a marine gastropod on the South American coast based on genomic markers associated with niche modelling.

Population Genetic Structure and Contribution of Philippine Chickens to the Pacific Chicken Diversity Inferred From Mitochondrial DNA

The Philippines is considered one of the biodiversity hotspots for animal genetic resources. In spite of this, population genetic structure, genetic diversity, and past population history of Philippine chickens are not well studied. In this study, phylogeny reconstruction and estimation of population genetic structure were based on 107 newly generated mitochondrial DNA (mtDNA) complete D-loop sequences and 37 previously published sequences of Philippine chickens, consisting of 34 haplotypes. Philippine chickens showed high haplotypic diversity (Hd = 0.915 ± 0.011) across Southeast Asia and Oceania. The phylogenetic analysis and median-joining (MJ) network revealed predominant maternal lineage haplogroup D classified throughout the population, while support for Philippine-Pacific subclade was evident, suggesting a Philippine origin of Pacific chickens. Here, we observed Philippine red junglefowls (RJFs) at the basal position of the tree within haplogroup D indicating an earlier introduction into the Philippines potentially via mainland Southeast Asia (MSEA). Another observation was the significantly low genetic differentiation and high rate of gene flow of Philippine chickens into Pacific chicken population. The negative Tajima's D and Fu's Fs neutrality tests revealed that Philippine chickens exhibited an expansion signal. The analyses of mismatch distribution and neutrality tests were consistent with the presence of weak phylogeographic structuring and evident population growth of Philippine chickens (haplogroup D) in the islands of Southeast Asia (ISEA). Furthermore, the Bayesian skyline plot (BSP) analysis showed an increase in the effective population size of Philippine chickens, relating with human settlement, and expansion events. The high level of genetic variability of Philippine chickens demonstrates conservation significance, thus, must be explored in the future.

Population Genetics of the Black Citrus Aphid Aphis aurantii (Hemiptera, Aphididae) in China

The black citrus aphid, Aphis aurantiiBoyer de Fonscolombe, 1841, is one of the most destructive pests in commercial tea plantations and gardens in China. In this study, we investigated the population genetic structure of A. aurantii based on the concatenated sequences of two mitochondrial genes, cytochrome c oxidase I (cox1) and cytochrome b (cytb). A total of 166 haplotypes were identified from 177 individuals collected at 11 locations in China. The whole Chinese A. aurantii population showed a low nucleotide diversity (0.00968) and a high population diversity (haplotype diversity; 0.9991). The haplotypes of the 11 local populations were widely distributed in the neighbor-joining phylogenetic tree and haplotype network diagram, whereas no apparent lineages were detected. Gene flow analysis showed gene exchanges among local populations. The pairwise Fst values revealed a certain amount of genetic difference among local populations. Analysis of molecular variance (AMOVA) reflected genetic differences both within and among populations. The isolation by distance (IBD) analysis revealed a high positive correlation between the geographic distance and genetic distance of the different populations. Neutral test and mismatch distribution suggested that A. aurantii may have experienced recent population expansion events.

Cytotype distribution and chloroplast phylogeography of the Actinidia chinensis complex

BACKGROUND

Plant phylogeographic studies of species in subtropical China have mainly focused on rare and endangered species, whereas few studies have been conducted on taxa with relatively wide distribution, especially polyploid species. We investigated the cytotype and haplotype distribution pattern of the Actinidia chinensis complex, a widespread geographically woody liana with variable ploidy in subtropical China comprising two varieties, with three chloroplast fragments DNA (ndhF-rpl132, rps16-trnQ and trnE-trnT). Macroevolutionary, microevolutionary and niche modeling tools were also combined to disentangle the origin and the demographic history of the species or cytotypes.

RESULTS

The ploidy levels of 3338 individuals from 128 populations sampled throughout the species distribution range were estimated with flow cytometry. The widespread cytotypes were diploids followed by tetraploids and hexaploids, whereas triploids and octoploids occurred in a few populations. Thirty-one chloroplast haplotypes were detected. The genetic diversity and genetic structure were found to be high between varieties (or ploidy races) chinensis and deliciosa. Our results revealed that these two varieties inhabit significantly different climatic niche spaces. Ecological niche models (ENMs) indicate that all varieties' ranges contracted during the Last Inter Glacial (LIG), and expanded eastward or northward during the Last Glacial Maximum (LGM).

CONCLUSIONS

Pliocene and Plio-Pleistocene climatic fluctuations and vicariance appear to have played key roles in shaping current population structure and historical demography in the A. chinensis complex. The polyploidization process also appears to have played an important role in the historical demography of the complex through improving their adaptability to environmental changes.

High-Level Gene Flow Restricts Genetic Differentiation in Dairy Cattle Populations in Thailand: Insights from Large-Scale Mt D-Loop Sequencing

Domestication and artificial selection lead to the development of genetically divergent cattle breeds or hybrids that exhibit specific patterns of genetic diversity and population structure. Recently developed mitochondrial markers have allowed investigation of cattle diversity worldwide; however, an extensive study on the population-level genetic diversity and demography of dairy cattle in Thailand is still needed. Mitochondrial D-loop sequences were obtained from 179 individuals (hybrids of Bos taurus and B. indicus) sampled from nine different provinces. Fifty-one haplotypes, of which most were classified in haplogroup "I", were found across all nine populations. All sampled populations showed severely reduced degrees of genetic differentiation, and low nucleotide diversity was observed in populations from central Thailand. Populations that originated from adjacent geographical areas tended to show high gene flow, as revealed by patterns of weak network structuring. Mismatch distribution analysis was suggestive of a stable population, with the recent occurrence of a slight expansion event. The results provide insights into the origins and the genetic relationships among local Thai cattle breeds and will be useful for guiding management of cattle breeding in Thailand.

Studying genetic population structure to shed light on the demographic explosion of the rare species Barbitistes vicetinus (Orthoptera, Tettigoniidae)

Insect outbreaks usually involve important ecological and economic consequences for agriculture and forestry. The short-winged bush-cricket Barbitistes vicetinus Galvagni & Fontana, 1993 is a recently described species that was considered rare until ten years ago, when unexpected population outbreaks causing severe defoliations across forests and crops were observed in north-eastern Italy. A genetic approach was used to analyse the origin of outbreak populations. The analysis of two mitochondrial regions (Cytochrome Oxidase I and II and 12S rRNA-Control Region) of 130 samples from the two disjunct ranges (Euganean and Berici Hills) showed high values of haplotype diversity and revealed a high geographical structure among populations of the two ranges. The high genetic variability observed supports the native origin of this species. In addition, results suggest that unexpected outbreaks are not a consequence of a single or few pestiferous haplotypes but rather the source of outbreaks are local populations which have experienced an increase in each area. The recent outbreaks have probably appeared independently of the genetic haplotypes whereas environmental conditions could have affected the outbreak populations. These findings contribute to a growing understanding of the status and evolutionary history of the pest that would be useful for developing and implementing biological control strategies for example by maximizing efforts to locate native natural enemies.

Genetic structure and population demographics of Hypnoidus bicolor (Coleoptera: Elateridae) in the Canadian Prairies

BACKGROUND

Following banning of the pesticide lindane in most counties, wireworms (i.e., the soil-living larval stages of click beetles) have become major pests of a variety of economically important field crops. Hypnoidus bicolor is a common pest species in the Canadian Prairies. However, little is known about its life history, which impedes the development of effective integrated pest management (IPM) strategies. Population genetic approaches have the potential to assist in the development of IPM.

RESULTS

We sequenced a 622-bp fragment of the COX1 gene from 326 H. bicolor wireworm and click beetles collected from 13 localities on the Canadian Prairies. Two genetically distinct (>4.66% sequence divergence) clades were identified, suggesting that they may be part of a species complex. Clade A predominated and increased in prevalence the further east samples were collected, whereas the opposite was true for clade B. Clade B appears to be comprised of two mitochondrial DNA groups, however, one group was represented by only one haplotype. Both clades were characterized by uneven gene flow among populations with low levels of regional genetic structuring. Clade A appeared to have undergone population and range expansions, which may coincide with the advent of intensive agriculture practices in the prairies.

CONCLUSION

Knowledge of species composition and population structure is important for the development of effective IPM strategies but is often lacking for wireworms. Our study fills these knowledge gaps for a predominant pest species in the prairies, H. bicolor, by providing robust evidence for cryptic forms and characterizing its dispersal patterns and population dynamics. © 2021 Society of Chemical Industry.

Late Pleistocene Expansion of Small Murid Rodents across the Palearctic in Relation to the Past Environmental Changes

We investigated the evolutionary history of the striped field mouse to identify factors that initiated its past demographic changes and to shed light on the causes of its current genetic structure and trans-Eurasian distribution. We sequenced mitochondrial cyt b from 184 individuals, obtained from 35 sites in central Europe and eastern Mongolia. We compared genetic analyses with previously published historical distribution models and data on environmental and climatic changes. The past demographic changes displayed similar population trends in the case of recently expanded clades C1 and C3, with the glacial (MIS 3–4) expansion and postglacial bottleneck preceding the recent expansion initiated in the late Holocene and were related to environmental changes during the upper Pleistocene and Holocene. The past demographic trends of the eastern Asian clade C3 were correlated with changes in sea level and the formation of new land bridges formed by the exposed sea shelf during the glaciations. These data were supported by reconstructed historical distribution models. The results of our genetic analyses, supported by the reconstruction of the historical spatial distributions of the distinct clades, confirm that over time the local populations mixed as a consequence of environmental and climatic changes resulting from cyclical glaciation and the interglacial period during the Pleistocene.

Demography reveals populational expansion of a recently extinct Iberian ungulate

Reconstructing the demographic history of endangered taxa is paramount to predict future fluctuations and disentangle the contributing factors. Extinct taxa or populations might also provide key insights in this respect by means of the DNA extracted from museum specimens. Nevertheless, the degraded status of biological material and the limited number of records may pose some constraints. For this reason, identifying all available sources, including private and public biological collections, is a crucial step forward. In this study, we reconstructed the demographic history based on cytochrome-b sequence data of the Pyrenean ibex (Capra pyrenaica pyrenaica), a charismatic taxon of the European wildlife that became extinct in the year 2000. Moreover, we built a database of the museum specimens available in public biological collections worldwide and genotyped a privately owned 140-year-old trophy from the Spanish Pyrenees to confirm its origin. We found that the population of the Pyrenean ibex underwent a recent expansion approximately 20,000 years ago, after which trophy hunting and epizootics triggered a relentless population decline. Our interpretations, based on the genetic information currently available in public repositories, provide a solid basis for more exhaustive analyses relying on all the new sources identified. In particular, the adoption of a genome-wide approach appears a fundamental prerequisite to disentangle the multiple contributing factors associated with low genetic diversity, including inbreeding depression, acting as extinction drivers.

Population Genetic Structure of Ornate Threadfin Bream (Nemipterus hexodon) in Thailand

Ornate threadfin bream (Nemipterus hexodon) is an economically important fishery species in Southeast Asia. In Thailand, N. hexodon decreased dramatically due to overexploitation for commercial purposes. To construct an effective sustainable management plan, genetic information is necessary. Thus, in our study, the population genetic structure and demographic history of N. hexodon were investigated using 419 bp of the mitochondrial DNA sequence in cytochrome oxidase subunit I gene (mtDNA COI). A total of 142 samples was collected from nine localities in the Gulf of Thailand (Chonburi, Samut Songkhram, Surat Thani, Nakhon Si Thammarat, Songkhla), and the Andaman Sea (Satun, Trang, Krabi, Phang Nga). Fourteen polymorphic sites defined 18 haplotypes, revealing a high haplotype diversity and low nucleotide diversity among nine localities. The analysis of molecular variance (AMOVA) analysis, pairwise F _ST , and minimum spanning network result revealed that the genetic structure of N. hexodon was separated into two populations: the Gulf of Thailand and the Andaman Sea population. The genetic structure of N. hexodon can be explained by a disruption of gene flow from the geographic barrier and the Pleistocene isolation of the marine basin hypothesis. Neutrality tests, Bayesian skyline analysis, mismatch distribution, and the estimated values of population growth suggested that N. hexodon had experienced a population expansion. The genetic information would certainly help us gain insight into the population genetic structure of N. hexodon living on the coast of Thailand.

Phylogeography of Cedros and Tiburón Island Mule Deer in North America's Desert Southwest

Though mule deer (Odocoileus hemionus) persist in robust populations throughout most of their North American distribution, habitat loss, unregulated hunting, and other factors have reduced their historical range in México. Two of the 6 putative subspecies inhabiting México's deserts and Baja California peninsula are of conservation concern, occupying islands in the Pacific Ocean (Odocoileus hemionus cerrosensis on Cedros Island: endangered) and Sea of Cortés (Odocoileus hemionus sheldoni on Tiburón Island: threatened). Focusing on the desert southwest (n = 448), we sampled Tiburón (n = 22) and Cedros (n = 15) Island mule deer using contemporary samples and natural history museum specimens to complete a phylogeographic evaluation of the species complex, and assess the phylogeography of these insular subspecies. Both insular subspecies formed endemic haplotype lineages, consistent with island biogeographic theory. Bayesian skyline plots were consistent with Holocene demographic expansion. Cedros Island deer were genetically most similar to adjacent mainland Baja California deer, but exhibited a suite of unique haplotypes and reduced genetic variation. Tiburón Island deer haplotypes unexpectedly nested within a mainland lineage found in distant New Mexico, rather than the adjacent mainland Sonoran lineage. Such findings suggest the importance of postglacial climate fluctuations and biotic community turnover in the phylogeographic history of mule deer in the desert southwest. Our genetic data corroborates cultural, archaeological, and phenotypic evidence supporting Cedros and Tiburón deer endemicity and subspecies status. Reduced genetic variation, divergence from mainland populations, and demographic trends on both islands indicate that conservation, monitoring, and management are critical to ensure persistence of these endemic insular subspecies.

Influence of Quaternary environmental changes on mole populations inferred from mitochondrial sequences and evolutionary rate estimation

Quaternary environmental changes fundamentally influenced the genetic diversity of temperate-zone terrestrial animals, including those in the Japanese Archipelago. The genetic diversity of present-day populations is taxon- and region-specific, but its determinants are poorly understood. Here, we analyzed cytochrome b gene (Cytb) sequences (1140 bp) of mitochondrial DNA (mtDNA) to elucidate the factors determining the genetic variation in three species of large moles: Mogera imaizumii and Mogera wogura, which occur in central and southern mainland Japan (Honshu, Shikoku, and Kyushu), and Mogera robusta, which occurs on the nearby Asian continent. Network construction with the Cytb sequences revealed 10 star-shaped clusters with apparent geographic affinity. Mismatch distribution analysis showed that modes of pairwise nucleotide differences (τ values) were grouped into five classes in terms of the level, implying the occurrence of five stages for rapid expansion. It is conceivable that severe cold periods and subsequent warm periods during the late Quaternary were responsible for the population expansion events. The first and third oldest events included island-derived haplotypes, indicative of the involvement of land bridge formation between remote islands, hence suggesting an association of the ends of the penultimate (PGM, ca. 130,000 years ago) and last (LGM, ca. 15,000 years ago) glacial maxima, respectively. Since the third event was followed by the fourth, it is plausible that the termination of the Younger Dryas and subsequent abrupt warming ca. 11,500 years ago facilitated the fourth expansion event. The second event most likely corresponded to early marine isotope stage (MIS) 3 (ca. 53,000 years ago) when the glaciation and subsequent warming period were predicted to have influenced biodiversity. Utilization of the critical times of 130,000, 53,000, 15,000, and 11,500 years ago as calibration points yielded evolutionary rates of 0.03, 0.045, 0.10 and 0.10 substitutions/site/million years, respectively, showing a time-dependent manner whose pattern was similar to that seen in small rodents reported in our previous studies. The age of the fifth expansion event was calculated to be 5800 years ago with a rate of 0.10 substitutions/site/million years ago during the mid-Holocene, suggestive of the influence of humans or other unspecified reasons, such as the Jomon marine transgression.

Origin, genetic diversity and evolution of Andaman local duck, a native duck germplasm of an insular region of India

Domestic ducks are of paramount importance as a cheap source of protein in rural India. Andaman local duck (ALD) is an indigenous avian genetic resource of Andaman and Nicobar islands (ANI) and is mainly distributed in Middle and Northern parts of these islands. Negligence has brought this breed on the edge of extinction necessitating immediate conservation efforts. Here, we report the genetic diversity, population structure and matrilineal genetic root of ALD. Partial mtDNA D-loop sequences were analyzed in 71 ALD samples and analysis revealed 19 polymorphic sites and 13 haplotypes. Estimated haplotype (Hd ± SD) and nucleotide diversity (π ± SD) were 0.881 ± 0.017 and 0.00897 ± 0.00078 respectively. The high genetic diversity of ALD indicates introgression of genetic material from other local duck breeds. In addition, it can be postulated that ALD bearing high genetic diversity has strong ability to adapt to environmental changes and can withstand impending climate change. Phylogenetic and network analysis indicate that ALD falls under Eurasian clade of mallard and ALD forms three clusters; one cluster is phylogenetically close to Southeast Asian countries, one close to Southern part of mainland India and the third one forms an independent cluster. Therefore, ALD might have migrated either from Southeast Asian countries which enjoy a close cultural bondage with ANI from time immemorial or from Southern part of India. The independent cluster may have evolved locally in these islands and natural selection pressure imposed by environmental conditions might be the driving force for evaluation of these duck haplotypes; which mimics Darwin’s theory of natural selection. The results of the study will be beneficial for formulating future breeding programme and conservation strategy towards sustainable development of the duck breed.

Genetic diversity and phylogenetic analysis of blackbuck (Antilope cervicapra) in southern India

Blackbuck (Antilope cervicapra) is a threatened species endemic to the Indian subcontinent. Many populations of blackbuck are found in southern India. Populations of blackbuck are negatively affected in many places for various reasons, such as habitat destruction and poaching. Their range decreased sharply during the 20th century. There is very limited information available on the population dynamics of blackbuck in southern India. For the phylogenetic and genetic diversity analyses of blackbuck populations among different distribution ranges in southern India, we sequenced mt DNA of cytochrome b (Cyt b) for 120, cytochrome c oxidase subunit-1 (COI) for 137 and the control region (CR) for 137 fecal pellets from eleven different locations in southern India. We analyzed the genetic structure of three mitochondrial markers, the CR, Cyt b and the COI region, separately and in a combined dataset. The haplotype diversity and nucleotide diversity of CR were 0.969 and 0.047, respectively, and were higher than those of Cyt b and COI. A Bayesian phylogeny and an MJ network based on the CR and combined dataset (105 sequences) signified several distinct haplotype clusters within blackbuck, whereas no clusters were identified with the Cyt b and COI phylogenetic analyses. The analysis of molecular variance of the combined data set revealed 52.46% genetic variation within the population. Mismatch distribution analysis revealed that blackbuck populations underwent complex changes with analysis of the combined dataset in each population and analysis of each marker separately in the overall population. The results provide evidence that blackbuck in different geographic locations has a distinct population structure due to habitat fragmentation after the formation of the Western and Eastern Ghats.

House mouse Mus musculus dispersal in East Eurasia inferred from 98 newly determined complete mitochondrial genome sequences

The Eurasian house mouse Mus musculus is useful for tracing prehistorical human movement related to the spread of farming. We determined whole mitochondrial DNA (mtDNA) sequences (ca. 16,000 bp) of 98 wild-derived individuals of two subspecies, M. m. musculus (MUS) and M. m. castaneus (CAS). We revealed directional dispersals reaching as far as the Japanese Archipelago from their homelands. Our phylogenetic analysis indicated that the eastward movement of MUS was characterised by five step-wise regional extension events: (1) broad spatial expansion into eastern Europe and the western part of western China, (2) dispersal to the eastern part of western China, (3) dispersal to northern China, (4) dispersal to the Korean Peninsula and (5) colonisation and expansion in the Japanese Archipelago. These events were estimated to have occurred during the last 2000-18,000 years. The dispersal of CAS was characterised by three events: initial divergences (ca. 7000-9000 years ago) of haplogroups in northernmost China and the eastern coast of India, followed by two population expansion events that likely originated from the Yangtze River basin to broad areas of South and Southeast Asia, including Sri Lanka, Bangladesh and Indonesia (ca. 4000-6000 years ago) and to Yunnan, southern China and the Japanese Archipelago (ca. 2000-3500). This study provides a solid framework for the spatiotemporal movement of the human-associated organisms in Holocene Eastern Eurasia using whole mtDNA sequences, reliable evolutionary rates and accurate branching patterns. The information obtained here contributes to the analysis of a variety of animals and plants associated with prehistoric human migration.

Population differentiation and phylogeography in Lycianthes moziniana (Solanaceae: Capsiceae), a perennial herb endemic to the Mexican Transition Zone

Lycianthes moziniana (Solanaceae: Capsiceae) is a perennial herb with edible fruits that is endemic to Mexico. Three varieties are recognized, all known in the Mexican Transition Zone. Lycianthes moziniana var. margaretiana grows in the Sierra Madre Oriental, whereas L. moziniana var. moziniana is common along the Trans-Mexican Volcanic Belt and in the Sierra Madre Occidental. Lycianthes moziniana var. oaxacana is found exclusively in the Sierra Madre del Sur. The Mexican Transition Zone is a complex geological, climatic and biogeographical area, the result of tectonic and volcanic activity that has promoted genetic divergence and speciation. We determined the genetic variation and structure of L. moziniana. Using phylogeographical approaches, we described the demographic history and evolutionary processes leading its divergence. The intergenic spacers rpl32-trnL and ycf1 were sequenced for 133 individuals pertaining to 15 populations. The genealogical relationships were analysed using haplotype networks. Finally, based on ecological niche models, we inferred the palaeodistribution of L. moziniana during the Pleistocene. The genetic differences and the haplogroups matched the three described varieties. Geological and climatic events of the Mexican Transition Zone facilitated these results. The Trans-Mexican Volcanic Belt isolated the populations of the Sierra Madre Oriental and the Sierra Madre del Sur, while allowing the migration to the Sierra Madre Occidental, during the middle Holocene.

Demographic history and genetic differentiation of an endemic and endangered Ulmus lamellosa (Ulmus)

BACKGROUND

Ulmus lamellosa (one of the ancient species of Ulmus) is an endemic and endangered plant that has undergone climatic oscillations and geographical changes. The elucidation of its demographic history and genetic differentiation is critical for understanding the evolutionary process and ecological adaption to forests in Northern China.

RESULTS

Polymorphic haplotypes were detected in most populations of U. lamellosa via DNA sequencing. All haplotypes were divided into three phylogeographic clades fundamentally corresponding to their geographical distribution, namely THM (Taihang Mountains), YM (Yinshan Mountains), and YSM (Yanshan Mountains) groups. The YSM group, which is regarded as ancestral, possessed higher genetic diversity and significant genetic variability in contrast to the YSM and YM groups. Meanwhile, the divergence time of intraspecies haplotypes occurred during the Miocene-Pliocene, which was associated with major Tertiary geological and/or climatic events. Different degrees of gene exchanges were identified between the three groups. During glaciation, the YSM and THM regions might have served as refugia for U. lamellosa. Based on ITS data, range expansion was not expected through evolutionary processes, except for the THM group. A series of mountain uplifts (e.g., Yanshan Mountains and Taihang Mountains) following the Miocene-Pliocene, and subsequently quaternary climatic oscillations in Northern China, further promoted divergence between U. lamellosa populations.

CONCLUSIONS

Geographical topology and climate change in Northern China played a critical role in establishing the current phylogeographic structural patterns of U. lamellosa. These results provide important data and clues that facilitate the demographic study of tree species in Northern China.

The Zagros Mountains acting as a natural barrier to gene flow in the Middle East: more evidence from the evolutionary history of spiny-tailed lizards (Uromasticinae: Saara)

We investigated the phylogeographic structure of the genus Saara and studied whether the Zagros Mountain ranges influence the gene flow and dispersal at the landscape scale between the Iranian Saara species, including S. asmussi in the central Iranian plateau and S. loricata in the Mesopotamian lowlands. Phylogenetic analyses clearly show three well-supported species, including S. loricata, S. asmussi and S. hardwickii, that are distinct from Uromastyx species. The S-DIVA and BBM analyses demonstrate that species of Saara originated from an ancestor somewhere in the Iranian Plateau and then dispersed to its current geographical range. Our results indicate that the separation of S. loricata from S. asmussi coincides with the orogenic events of the Zagros Mountains during the Pliocene, and thus it may show a vicariance event. Diversification within populations of S. loricata and S. assmussi are estimated to have occurred during the Pleistocene. The haplotype network indicates one haplogroup for each of the Iranian Saara species. Population genetic analyses shows signals of demographic expansions at the beginning of the Holocene for S. loricata. Our results support the hypothesis that the Zagros Mountains act as a barrier for gene flow and Quaternary climatic oscillations affected intraspecific genetic divergences of S. loricata and S. asmussi.

Late Quaternary Environmental and Human Impacts on the Mitochondrial DNA Diversity of Four Commensal Rodents in Myanmar

We addressed the spatiotemporal characteristics of four commensal rodent species occurring in Myanmar in comparison with other areas of the Indo-Malayan region. We examined sequence variations of the mitochondrial cytochrome b gene (Cytb) in the Pacific rat (Rattus exulans), roof rat (Rattus rattus complex, RrC), lesser bandicoot rat (Bandicota bengalensis), and house mouse (Mus musculus) using the recently developed time-dependent evolutionary rates of mtDNA. The Cytb sequences of RrC from Myanmar were shown to belong to RrC Lineage II, and their level of genetic diversity was relatively high compared to those of the other three species. RrC was found to have experienced bottleneck and rapid expansion events at least twice in the late Pleistocene period in Myanmar and a nearby region. Accordingly, paleoclimatic environmental fluctuations were shown to be an important factor affecting rodents in the subtropics of the Indo-Malayan region. Our results show that human activities during the last 10,000 years of the Holocene period affected the population dynamics of the rodent species examined, including introducing them to Myanmar from neighboring countries. Further study of these four commensal rodents in other geographic areas of the Indo-Malayan region would allow us to better understand the factors that drove their evolution and their ecological trends.

Genetic diversity and population structure of Terapon jarbua (Forskål, 1775) (Teleostei, Terapontidae) in Malaysian waters

A background study is important for the conservation and stock management of a species. Terapon jarbua is a coastal Indo-Pacific species, sourced for human consumption. This study examined 134 samples from the central west and east coasts of Peninsular (West) Malaysia and East Malaysia. A 1446-bp concatenated dataset of mtDNA COI and Cyt b sequences was used in this study and 83 haplotypes were identified, of which 79 are unique haplotypes and four are shared haplotypes. Populations of T. jarbua in Malaysia are genetically heterogenous as shown by the high level of haplotype diversity ranging from 0.9167-0.9952, low nucleotide diversity ranging from 0.0288-0.3434, and high FST values (within population genetic variation). Population genetic structuring is not distinct as shown by the shared haplotypes between geographic populations and mixtures of haplotypes from different populations within the same genetic cluster. The gene flow patterns and population structuring observed among these regions are likely attributed to geographical distance, past historical events, allopatric speciation, dispersal ability and water currents. For instance, the mixture of haplotypes revealed an extraordinary migration ability of T. jarbua (>1200 km) via ancient river connectivity. The negative overall value of the neutrality test and a non-significant mismatch distribution are consistent with demographic expansion(s) in the past. The median-joining network concurred with the maximum likelihood haplotype tree with three major clades resolved. The scarcity of information on this species is an obstacle for future management and conservation purposes. Hence, this study aims to contribute information on the population structure, genetic diversity, and historical demography of T. jarbua in Malaysia.

Sky islands as foci for divergence of fig trees and their pollinators in southwest China

The dynamics of populations and their divergence over time have shaped current levels of biodiversity and in the case of the "sky islands" of mountainous southwest (SW) China have resulted in an area of exceptional botanical diversity. Ficus tikoua is a prostrate fig tree subendemic to the area that displays unique intraspecific diversity, producing figs typical of different pollination modes in different parts of its range. By combining climate models, genetic variation in populations of the tree's obligate fig wasp pollinators and distributions of the different plant phenotypes, we examined how this unusual situation may have developed. We identified three genetically distinct groups of a single Ceratosolen pollinator species that have largely parapatric distributions. The complex topography of the region contributed to genetic divergence among the pollinators by facilitating geographical isolation and providing refugia. Migration along elevations in response to climate oscillations further enhanced genetic differentiation of the three pollinator groups. Their distributions loosely correspond to the distributions of the functionally significant morphological differences in the male figs of their host plants, but postglacial expansion of one group has not been matched by spread of its associated plant phenotype, possibly due to a major river barrier. The results highlight how interplay between the complex topography of the "sky island" complex and climate change has shaped intraspecies differentiation and relationships between the plant and its pollinator. Similar processes may explain the exceptional botanical diversity of SW China.

Genetic diversity of a widespread annual killifish from coastal Tanzania

Background

African annual killifishes (Nothobranchius spp.) are adapted to seasonally desiccating habitats (ephemeral pools), surviving dry periods as dormant eggs. Given their peculiar life history, geographic aspects of their diversity uniquely combine patterns typical for freshwater taxa (river basin structure and elevation gradient) and terrestrial animals (rivers acting as major dispersal barriers). However, our current knowledge on fine-scale inter-specific and intra-specific genetic diversity of African annual fish is limited to a single, particularly dry region of their distribution (subtropical Mozambique). Using a widespread annual killifish from coastal Tanzania and Kenya, we tested whether the same pattern of genetic divergence pertains to a wet equatorial region in the centre of Nothobranchius distribution.

Results

In populations of Nothobranchius melanospilus species group across its range, we genotyped a part of mitochondrial cytochrome oxidase subunit 1 (COI) gene (83 individuals from 22 populations) and 10 nuclear microsatellite markers (251 individuals from 16 populations). We found five lineages with a clear phylogeographic structure but frequent secondary contact. Mitochondrial lineages were largely congruent with main population genetic clusters identified on microsatellite markers. In the upper Wami basin, populations are isolated as a putative Nothobranchius prognathus, but include also a population from a periphery of the middle Ruvu basin. Other four lineages (including putative Nothobranchius kwalensis) coexisted in secondary contact zones, but possessed clear spatial pattern. Main river channels did not form apparent barriers to dispersal. The most widespread lineage had strong signal of recent population expansion.

Conclusions

We conclude that dispersal of a Nothobranchius species from a wet part of the genus distribution (tropical lowland) is not constrained by main river channels and closely related lineages frequently coexist in secondary contact zones. We also demonstrate contemporary connection between the Ruvu and Rufiji river basins. Our data do not provide genetic support for existence of recently described cryptic species from N. melanospilus complex, but cannot resolve this issue.

The Anatolian glacial refugium and human-mediated colonization: a phylogeographical study of the stone marten (Martes foina) in Turkey

The Anatolian Peninsula, comprising most of modern Turkey, is situated at the junction of Europe, the Middle East and Asia. Together with its complex geomorphological and climatic history, this has given rise to a rich fauna and flora, which exhibits a wide range of historical biogeographical patterns. The stone marten (Martes foina) is a small carnivore that is found across the temperate Palaearctic region including Anatolia, where it is often associated with habitats modified by humans, but few genetic data exist for this species. We sequenced a 1840-bp region of the mitochondrial genome from 97 martens sampled across the peninsula and intron 7 of the nuclear β-fibrinogen gene from 53 of these. Two mitochondrial lineages were recovered, with overlapping eastern and western distributions, but there was no geographical structure for the autosomal marker. Coalescent analyses indicated that both of the lineages originated during the Last Glacial Maximum, one of them within an eastern Anatolian refugium and the other in a western Anatolian or Balkan refugium. The western lineage colonized most of Europe in the Holocene, while the eastern lineage may be endemic to Anatolia, from where it colonized the Iberian Peninsula via human translocation. The presence of at least one refugial stone marten population highlights the importance of Anatolia to the preservation of genetic variation and biodiversity.

Contrasting patterns of genetic structure and phylogeography in the marine agarophytes Gelidiophycus divaricatus and G. freshwateri (Gelidiales, Rhodophyta) from East Asia

The evolutionary and population demographic history of marine red algae in East Asia is poorly understood. Here, we reconstructed the phylogeographies of two upper intertidal species endemic to East Asia, Gelidiophycus divaricatus and G. freshwateri. Phylogenetic and phylogeographic inferences of 393 mitochondrial cox1, 128 plastid rbcL, and 342 nuclear ITS2 sequences were complemented with ecological niche models. Gelidiophycus divaricatus, a southern species adapted to warm water, is characterized by a high genetic diversity and a strong geographical population structure, characteristic of stable population sizes and sudden reduction to recent expansion. In contrast, G. freshwateri, a northern species adapted to cold temperate conditions, is genetically relatively homogeneous with a shallow population structure resulting from steady population growth and recent equilibrium. The overlap zone of the two species roughly matches summer and winter isotherms, indicating that surface seawater temperature is a key feature influencing species range. Unidirectional genetic introgression was detected at two sites on Jeju Island where G. divaricatus was rare while G. freshwateri was common, suggesting the occurrence of asymmetric natural hybrids, a rarely reported event for rhodophytes. Our results illustrate that Quaternary climate oscillations have left strong imprints on the current day genetic structure and highlight the importance of seawater temperature and sea level change in driving speciation in upper intertidal seaweed species.

Population genetics of Aedes albopictus (Diptera: Culicidae) in its native range in Lao People's Democratic Republic

BACKGROUND

The Asian tiger mosquito, Aedes (Stegomyia) albopictus (Skuse) is an important worldwide invasive species and can be a locally important vector of chikungunya, dengue and, potentially, Zika. This species is native to Southeast Asia where populations thrive in both temperate and tropical climates. A better understanding of the population structure of Ae. albopictus in Lao PDR is very important in order to support the implementation of strategies for diseases prevention and vector control. In the present study, we investigated the genetic variability of Ae. albopictus across a north-south transect in Lao PDR.

METHODS

We used variability in a 1337-bp fragment of the mitochondrial cytochrome c oxidase subunit 1 gene (cox1), to assess the population structure of Ae. albopictus in Lao PDR. For context, we also examined variability at the same genetic locus in samples of Ae. albopictus from Thailand, China, Taiwan, Japan, Singapore, Italy and the USA.

RESULTS

We observed very high levels of genetic polymorphism with 46 novel haplotypes in Ae. albopictus from 9 localities in Lao PDR and Thailand populations. Significant differences were observed between the Luangnamtha population and other locations in Lao PDR. However, we found no evidence of isolation by distance. There was overall little genetic structure indicating ongoing and frequent gene flow among populations or a recent population expansion. Indeed, the neutrality test supported population expansion in Laotian Ae. albopictus and mismatch distribution analyses showed a lack of low frequency alleles, a pattern often seen in bottlenecked populations. When samples from Lao PDR were analyzed together with samples from Thailand, China, Taiwan, Japan, Singapore, Italy and the USA, phylogenetic network and Bayesian cluster analysis showed that most populations from tropical/subtropical regions are more genetically related to each other, than populations from temperate regions. Similarly, most populations from temperate regions are more genetically related to each other, than those from tropical/subtropical regions.

CONCLUSIONS

Aedes albopictus in Lao PDR are genetically related to populations from tropical/subtropical regions (i.e. Thailand, Singapore, and California and Texas in the USA). The extensive gene flow among locations in Lao PDR indicates that local control is undermined by repeated introductions from untreated sites.

Population genetic structure of Texas horned lizards: implications for reintroduction and captive breeding

The Texas horned lizard (Phrynosoma cornutum) inhabits much of the southern Great Plains of North America. Since the 1950s, this species has been extirpated from much of its eastern range and has suffered declines and local extinctions elsewhere, primarily due to habitat loss. Plans are underway to use captive breeding to produce large numbers of Texas horned lizards for reintroduction into areas that were historically occupied by this species and that currently have suitable habitat. We used mitochondrial markers and nuclear microsatellite markers to determine levels of genetic diversity and population structure in 542 Texas horned lizards sampled from across Texas and some neighboring states to help inform these efforts. Texas horned lizards still retain high genetic diversity in many parts of their current range. We found two highly divergent mitochondrial clades (eastern and western) and three major genetic groupings at nuclear microsatellite loci: a west group corresponding to the western mitochondrial clade and north and south groups within the eastern mitochondrial clade. We also found some evidence for human-mediated movement between these genetic clusters that is probably related to the historical importance of this species in the pet trade and as an iconic symbol of the southwestern United States. We do not know, however, if there are fitness costs associated with admixture (especially for the western and eastern clades) or if there are fitness costs to moving these lizards into habitats that are distinctly different from their ancestral areas. If present, either one or both of these fitness costs would decrease the effectiveness of reintroduction efforts. We therefore recommend that reintroduction efforts should maintain current genetic structure by restricting breeding to be between individuals within their respective genetic clusters, and by reintroducing individuals only into those areas that encompass their respective genetic clusters. This cautionary approach is based on the strong divergence between genetic groupings and their correspondence to different ecoregions.

Study on population genetics of Sillago aeolus (Perciformes: Sillaginidae) in the Coast of China

Sillago aeolus is a species from Sillaginidae, which is a widely distributed species with important scientific and economic value in the coast of China. Its population genetics have not been studied. This study investigated the population genetics of S. aeolus in the eastern and southern coast of China based on the mitochondrial control region markers obtained from 248 individuals of 9 locations. The population was characterized by a high population diversity with a low nucleotide diversity. There were no branches corresponding to the sampling sites according to the haplotype network and NJ tree. Recent asymmetric gene flow exchanges and significant genetic differences were detected between the Haikou population and the other populations. AMOVA result indicated slight genetic structures with homogeneity suggested. The neutral test and the mismatch distribution revealed a recent population expansion event. Historical geographic events may be the reason for the homogeneity within the population.

Phylogeography of Dendrolimus punctatus (Lepidoptera: Lasiocampidae): Population differentiation and last glacial maximum survival

Although the Masson pine moth, Dendrolimus punctatus, is one of the most destructive forest pest insects and is an endemic condition in China, we still do not fully understand the patterns of how its distribution range varies in response to Quaternary climatic oscillations. Here, we sequenced one maternally inherited mitochondrial gene (COI) and biparentally inherited nuclear data (ITS1 and ITS2) among 23 natural populations across the entire range of the species in China. A total of 51 mitotypes and 38 ribotypes were separately obtained using mtDNA and ITS1 data. Furthermore, significant phylogeographical structure (N ST > G ST, p < 0.01) were detected. The spatial distribution of mitotypes implied that two distinct groups existed in the species: one in the southwest distribution, including 10 locations, and the other located in the northeast region of China. It is suggested, therefore, that each group was derived from ancestors that occupied different isolated refugia during previous periods, possibly last glacial maximum. Mismatch distribution and Bayesian population dynamics analysis suggested the population size underwent sudden expansion, which is consistent with the results of ecological niche modeling. As a typical phytophagous insect, the history of population expansion was in accordance with the host plants, providing abundant food resources and habitat. Intraspecific success rate of barcoding identification was lower than interspecific ones, indicating a level of difficulty in barcoding individuals from different populations. However, it still provides an early insight into the pattern of genetic diversity within a species.

OPEN RESEARCH BADGES

This article has been awarded an Open Data and Open Materials. All materials and data are publicly accessible via the Open Science Framework at https://doi.org/10.5061/dryad.2df87g2. Learn more about the Open Practices badges from the Center for Open Science: https://osf.io/tvyxz/wiki.

Cryptic diversity and limited connectivity in octopuses: Recommendations for fisheries management

The market demand for octopus grows each year, but landings are decreasing, and prices are rising. The present study investigated (1) diversity of Octopodidae in the Western Indian Ocean (WIO) and (2) connectivity and genetic structure of Octopus cyanea and O. vulgaris populations in order to obtain baseline data for management plans. A fragment of the cytochrome C oxidase subunit 1 (COI) gene was sequenced in 275 octopus individuals from Madagascar, Kenya and Tanzania. In addition, 41 sequences of O. vulgaris from South Africa, Brazil, Amsterdam Island, Tristan da Cunha, Senegal and Galicia were retrieved from databases and included in this study. Five different species were identified using DNA barcoding, with first records for O. oliveri and Callistoctopus luteus in the WIO. For O. cyanea (n = 229, 563 bp), 22 haplotypes were found, forming one haplogroup. AMOVA revealed shallow but significant genetic population structure among all sites (ϕST = 0.025, p = 0.02), with significant differentiation among: (1) Kanamai, (2) southern Kenya, Tanzania, North and West Madagascar, (3) Southwest Madagascar and (4) East Madagascar (ϕCT = 0.035, p = 0.017). For O. vulgaris (n = 71, 482 bp), 15 haplotypes were identified, forming three haplogroups. A significant genetic population structure was found among all sites (ϕST = 0.82, p ≤ 0.01). Based on pairwise ϕST-values and hierarchical AMOVAs, populations of O. vulgaris could be grouped as follows: (1) Brazil, (2) Madagascar and (3) all other sites. A significant increase in genetic distance with increasing geographic distance was found (Z = 232443, 81 r = 0.36, p = 0.039). These results indicate that for O. cyanea four regions should be considered as separate management units in the WIO. The very divergent haplogroups in O. vulgaris from Brazil and Madagascar might be evolving towards speciation and therefore should be considered as separate species in FAO statistics.