Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy.

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

Genetic Population Structure of Summer Flounder Paralichthys dentatus using Microsatellite DNA Analysis.

Summer flounder Paralichthys dentatus supports one of the most valuable commercial and recreational fisheries along the Atlantic Coast of the U.S. However, in recent decades the management of this species has proven to be one of the most contentious for any exploited marine resource in the region. A coastwide catch quota is imposed annually for summer flounder of which 60% is allocated to the commercial fishery and 40% to the recreational fishery. The allocation is further divided among the individual coastal states from North Carolina to Massachusetts based on their landings in the 1980s. This process, based on political jurisdictions, does not consider the species' biological stock structure. Previous genetic studies (allozyme, mtDNA, and SNPs) provided contradictory results regarding the possible population structure of summer. To address this issue, we used DNA microsatellite analysis at 9 loci to define the coastwide population structure of summer flounder. In total, 1,182 specimens were analyzed from 18 collection sites. Most collections were from the continental shelf during the fall-winter spawning season. These were supplemented with additional samples from inshore waters from North Carolina to Florida, and inshore sites which support significant recreational fisheries at Nantucket Shoals, Massachusetts and Fire Island, New York. The overall level of genetic differentiation in pairwise comparison between collections was very low, mean F ST = 0.001. There was no evidence of genetic differentiation between collections from north and south of Cape Hatteras. Our microsatellite results are consistent with an earlier SNP study which failed to find significant allelic heterogeneity among coastwide collections of summer flounder. However, a subset of pairwise F ST comparisons between some collections proved statistically significant. Furthermore, in STRUCTURE analysis we found evidence of two genetic clusters within the species' northern landings area, however, this finding was not supported by DPAC analysis. We conclude that summer flounder most likely constitute a single population along their entire Atlantic Coast distribution.

Genetic population structure of the agaric Blastosporella zonata (Lyophyllaceae) reveals cryptic species and different roles for sexual and asexual spores in dispersal.

Blastosporella zonata is one of the few basidiomycete fungi that produce asexual spores (conidia) on the mushroom. The role of these conidia in the fungal lifecycle is not known. We tested whether conidia are being utilized in local dispersal by looking for signatures of clonality in 21 samples from three localities separated by about three kilometres in Murillo, Colombia. To identify clonally related individuals, we sequenced three polymorphic markers at two unlinked loci (nuclear rRNA: ITS and LSU, and TEFIα) for all collections plus three herbarium samples. We identified two sets of clonally related individuals growing closely together in one of the three localities, and only one pair shared between localities. In all three localities we observed multiple non-clonally related dikaryons showing that sexual reproduction is also important. Our results indicate that the conidia on the mushroom are primarily important for local dispersal. Unexpectedly, our results also indicate two reproductively isolated populations, possibly representing cryptic biological species. Citation: Van de Peppel LJJ, Baroni TJ, Franco-Molano AE, et al. 2022. Genetic population structure of the agaric Blastosporella zonata (Lyophyllacea) reveals cryptic species and different roles for sexual and asexual spores in dispersal. Persoonia 49: 195-200. https://doi.org/10.3767/persoonia.2022.49.06.

Pleistocene Ice Ages Created New Evolutionary Lineages, but Limited Speciation in Northeast Pacific Winged Kelp.

The extent that Pleistocene climate variability promoted speciation has been much debated. Here, we surveyed genetic markers in winged kelp Alaria in the Gulf of Alaska, Northeast Pacific Ocean to understand how paleoclimates may have influenced diversity in this kelp. The study included wide geographic sampling over 2800 km and large sample sizes compared to previous studies of this kelp. Mitochondrial 5'-COI (664 bp), plastid rbcL-3' (740 bp) and 8 microsatellite markers in 16 populations resolved 5 well-defined lineages. COI-rbcL haplotypes were distributed chaotically among populations around the Gulf of Alaska. Principal Coordinates Analysis of microsatellite genotypes grouped plants largely by organellar lineage instead of geography, indicating reproductive isolation among lineages. However, microsatellite markers detected hybrids at 3 sites where lineages co-occurred. Local adaptation on various time scales may be responsible for some genetic differences between populations located along wave-energy and salinity gradients, but the chaotic pattern of variability over hundreds of kilometers is likely due to isolations in northern refugia during Pleistocene ice ages. The range of divergences between populations indicates that episodic glaciations led to the creation of new lineages, but population turnover (local extinctions and recolonizations) limited the formation of new species in the Northeastern Pacific Ocean.

A Vietnamese human genetic variation database.

Large scale human genome projects have created tremendous human genome databases for some well-studied populations. Vietnam has about 95 million people (the 14th largest country by population in the world) of which more than 86% are Kinh people. To date, genetic studies for Vietnamese people mostly rely on genetic information from other populations. Building a Vietnamese human genetic variation database is a must for properly interpreting Vietnamese genetic variants. To this end, we sequenced 105 whole genomes and 200 whole exomes of 305 unrelated Kinh Vietnamese (KHV) people. We also included 101 other previously published KHV genomes to build a Vietnamese human genetic variation database of 406 KHV people. The KHV database contains 24.81 million variants (22.47 million single nucleotide polymorphisms (SNPs) and 2.34 million indels) of which 0.71 million variants are novel. It includes more than 99.3% of variants with a frequency of >1% in the KHV population. Noticeably, the KHV database revealed 107 variants reported in the human genome mutation database as pathological mutations with a frequency above 1% in the KHV population. The KHV database (available at https://genomes.vn) would be beneficial for genetic studies and medical applications not only for the Vietnamese population but also for other closely related populations.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium.

BACKGROUND: The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining. RESULTS: We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination. CONCLUSIONS: The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

The genetic population structure and temporal genetic stability of gilthead sea bream Sparus aurata populations in the Aegean and Ionian Seas, using microsatellite DNA markers.

We examined 662 gilthead sea bream Sparus aurata from wild samples of the species in the Aegean and Ionian Seas, using 20 EST-linked microsatellite markers, in three multiplex panels, as well as seven anonymous loci. Most of the markers were revealed to be highly polymorphic. We found low genetic differentiation between the sampling stations/areas with total FST 0.002 (P < 0.05). Based on comparison of five temporal samples, our results indicate genetic data consistency over time for all tested samples, pointing to stable populations, despite reported repeated escape events. Our results confirm the genetic population structure previously observed in these specific areas, using by far more markers than in previous studies in both coding and non-coding DNA loci. The limited genetic structure and the temporal genetic stability indicate neither major genetic differentiation of local populations by geographic isolation nor influence from anthropogenic factors. These results provide a baseline for future reference in any management programme of both wild and farmed population of S. aurata as well as of other aquaculture species with a potential introgression among farmed and wild populations.

Spatial variation in anthropogenic mortality induces a source-sink system in a hunted mesopredator.

Lethal carnivore management is a prevailing strategy to reduce livestock predation. Intensity of lethal management varies according to land-use, where carnivores are more intensively hunted on farms relative to reserves. Variations in hunting intensity may result in the formation of a source-sink system where carnivores disperse from high-density to low-density areas. Few studies quantify dispersal between supposed sources and sinks-a fundamental requirement for source-sink systems. We used the black-backed jackal (Canis mesomelas) as a model to determine if heterogeneous anthropogenic mortality induces a source-sink system. We analysed 12 microsatellite loci from 554 individuals from lightly hunted and previously unhunted reserves, as well as heavily hunted livestock- and game farms. Bayesian genotype assignment showed that jackal populations displayed a hierarchical population structure. We identified two genetically distinct populations at the regional level and nine distinct subpopulations at the local level, with each cluster corresponding to distinct land-use types separated by various dispersal barriers. Migration, estimated using Bayesian multilocus genotyping, between reserves and farms was asymmetric and heterogeneous anthropogenic mortality induced source-sink dynamics via compensatory immigration. Additionally some heavily hunted populations also acted as source populations, exporting individuals to other heavily hunted populations. This indicates that heterogeneous anthropogenic mortality results in the formation of a complex series of interconnected sources and sinks. Thus, lethal management of mesopredators may not be an effective long-term strategy in reducing livestock predation, as dispersal and, more importantly, compensatory immigration may continue to affect population reduction efforts as long as dispersal from other areas persists.

Population Connectivity in the Common Reef Zoantharian Zoanthus sansibaricus (Anthozoa: Hexacorallia) in Southern Japan.

Tropical and subtropical shallow benthic marine communities are highly diverse and balanced systems that constitute an important natural resource. Knowledge of the genetic diversity, connectivity and reproduction mode of each population is critical to understanding the fate of whole assemblages in times of disturbances. Importantly, the capability of populations to adapt to environmental challenges will be crucial to determining their survival. Here, we report on the population structure of the common reef zoantharian Zoanthus sansibaricus in the northwestern Pacific, by examining populations at three different locations in southern Japan using five highly variable microsatellite markers. Analyses of a population at the species' northern distribution limit combined with analyses of two subtropical populations suggest that habitat characteristics and ocean currents influence the connectivity and genetic diversity of this species. Our findings emphasize the adaptive ability of Z. sansibaricus to different environmental conditions and may help explain the wide distribution and generalist nature of this species.

Global population genetic dynamics of a highly migratory, apex predator shark.

Knowledge of genetic connectivity dynamics in the world's large-bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top-down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo-Pacific ocean basins (microsatellite FST  > 0.129; CR ΦST  > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex-biased gene flow. Coalescent- and genetic distance-based estimates of divergence from CR sequences were largely congruent (dcorr  = 0.0015-0.0050), indicating a separation of Indo-Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller-than-anticipated spatial scales.

Conservation genomics of natural and managed populations: building a conceptual and practical framework.

The boom of massive parallel sequencing (MPS) technology and its applications in conservation of natural and managed populations brings new opportunities and challenges to meet the scientific questions that can be addressed. Genomic conservation offers a wide range of approaches and analytical techniques, with their respective strengths and weaknesses that rely on several implicit assumptions. However, finding the most suitable approaches and analysis regarding our scientific question are often difficult and time-consuming. To address this gap, a recent workshop entitled 'ConGen 2015' was held at Montana University in order to bring together the knowledge accumulated in this field and to provide training in conceptual and practical aspects of data analysis applied to the field of conservation and evolutionary genomics. Here, we summarize the expertise yield by each instructor that has led us to consider the importance of keeping in mind the scientific question from sampling to management practices along with the selection of appropriate genomics tools and bioinformatics challenges.

Update of phylogenetic and genetic diversity of Sporothrix schenckii sensu lato.

Sporothrix schenckii sensu lato causes subcutaneous mycosis. In this article we analysed its phylogeny and genetic diversity using calmodulin DNA sequences deposited in GenBank database. Population genetics indices were calculated, plus phylogenetic and haplotype network trees were built. Five clades with high values of posterior probability, 47 haplotypes and high diversity in the complex were found. Analysis of partial calmodulin sequences alignment revealed conserved and polymorphic regions that could be used as reference for taxonomic identification. The use of population genetics analysis allowed understanding the phylogenetic proximity of S. schenckii s. str. and S. brasiliensis; scarce genetic flow among them with low migration index and high ancestry coefficient was found. Similarly, S. globosa, S. mexicana and S. pallida sequences showed highly differentiated species with no genetic exchange. The phylogenetic tree suggests that S. mexicana shared a common ancestor with S. pallida; while S. globosa and S. brasiliensis are more related to S. schenckii s. str. and showed less haplotype diversity and restrictions in geographic distribution. In the haplotype network tree S. schenckii s. str. species displayed worldwide distribution without dispersion centres; while S. brasiliensis and S. globosa, exhibited Brazil and Euro-Asia as dispersion centres, respectively. Our data suggest that S. schenckii complex has been submitted to a divergent evolution process, probably due to the pressure of the environment and of the host. In contrast, S. brasiliensis could have been submitted to purifying selection or expansion process.

New evidence for hybrid zones of forest and savanna elephants in Central and West Africa.

The African elephant consists of forest and savanna subspecies. Both subspecies are highly endangered due to severe poaching and habitat loss, and knowledge of their population structure is vital to their conservation. Previous studies have demonstrated marked genetic and morphological differences between forest and savanna elephants, and despite extensive sampling, genetic evidence of hybridization between them has been restricted largely to a few hybrids in the Garamba region of northeastern Democratic Republic of Congo (DRC). Here, we present new genetic data on hybridization from previously unsampled areas of Africa. Novel statistical methods applied to these data identify 46 hybrid samples--many more than have been previously identified--only two of which are from the Garamba region. The remaining 44 are from three other geographically distinct locations: a major hybrid zone along the border of the DRC and Uganda, a second potential hybrid zone in Central African Republic and a smaller fraction of hybrids in the Pendjari-Arli complex of West Africa. Most of the hybrids show evidence of interbreeding over more than one generation, demonstrating that hybrids are fertile. Mitochondrial and Y chromosome data demonstrate that the hybridization is bidirectional, involving males and females from both subspecies. We hypothesize that the hybrid zones may have been facilitated by poaching and habitat modification. The localized geography and rarity of hybrid zones, their possible facilitation from human pressures, and the high divergence and genetic distinctness of forest and savanna elephants throughout their ranges, are consistent with calls for separate species classification.

Genetic variation and population structure of American mink Neovison vison from PCB-contaminated and non-contaminated locales in eastern North America.

American mink Neovison vison may be particularly vulnerable to toxicities of persistent contaminants such as PCBs because of their aquatic-based diet, position near the top of the food web, and small deme sizes. Furthermore, ranched mink are sensitive to reproductive toxicities of fish diets from PCB-polluted sites. The upper Hudson River is highly contaminated with PCBs and previous studies have shown elevated hepatic burdens of total and coplanar PCBs in mink collected near the river compared with those from more distant locales in New York and elsewhere. We hypothesized that bioaccumulation of PCBs in Hudson River mink has reduced their levels of genetic diversity or altered their genetic population structure. To address this, we conducted microsatellite DNA analysis on collections made in proximity to and from more distant locales in the Hudson River watershed, elsewhere in New York State, and at other sites in eastern North America including New Brunswick, four locales in Ontario, multiple drainages in Maine, and two ecoregions in Rhode Island. We did not find reduced genetic diversity at the individual or population levels in mink collected near (<6 km) to PCB hotspots in the Hudson River nor evidence of altered population structure. Consistent with their distribution in small localized and isolated demes, we did find significant genetic population structure among many mink collections in New York State and elsewhere. Depending on the analytical approach used, genetically distinct populations numbered between 16 when using STRUCTURE to 19-20 when using Exact G tests, F ST, or AMOVA analyses. Genetically distinct population units were found among major ecoregions and minor ecoregions in New York State, among different hydrologic subunits within the Hudson River watershed, among spatially separate locales in Ontario, and among most watersheds in Maine. However, despite this localization and potential heightened impact of stressors, genetic diversity and genetic population structure in mink does not seem to be affected by their bioaccumulation of high levels of PCBs of Hudson River origin.

Wild yeast isolation by middle school students reveals features of North American oak populations of Saccharomyces cerevisiae and Kluyveromyces lactis.

Features of the natural life cycle of the budding yeast Saccharomyces cerevisiae were crucial to its domestication as a laboratory experimental model, especially the ability to maintain stable haploid clones and cross them at will to combine alleles via meiosis. Stable haploidy results from mutations in HO, which encodes an endonuclease required for haploid-specific mating-type switching. Previous studies found an unexpected diversity of HO alleles among natural isolates within a small geographic area. We developed a hands-on field and laboratory activity for middle school students in Denver, Colorado, USA to isolate wild yeast from oak bark, identify species via DNA sequencing, and sequence HO from S. cerevisiae isolates. We find limited HO diversity in North American oak isolates, pointing to efficient, continuous dispersal across the continent. By contrast, we isolated the "dairy yeast", Kluyveromyces lactis, from a tree <10 m away and found that it represents a new population distinct from an oak population in an adjacent state, pointing to high genetic diversity. The outreach activity partnered middle school, high school, and university students in making scientific discoveries and can be adapted to other locations and natural yeast habitats. Indeed, a pilot sampling activity in southeast Texas yielded S. cerevisiae oak isolates with a new allele of HO and, from a nearby prickly pear cactus, a heat-tolerant isolate of Saccharomyces paradoxus.

Phylodynamics unveils invading and diffusing patterns of dengue virus serotype-1 in Guangdong, China from 1990 to 2019 under a global genotyping framework.

BACKGROUND: The strong invasiveness and rapid expansion of dengue virus (DENV) pose a great challenge to global public health. However, dengue epidemic patterns and mechanisms at a genetic scale, particularly in term of cross-border transmissions, remain poorly understood. Importation is considered as the primary driver of dengue outbreaks in China, and since 1990 a frequent occurrence of large outbreaks has been triggered by the imported cases and subsequently spread to the western and northern parts of China. Therefore, this study aims to systematically reveal the invasion and diffusion patterns of DENV-1 in Guangdong, China from 1990 to 2019. METHODS: These analyses were performed on 179 newly assembled genomes from indigenous dengue cases in Guangdong, China and 5152 E gene complete sequences recorded in Chinese mainland. The genetic population structure and epidemic patterns of DENV-1 circulating in Chinese mainland were characterized by phylogenetics, phylogeography, phylodynamics based on DENV-1 E-gene-based globally unified genotyping framework. RESULTS: Multiple serotypes of DENV were co-circulating in Chinese mainland, particularly in Guangdong and Yunnan provinces. A total of 189 transmission clusters in 38 clades belonging to 22 subgenotypes of genotype I, IV and V of DENV-1 were identified, with 7 Clades of Concern (COCs) responsible for the large outbreaks since 1990. The epidemic periodicity was inferred from the data to be approximately 3 years. Dengue transmission events mainly occurred from Great Mekong Subregion-China (GMS-China), Southeast Asia (SEA), South Asia Subcontinent (SASC), and Oceania (OCE) to coastal and land border cities respectively in southeastern and southwestern China. Specially, Guangzhou was found to be the most dominant receipting hub, where DENV-1 diffused to other cities within the province and even other parts of the country. Genome phylogeny combined with epidemiological investigation demonstrated a clear local consecutive transmission process of a 5C1 transmission cluster (5C1-CN4) of DENV-1 in Guangzhou from 2013 to 2015, while the two provinces of Guangdong and Yunnan played key roles in ongoing transition of dengue epidemic patterns. In contextualizing within Invasion Biology theories, we have proposed a derived three-stage model encompassing the stages of invasion, colonization, and dissemination, which is supposed to enhance our understanding of dengue spreading patterns. CONCLUSIONS: This study demonstrates the invasion and diffusion process of DENV-1 in Chinese mainland within a global genotyping framework, characterizing the genetic diversities of viral populations, multiple sources of importation, and periodic dynamics of the epidemic. These findings highlight the potential ongoing transition trends from epidemic to endemic status offering a valuable insight into early warning, prevention and control of rapid spreading of dengue both in China and worldwide.

The genetic structure of populations of Isthmiophora melis (Schrank, 1788) (Digenea: Echinostomatidae). Does the host's diet matter?

BACKGROUND: Here we provide a comparative analysis of the genetic structure of populations (based on nad1 mtDNA) of Isthmiophora melis isolated from the American mink (Neogale vison), an introduced invasive species, commonly occurring in the territory of Poland, and from the striped field mouse (Apodemus agrarius). METHODS: A total of 133 specimens of I. melis were obtained from naturally infected N. vison collected from six localities in Poland (108 samples) and 25 individuals of I. melis from A. agrarius. All sequences of the nad1 gene obtained during the present study were assembled and aligned. The standard statistics for haplotype composition, i.e., the number of haplotypes, haplotype diversity, nucleotide diversity, and average number of nucleotide differences, were calculated. Haplotype analysis and visualization of haplotype frequency among populations were performed using a median-joining network. RESULTS: Based on the samples collected from different localities in Poland, our study revealed that the overall genetic diversity of I. melis isolated from the American mink and of the striped field mouse do not differ significantly. The median-joining network showed that the three main haplotypes are in the centre of a star-like structure, with the remaining haplotypes as the satellites, reflecting the recent expansion of the populations. CONCLUSIONS: The overall genetic diversity of I. melis isolated from the American mink and striped field mouse reveals a high level of homogeneity. Moreover, regional differences in the food composition of the definitive hosts play an important role in shaping the genetic structure of the trematode populations.

ddRAD Sequencing Identifies Pesticide Resistance-Related Loci and Reveals New Insights into Genetic Structure of Bactericera cockerelli as a Plant Pathogen Vector.

(1) Background: Many hemipteran insects transmit plant pathogens that cause devastating crop diseases, while pest management frequently relies primarily on insecticide applications. These intense insecticide applications lead to the development of insecticide resistance, as was the case for potato psyllid, Bactericera cockerelli (Hemiptera: Triozidae), a vector of Candidatus Liberibacter solanacearum, which causes zebra chip disease in potato. (2) Methods: Here, we use double-digest restriction site-associated DNA (ddRAD) to genotype eight psyllid populations (one susceptible and seven resistant to neonicotinoid insecticides). (3) Results: Association tests identified over 400 loci that were strongly segregated between susceptible and resistant populations. Several loci were located within genes involved in insecticide resistance, gene regulation, fertility, and development. Moreover, we explored the genetic structure of these eight populations and discovered that routinely utilized haplotyping was not an accurate predictor of population structure. Pairwise comparisons of the fixation index (FST) of populations of the same haplotype were not different from pairwise FST of populations that belonged to different haplotypes. (4) Conclusions: Our findings suggest that neonicotinoid insecticide resistance has a genetic basis, most likely as a result of similar selection pressure. Furthermore, our results imply that using a single maternally inherited gene marker to designate genetic lineages for potato psyllids should be re-evaluated.

Population Structure of German Cockroaches (Blattodea: Ectobiidae) in an Urban Environment Based on Single Nucleotide Polymorphisms.

German cockroaches (Blattella germanica L.) harbor and disperse medically important pathogens and are a source of allergens that impact human health and wellbeing. Management of this pest requires an understanding of their distribution and dispersal. In this study, we collected German cockroaches from three apartment buildings in New Jersey, USA. We identified single-nucleotide polymorphisms (SNPs) from DNA extractions using next generation sequencing. We analyzed the SNPs and characterized cockroach population genetic structure using Fst, principal component, phylogenetic, and STRUCTURE analyses. We found significant differences in German cockroach population structure among the buildings. Within buildings, we found variable population structure that may be evidence for multiple colonization events. This study shows that SNPs derived from next generation sequencing provide a powerful tool for analyzing the genetic population structure of these medically important pests.