Geographic and haplotype structure of candidate type 2 diabetes susceptibility variants at the calpain-10 locus

Recently, a positional cloning study proposed that haplotypes at the calpain-10 locus (CAPN10) are associated with increased risk of type 2 diabetes, or non-insulin-dependent diabetes mellitus, in Mexican Americans, Finns, and Germans. To inform the interpretation of the original mapping results and to look for evidence for the action of natural selection on CAPN10, we undertook a population-based genotyping survey of the candidate susceptibility variants. First, we genotyped sites 43, 19, and 63 (the haplotype-defining variants previously proposed) and four closely linked SNPs, in 561 individuals from 11 populations from five continents, and we examined the linkage disequilibrium among them. We then examined the ancestral state of these sites by sequencing orthologous portions of CAPN10 in chimpanzee and orangutan (the identity of sites 43 and 19 was further investigated in a limited sample of other great apes and Old World and New World monkeys). Our survey suggests larger-than-expected differences in the distribution of CAPN10 susceptibility variants between African and non-African populations, with common, derived haplotypes in European and Asian samples (including one of two proposed risk haplotypes) being rare or absent in African samples. These results suggest a history of positive natural selection at the locus, resulting in significant geographic differences in polymorphism frequencies. The relationship of these differences to disease risk is discussed.     

Human skin color diversity is highest in sub-Saharan African populations

Previous studies of genetic and craniometric traits have found higher levels of within-population diversity in sub-Saharan Africa compared to other geographic regions. This study examines regional differences in within-population diversity of human skin color. Published data on skin reflectance were collected for 98 male samples from eight geographic regions: sub-Saharan Africa, North Africa, Europe, West Asia, Southwest Asia, South Asia, Australasia, and the New World. Regional differences in local within-population diversity were examined using two measures of variability: the sample variance and the sample coefficient of variation. For both measures, the average level of within-population diversity is higher in sub-Saharan Africa than in other geographic regions. This difference persists even after adjusting for a correlation between within-population diversity and distance from the equator. Though affected by natural selection, skin color variation shows the same pattern of higher African diversity as found with other traits.     

Linkage disequilibrium among commonly genotyped SNP variants detected from bull sequence,

Genomic prediction utilizing causal variants could increase selection accuracy above that achieved with SNPs genotyped by currently available arrays used for genomic selection. A number of variants detected from sequencing influential sires are likely to be causal, but noticeable improvements in prediction accuracy using imputed sequence variant genotypes have not been reported. Improvement in accuracy of predicted breeding values may be limited by the accuracy of imputed sequence variants. Using genotypes of SNPs on a high‐density array and non‐synonymous SNPs detected in sequence from influential sires of a multibreed population, results of this examination suggest that linkage disequilibrium between non‐synonymous and array SNPs may be insufficient for accurate imputation from the array to sequence. In contrast to 75% of array SNPs being strongly correlated to another SNP on the array, less than 25% of the non‐synonymous SNPs were strongly correlated to an array SNP. When correlations between non‐synonymous and array SNPs were strong, distances between the SNPs were greater than separation that might be expected based on linkage disequilibrium decay. Consistently near‐perfect whole‐genome linkage disequilibrium between the full array and each non‐synonymous SNP within the sequenced bulls suggests that whole‐genome approaches to infer sequence variants might be more accurate than imputation based on local haplotypes. Opportunity for strong linkage disequilibrium between sequence and array SNPs may be limited by discrepancies in allele frequency distributions, so investigating alternate genotyping approaches and panels providing greater chances of frequency‐matched SNPs strongly correlated to sequence variants is also warranted. Genotypes used for this study are available from https://www.animalgenome.org/repository/pub/ ;USDA2017.0519/.     

Fixation of the human-specific CMP-N-acetylneuraminic acid hydroxylase pseudogene and implications of haplotype diversity for human evolution

The human CMP-N-acetylneuraminic acid hydroxylase gene (CMAH) suffered deletion of an exon that encodes an active center for the enzyme approximately 3.2 million years ago (MYA). We analyzed a 7.3-kb intronic region of 132 CMAH genes to explore the fixation process of this pseudogene and the demographic implication of its haplotype diversity. Fifty-six variable sites were sorted into 18 different haplotypes with significant linkage disequilibrium. Despite the rather low nucleotide diversity, the most recent common ancestor at CMAH dates to 2.9 MYA. This deep genealogy follows shortly after the original exon deletion, indicating that the deletion has fixed in the population, although whether this fixation was facilitated by natural selection remains to be resolved. Remarkable features are exceptionally long persistence of two lineages and the confinement of one lineage in Africa, implying that some African local populations were in relative isolation while others were directly involved in multiple African exoduses of the genus Homo. Importantly, haplotypes found in Eurasia suggest interbreeding between then-contemporaneous human species. Although population structure within Africa complicates the interpretation of phylogeographic information of haplotypes, the data support a single origin of modern humans, but not with complete replacement of archaic inhabitants by modern humans.     

Recombination networks as genetic markers in a human variation study of the Old World

We have analyzed human genetic diversity in 33 Old World populations including 23 populations obtained through Genographic Project studies. A set of 1,536 SNPs in five X chromosome regions were genotyped in 1,288 individuals (mostly males). We use a novel analysis employing subARG network construction with recombining chromosomal segments. Here, a subARG is constructed independently for each of five gene-free regions across the X chromosome, and the results are aggregated across them. For PCA, MDS and ancestry inference with STRUCTURE, the subARG is processed to obtain feature vectors of samples and pairwise distances between samples. The observed population structure, estimated from the five short X chromosomal segments, supports genome-wide frequency-based analyses: African populations show higher genetic diversity, and the general trend of shared variation is seen across the globe from Africa through Middle East, Europe, Central Asia, Southeast Asia, and East Asia in broad patterns. The recombinational analysis was also compared with established methods based on SNPs and haplotypes. For haplotypes, we also employed a fixed-length approach based on information-content optimization. Our recombinational analysis suggested a southern migration route out of Africa, and it also supports a single, rapid human expansion from Africa to East Asia through South Asia.     

Linkage disequilibrium and recent selection at three immunity receptor loci in Drosophila simulans

Immune system genes in a California population sample of Drosophila simulans were shown to bear several hallmarks of the effects of past directional selection. One potential effect of directional selection is an increase in linkage disequilibrium among the polymorphic sites that are linked to the site under selection. In this study, we focus on three D. simulans immunity loci, Hmu, Sr-CI/Sr-CIII, and Tehao, for which the polymorphic sites are in nearly perfect linkage disequilibrium, an unusual finding even with respect to other immunity genes sampled from the same lines. The most likely explanation for this finding is that, at each locus, two divergent alleles have been selected to intermediate frequencies in the recent past. The extent to which the linkage disequilibrium extends to the flanks of each of the immunity genes is minimal, suggesting that the favored mutations actually occurred within the immunity genes themselves. Furthermore, the excess linkage disequilibrium found in the California population is not found in an African D. simulans population sample and may be a result of novel pathogen-mediated selection pressures encountered during establishment of non-African populations.     

Natural selection and population history in the human angiotensinogen gene (AGT): 736 complete AGT sequences in chromosomes from around the world

Several lines of evidence suggest that patterns of genetic variability in the human angiotensinogen gene (AGT) contribute to phenotypic variability in human hypertension. The A(-6) promoter variant of AGT is associated with higher plasma angiotensinogen levels and increased risk of essential hypertension. The geographic distribution of the A(-6) variant leads to the intriguing hypothesis that the G(-6) promoter variant has been selectively advantageous outside Africa. To test these hypotheses, we investigated the roles of population history and natural selection in shaping patterns of genetic diversity in AGT, by sequencing the entire AGT gene (14400 bp) in 736 chromosomes from Africa, Asia, and Europe. We found that the A(-6) variant is present at higher frequency in African populations than in non-African populations. Neutrality tests found no evidence of a departure from selective neutrality, when whole AGT sequences were compared. However, tests restricted to sites in the vicinity of the A(-6)G polymorphism found evidence of a selective sweep. Sliding-window analyses showed that evidence of the sweep is restricted to sites in tight linkage disequilibrium (LD) with the A(-6)G polymorphism. Further, haplotypes carrying the G(-6) variant showed elevated levels of LD, suggesting that they have risen recently to high frequency. Departures from neutral expectation in some but not all regions of AGT indicate that patterns of diversity in the gene cannot be accounted for solely by population history, which would affect all regions equally. Taken together, patterns of genetic diversity in AGT suggest that natural selection has generally favored the G(-6) variant over the A(-6) variant in non-African populations. However, important localized effects may also be present.     

Recombination gives a new insight in the effective population size and the history of the old world human populations

The information left by recombination in our genomes can be used to make inferences on our recent evolutionary history. Specifically, the number of past recombination events in a population sample is a function of its effective population size (Ne). We have applied a method, Identifying Recombination in Sequences (IRiS), to detect specific past recombination events in 30 Old World populations to infer their Ne. We have found that sub-Saharan African populations have an Ne that is approximately four times greater than those of non-African populations and that outside of Africa, South Asian populations had the largest Ne. We also observe that the patterns of recombinational diversity of these populations correlate with distance out of Africa if that distance is measured along a path crossing South Arabia. No such correlation is found through a Sinai route, suggesting that anatomically modern humans first left Africa through the Bab-el-Mandeb strait rather than through present Egypt.     

A family-based association study identified CYP17 as a candidate gene for obesity susceptibility in Caucasians

The cytochrome P450c17α gene (CYP17) encodes a key biosynthesis enzyme of estrogen, which is critical in regulating adipogenesis and adipocyte development in humans. We therefore hypothesized that CYP17 is a candidate gene for predicting obesity. In order to test this hypothesis, we performed a family-based association test to investigate the relationship between the CYP17 gene and obesity phenotypes in a large sample comprising 1873 subjects from 405 Caucasian nuclear families of European origin recruited by the Osteoporosis Research Center of Creighton University, USA. Both single SNPs and haplotypes were tested for associations with obesity-related phenotypes, including body mass index (BMI) and fat mass. We identified three SNPs to be significantly associated with BMI, including rs3740397, rs6163, and rs619824. We further characterized the linkage disequilibrium structure for CYP17 and found that the whole CYP17 gene was located in a single-linkage disequilibrium block. This block was observed to be significantly associated with BMI. A major haplotype in this block was significantly associated with both BMI and fat mass. In conclusion, we suggest that the CYP17 gene has an effect on obesity in the Caucasian population. Further independent studies will be needed to confirm our findings.     

Genomic organization and sequence variation of the human integrin subunit alpha8 gene (ITGA8).

The integrin alpha8 is highly expressed during kidney and lung development. alpha8-deficient mice display abnormal renal development suggesting that alpha8 plays a critical role in organogenesis. Therefore, it would be of considerable interest to understand the genomic structure, localization and sequence variation of the alpha8 gene. Using FISH and genomic database analysis, we show that alpha8 gene maps to chromosome 10p13 and consists of >200 kbp organized into 30 exons. Examination of 47 individuals from two different ethnic groups (European and African descent) identified 286 varying sites. The diversity of alpha8 is comparable to that of other regions within the human genome. Eight of the varying sites were located in the coding regions: six resulted in nonsynonymous substitutions of which two lead to non-conservative changes in protein. None of the sites showed significant deviation from Hardy-Weinberg equilibrium. We mapped the coding region single nucleotide polymorphisms (SNPs) onto a model of the predicted alpha8 structure and found all the SNPs were located in the "calf" of the extracellular domain. In the European population, the linkage disequilibrium statistic D' showed three blocks of relatively non-recombinant regions in the alpha8 gene while the African population showed more evidence of recombination. The observed patterns of the linkage disequilibrium statistic R2 suggest that a large number of sites will need to be genotyped to ensure coverage of the entire gene for genetic association studies. Identification of the sequence variation will allow genetic association studies of alpha8 in kidney and lung disease.     

Genetic Variation in TLR Genes in Ugandan and South African Populations and Comparison with HapMap Data

Genetic epidemiological studies of complex diseases often rely on data from the International HapMap Consortium for identification of single nucleotide polymorphisms (SNPs), particularly those that tag haplotypes. However, little is known about the relevance of the African populations used to collect HapMap data for study populations conducted elsewhere in Africa. Toll-like receptor (TLR) genes play a key role in susceptibility to various infectious diseases, including tuberculosis. We conducted full-exon sequencing in samples obtained from Uganda (n = 48) and South Africa (n = 48), in four genes in the TLR pathway: TLR2, TLR4, TLR6, and TIRAP. We identified one novel TIRAP SNP (with minor allele frequency [MAF] 3.2%) and a novel TLR6 SNP (MAF 8%) in the Ugandan population, and a TLR6 SNP that is unique to the South African population (MAF 14%). These SNPs were also not present in the 1000 Genomes data. Genotype and haplotype frequencies and linkage disequilibrium patterns in Uganda and South Africa were similar to African populations in the HapMap datasets. Multidimensional scaling analysis of polymorphisms in all four genes suggested broad overlap of all of the examined African populations. Based on these data, we propose that there is enough similarity among African populations represented in the HapMap database to justify initial SNP selection for genetic epidemiological studies in Uganda and South Africa. We also discovered three novel polymorphisms that appear to be population-specific and would only be detected by sequencing efforts.     

Out of Africa and back again: nested cladistic analysis of human Y chromosome variation

We surveyed nine diallelic polymorphic sites on the Y chromosomes of 1,544 individuals from Africa, Asia, Europe, Oceania, and the New World. Phylogenetic analyses of these nine sites resulted in a tree for 10 distinct Y haplotypes with a coalescence time of approximately 150,000 years. The 10 haplotypes were unevenly distributed among human populations: 5 were restricted to a particular continent, 2 were shared between Africa and Europe, 1 was present only in the Old World, and 2 were found in all geographic regions surveyed. The ancestral haplotype was limited to African populations. Random permutation procedures revealed statistically significant patterns of geographical structuring of this paternal genetic variation. The results of a nested cladistic analysis indicated that these geographical associations arose through a combination of processes, including restricted, recurrent gene flow (isolation by distance) and range expansions. We inferred that one of the oldest events in the nested cladistic analysis was a range expansion out of Africa which resulted in the complete replacement of Y chromosomes throughout the Old World, a finding consistent with many versions of the Out of Africa Replacement Model. A second and more recent range expansion brought Asian Y chromosomes back to Africa without replacing the indigenous African male gene pool. Thus, the previously observed high levels of Y chromosomal genetic diversity in Africa may be due in part to bidirectional population movements. Finally, a comparison of our results with those from nested cladistic analyses of human mtDNA and beta-globin data revealed different patterns of inferences for males and females concerning the relative roles of population history (range expansions) and population structure (recurrent gene flow), thereby adding a new sex-specific component to models of human evolution.      

SNPs across time and space: population genomic signatures of founder events and epizootics in the House Finch (Haemorhous mexicanus)

Identifying genomic signatures of natural selection can be challenging against a background of demographic changes such as bottlenecks and population expansions. Here, we disentangle the effects of demography from selection in the House Finch (Haemorhous mexicanus) using samples collected before and after a pathogen‐induced selection event. Using ddRADseq, we genotyped over 18,000 SNPs across the genome in native pre‐epizootic western US birds, introduced birds from Hawaii and the eastern United States, post‐epizootic eastern birds, and western birds sampled across a similar time span. We found 14% and 7% reductions in nucleotide diversity, respectively, in Hawaiian and pre‐epizootic eastern birds relative to pre‐epizootic western birds, as well as elevated levels of linkage disequilibrium and other signatures of founder events. Despite finding numerous significant frequency shifts (outlier loci) between pre‐epizootic native and introduced populations, we found no signal of reduced genetic diversity, elevated linkage disequilibrium, or outlier loci as a result of the epizootic. Simulations demonstrate that the proportion of outliers associated with founder events could be explained by genetic drift. This rare view of genetic evolution across time in an invasive species provides direct evidence that demographic shifts like founder events have genetic consequences more widespread across the genome than natural selection.     

Haplotypes and linkage disequilibrium at the phenylalanine hydroxylase locus, PAH, in a global representation of populations

Because defects in the phenylalanine hydroxylase gene (PAH) cause phenylketonuria (PKU), PAH was studied for normal polymorphisms and linkage disequilibrium soon after the gene was cloned. Studies in the 1980s concentrated on European populations in which PKU was common and showed that haplotype-frequency variation exists between some regions of the world. In European populations, linkage disequilibrium generally was found not to exist between RFLPs at opposite ends of the gene but was found to exist among the RFLPs clustered at each end. We have now undertaken the first global survey of normal variation and disequilibrium across the PAH gene. Four well-mapped single-nucleotide polymorphisms (SNPs) spanning approximately 75 kb, two near each end of the gene, were selected to allow linkage disequilibrium across most of the gene to be examined. These SNPs were studied as PCR-RFLP markers in samples of, on average, 50 individuals for each of 29 populations, including, for the first time, multiple populations from Africa and from the Americas. All four sites are polymorphic in all 29 populations. Although all but 5 of the 16 possible haplotypes reach frequencies >5% somewhere in the world, no haplotype was seen in all populations. Overall linkage disequilibrium is highly significant in all populations, but disequilibrium between the opposite ends is significant only in Native American populations and in one African population. This study demonstrates that the physical extent of linkage disequilibrium can differ substantially among populations from different regions of the world, because of both ancient genetic drift in the ancestor common to a large regional group of modern populations and recent genetic drift affecting individual populations.     

Global Analysis of ATM Polymorphism Reveals Significant Functional Constraint

ATM, the gene that is mutated in ataxia-telangiectasia, is associated with cerebellar degeneration, abnormal proliferation of small blood vessels, and cancer. These clinically important manifestations have stimulated interest in defining the sequence variation in the ATM gene. Therefore, we undertook a comprehensive survey of sequence variation in ATM in diverse human populations. The protein-encoding exons of the gene (9,168 bp) and the adjacent intron and untranslated sequences (14,661 bp) were analyzed in 93 individuals from seven major human populations. In addition, the coding sequence was analyzed in one chimpanzee, one gorilla, one orangutan, and one Old World monkey. In human ATM, 88 variant sites were discovered by denaturing high-performance liquid chromatography, which is 96%-100% sensitive for detection of DNA sequence variation. ATM was compared to 14 other autosomal genes for nucleotide diversity. The noncoding regions of ATM had diversity values comparable to other genes, but the coding regions had very low diversity, especially in the last 29% of the protein sequence. A test of the neutral evolution hypothesis, through use of the Hudson/Kreitman/Aguadé statistic, revealed that this region of the human ATM gene was significantly constrained relative to that of the orangutan, the Old World monkey, and the mouse, but not relative to that of the chimpanzee or the gorilla. ATM displayed extensive linkage disequilibrium, consistent with suppression of meiotic recombination at this locus. Seven haplotypes were defined. Two haplotypes accounted for 82% of all chromosomes analyzed in all major populations; two others carrying the same D126E missense polymorphism accounted for 33% of chromosomes in Africa but were never observed outside of Africa. The high frequency of this polymorphism may be due either to a population expansion within Africa or to selective pressure.     

Extent of linkage disequilibrium between the androgen receptor gene CAG and GGC repeats in human populations: implications for prostate cancer risk

While studies have implicated alleles at the CAG and GGC trinucleotide repeats of the androgen receptor gene with high-grade, aggressive prostate cancer disease, little is known about the normal range of variation for these two loci, which are separated by about 1.1 kb. More importantly, few data exist on the extent of linkage disequilibrium (LD) between the two loci in different human populations. Here we present data on CAG and GGC allelic variation and LD in six diverse populations. Alleles at the CAG and GGC repeat loci of the androgen receptor were typed in over 1000 chromosomes from Africa, Asia, and North America. Levels of linkage disequilibrium between the two loci were compared between populations. Haplotype variation and diversity were estimated for each population. Our results reveal that populations of African descent possess significantly shorter alleles for the two loci than non-African populations (P<0.0001). Allelic diversity for both markers was higher among African Americans than any other population, including indigenous Africans from Sierra Leone and Nigeria. Analysis of molecular variance revealed that approx. 20% of CAG and GGC repeat variance could be attributed to differences between the populations. All non-African populations possessed the same common haplotype while the three populations of African descent possessed three divergent common haplotypes. Significant LD was observed in our sample of healthy African Americans. The LD observed in the African American population may be due to several reasons; recent migration of African Americans from diverse rural communities following urbanization, recurrent gene flow from diverse West African populations, and admixture with European Americans. This study represents the largest genotyping effort to be performed on the two androgen receptor trinucleotide repeat loci in diverse human populations.     

The extent of linkage disequilibrium caused by selection on G6PD in humans

The gene coding for glucose-6-phosphate dehydrogenase (G6PD) is subject to positive selection by malaria in some human populations. The G6PD A- allele, which is common in sub-Saharan Africa, is associated with deficient enzyme activity and protection from severe malaria. To delimit the impact of selection on patterns of linkage disequilibrium (LD) and nucleotide diversity, we resequenced 5.1 kb at G6PD and approximately 2-3 kb at each of eight loci in a 2.5-Mb region roughly centered on G6PD in a diverse sub-Saharan African panel of 51 unrelated men (including 20 G6PD A-, 11 G6PD A+, and 20 G6PD B chromosomes). The signature of selection is evident in the absence of genetic variation at G6PD and at three neighboring loci within 0.9 Mb from G6PD among all individuals bearing G6PD A- alleles. A genomic region of approximately 1.6 Mb around G6PD was characterized by long-range LD associated with the A- alleles. These patterns of nucleotide variability and LD suggest that G6PD A- is younger than previous age estimates and has increased in frequency in sub-Saharan Africa due to strong selection (0.1 < s < 0.2). These results also show that selection can lead to nonrandom associations among SNPs over great physical and genetic distances, even in African populations.     

Using whole-genome sequence data to predict quantitative trait phenotypes in Drosophila melanogaster

Predicting organismal phenotypes from genotype data is important for plant and animal breeding, medicine, and evolutionary biology. Genomic-based phenotype prediction has been applied for single-nucleotide polymorphism (SNP) genotyping platforms, but not using complete genome sequences. Here, we report genomic prediction for starvation stress resistance and startle response in Drosophila melanogaster, using ～2.5 million SNPs determined by sequencing the Drosophila Genetic Reference Panel population of inbred lines. We constructed a genomic relationship matrix from the SNP data and used it in a genomic best linear unbiased prediction (GBLUP) model. We assessed predictive ability as the correlation between predicted genetic values and observed phenotypes by cross-validation, and found a predictive ability of 0.239±0.008 (0.230±0.012) for starvation resistance (startle response). The predictive ability of BayesB, a Bayesian method with internal SNP selection, was not greater than GBLUP. Selection of the 5% SNPs with either the highest absolute effect or variance explained did not improve predictive ability. Predictive ability decreased only when fewer than 150,000 SNPs were used to construct the genomic relationship matrix. We hypothesize that predictive power in this population stems from the SNP-based modeling of the subtle relationship structure caused by long-range linkage disequilibrium and not from population structure or SNPs in linkage disequilibrium with causal variants. We discuss the implications of these results for genomic prediction in other organisms.     

Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms in five populations

Single-nucleotide polymorphisms (SNPs) may be extremely important for deciphering the impact of genetic variation on complex human diseases. The ultimate value of SNPs for linkage and association mapping studies depends in part on the distribution of SNP allele frequencies and intermarker linkage disequilibrium (LD) across populations. Limited information is available about these distributions on a genomewide scale, particularly for LD. Using 114 SNPs from 33 genes, we compared these distributions in five American populations (727 individuals) of African, European, Chinese, Hispanic, and Japanese descent. The allele frequencies were highly correlated across populations but differed by >20% for at least one pair of populations in 35% of SNPs. The correlation in LD was high for some pairs of populations but not for others (e.g., Chinese American or Japanese American vs. any other population). Regardless of population, average minor-allele frequencies were significantly higher for SNPs in noncoding regions (20%-25%) than for SNPs in coding regions (12%-16%). Interestingly, we found that intermarker LD may be strongest with pairs of SNPs in which both markers are nonconservative substitutions, compared to pairs of SNPs where at least one marker is a conservative substitution. These results suggest that population differences and marker location within the gene may be important factors in the selection of SNPs for use in the study of complex disease with linkage or association mapping methods.     

Population Structure and the Colonization Route of One of the Oldest North American Invasive Insects: Stories from the Worn Road of the Hessian Fly,Mayetiola destructor (Say)

An integral part to understanding the biology of an invasive species is determining its origin, particularly in pest species. As one of the oldest known invasive species, the goals of this study were to evaluate the evidence of a westward expansion of Hessian fly into North America, from a potential singular introduction event, and the population genetic structure of current populations. Levels of genetic diversity and population structure in the Hessian fly were compared across North America, Europe, North Africa, Western Asia, and New Zealand. Furthermore, Old World populations were evaluated as possible sources of introduction. We tested diversity and population structure by examining 18 microsatellite loci with coverage across all four Hessian fly chromosomes. Neither genetic diversity nor population genetic structure provided evidence of a westward movement from a single introduction in North America. Introduced populations in North America did not show identity or assignment to any Old World population, likely indicating a multiple introduction scenario with subsequent gene flow between populations. Diversity and selection were assessed on a chromosomal level, with no differences in diversity or selection between chromosomes or between native and introduced populations.