Invasive species’ Pleistocene history contains much information on its present population structure, dispersability and adaptability. In this study, the Pleistocene history of a global invasive pest (Brown Marmorated Stink Bug BMSB, Halyomorpha halys) was unveiled using the coupled approach of phylogeography and ecological niche modelling. Rangewide molecular data suggests that the Taiwan and other native populations had diverged in mid-Pleistocene. In mainland China, the native BMSB did not experience population contraction and divergence during last glacial, but persisted in interconnected populations. Combined Bayesian Skyline Plot (BSP) and niche modelling revealed a rapid expansion occurred during the transition of Last Inter Glacial (LIG) to Last Glacial Maximum (LGM). High genetic diversity and multi-reticular haplotypes network exist in the original sources populations of BMSB invasion in northern China. They were speculated to be colonized from the central China, with many derived haplotypes evolved to adapt the novel environment. The ENM future prediction suggest that BMSB may expand northward to higher latitudes in the US and Europe, because of its high invasive ability, together with the available suitable climate space there.

BackgroundOn account of repeated exposure and submergence of the East China Sea (ECS) land bridge, sea level fluctuation played an important role in shaping the population structure of many temperate species across the ECS during the glacial period. The flower bug Amphiareus obscuriceps (Poppius, 1909) (Hemiptera: Anthocoridae) is an invasive species native to the Sino-Japanese Region (SJR) of East Asia. We tested the hypothesis of the ECS land bridge acting as a dispersal corridor or filter for A. obscuriceps during the glacial period. Specifically, we tested whether and the extent to which dispersal ability and host plant habitat requirement influenced the genetic structure of A. obscuriceps during the exposure of the ECS land bridge.ResultsPhylogenetic and network analyses indicated that A. obscuriceps is composed of two major lineages, i.e., China and Japan. Divergence time on both sides of the ECS was estimated to be approximately 1.07 (0.79–1.32) Ma, which was about the same period that the sea level increased. No significant Isolation by Distance (IBD) relationship was found between Фst and Euclidean distances in the Mantel tests, which is consistent with the hypothesis that this species has a good dispersal ability. Our Last Glacial Maximum (LGM) niche modeling of plants that constitute preferred habitats for A. obscuriceps exhibited a similar habitat gap on the exposed ECS continental shelf between China and Japan, but showed a continuous distribution across the Taiwan Strait.ConclusionOur results suggest that ecological properties (habitat requirement and dispersal ability), together with sea level fluctuation during the Pleistocene across the ECS, have shaped the genetic structure and demographic history of A. obscuriceps in its native area. The host plant habitat requirement could also be a key to the colonization of the A. obscuriceps species during the exposure of the ECS land bridge. Our findings will shed light on the potential role of habitat requirement in the process of biological invasion in future studies.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-016-0748-3) contains supplementary material, which is available to authorized users.

Aim: Local environmental selection and lineage admixture have long been accepted as important adaptive mechanisms in adjusting widespread taxa to new environments.We studied a pond skater, Gerris latiabdominis, to explore the relative roles played by these two mechanisms in its process of adaptation to heterogeneous landscapes. Location: East Asia.Taxon: Gerris latiabdominis Miyamoto, 1958.Methods: Mitochondrial (COI, COII) and nuclear (ddRAD-seq) markers were sequenced from 202 individuals. Phylogenetic and phylogeographical analyses were conducted to reveal the population genetic structure. The demographical history was simulated by approximate Bayesian computation (ABC). The single-nucleotide polymorphisms (SNPs) that may have been subjected to natural divergent selection among the populations were assessed, and the recent migration rates were estimated. Finally, shifts in suitable habitat from the last interglacial (LIG) to the present were predicted through ecological niche modelling (ENM). Results:We found population structures inferred from both mitochondrial and nuclear genomes to be almost coincident and to correspond to temporal or spatial heterogeneity. The early Pliocene splitting event between the northern (NO) and southeastern (SE) groups is consistent with the 'arid belt hypothesis', followed by local environmental selection. Approximate Bayesian computation and gene flow estimation strongly suggested that the southwestern (SW) groups originated from a lineage admixture event between the SE and western (WE) groups and subsequently underwent ongoing gene flow from the SE group during the Pleistocene. The mitonuclear discordance pattern in the SW group was probably due to the ancient mitochondrial differentiation. Main conclusion:Our results suggest that both local environmental selection and lineage admixture acted as significant mechanisms of G. latiabdominis adaptation to | 1155 YE Et al. B I OS K E TCHZhen Ye is a researcher at College of Life Sciences, Nankai University. His research interests are phylogeography, speciation and niche model building of aquatic insects.

Gondwanan vicariance, long‐distance dispersal (LDD), and boreotropical migration have been proposed as alternative hypotheses explaining the pantropical distribution pattern of organisms. In this study, the historical biogeography of the pond skater genus Limnogonus was reconstructed to evaluate the impact of biogeographical scenarios in shaping their modern transoceanic disjunction. We sampled almost 65% of recognized Limnogonus species. Four DNA fragments including 69 sequences were used to reconstruct a phylogram. Divergence time was estimated using a Bayesian relaxed clock method and three fossil calibrations. Diversification dynamics and ancestral area reconstruction were investigated by using maximum likelihood and Bayesian approaches. Our results showed the crown group of Limnogonus originated and diversified in Africa in the early Eocene (49 Ma, HPD: 38–60 Ma), subsequently expanding into other regions via dispersal. The colonization of the New World originated from the Oriental Region probably via the Bering Land Bridge in the late Eocene. Two split events between the Old World and New World were identified: one between Neotropics and Oriental region around the middle Oligocene (30 Ma, HPD: 22–38 Ma), and the other between Neotropics and Africa during the middle Miocene (14 Ma, HPD: 8–21 Ma). The evolutionary history of Limnogonus involved two biogeographical processes. Gondwanan vicariance was not supported in our analyses. The diversification of Limnogonus among Africa, Oriental, and Neotropical regions corresponded with the age of land bridge connection and dispersed as a member associated with the broad boreotropical belt before local cooling (34 Ma). The current transoceanic disjunctions in Limnogonus could be better explained by the disruption of “mixed‐mesophytic” forest belt; however, the direct transoceanic LDD between the Neotropics and Africa could not be ruled out. In addition, the “LDD” model coupled with island hopping could be a reasonable explanation for the diversification of the Oriental and Australian regions during the Oligocene.

Objectives Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area. Methods Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information. Results M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%). Conclusion Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.