Explore Population Structure Research Articles

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant.

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Genetic diversity and wing geometric morphometrics among four populations of Aedes aegypti (Diptera: Culicidae) from Benin

BackgroundThe impact of the arbovirus vector Aedes aegypti is of major concern for global public health as the viruses that it transmits affect millions of people each year worldwide. Originating in Africa, Ae. aegypti has now spread throughout much of the world. While the genetic makeup of Ae. aegypti in the New World has been extensively studied, there is limited knowledge on its genetic diversity in Africa, particularly at a microgeographical level.MethodsWe investigated mitochondrial cytochrome oxidase I of four Ae. aegypti populations from Benin and employed wing morphometric analyses as a cost-effective and reliable tool to explore population structure. Our sampling encompassed various areas of Benin, from the southern to the northern borders of the country, and included urban, semi-urban, and sylvatic sites.ResultsWe observed a notable level of genetic diversity (haplotype diversity of 0.8333) and nucleotide diversity (0.00421986), and identified seven distinct haplotypes. Sylvatic and semi-urban sites exhibited a greater number of haplotypes compared to urban sites. Utilizing 18 wing landmarks, we calculated the centroid size, which revealed significant variation among the three landscape types. However, principal component analysis, employed to assess wing shape variation, did not demonstrate significant differences between populations based on landscape type.ConclusionsOur findings indicate substantial genetic and morphological diversity among Ae. aegypti populations in Benin, and provide insight into important biological characteristics of these populations with respect to their potential to transmit viruses. To the best of our knowledge, this is the first study undertaken in Africa to integrate genetics with morphology to analyse the population structure of the major arbovirus vector Ae. aegypti.Graphical

Marker-trait association analyses revealed major novel QTLs for grain yield and related traits in durum wheat.

The growing global demand for wheat for food is rising due to the influence of population growth and climate change. The dissection of complex traits by employing a genome-wide association study (GWAS) allows the identification of DNA markers associated with complex traits to improve the productivity of crops. We used GWAS with 10,045 single nucleotide polymorphism (SNP) markers to search for genomic regions associated with grain yield and related traits based on diverse panels of Ethiopian durum wheat. In Ethiopia, multi-environment trials of the genotypes were carried out at five locations. The genotyping was conducted using the 25k Illumina Wheat SNP array to explore population structure, linkage disequilibrium (LD), and marker-trait associations (MTAs). For GWAS, the multi-locus Fixed and Random Model Circulating Probability Unification (FarmCPU) model was applied. Broad-sense heritability estimates were high, ranging from 0.63 (for grain yield) to 0.97 (for thousand-kernel weight). The population structure based on principal component analysis, and model-based cluster analysis revealed two genetically distinct clusters with limited admixtures. The LD among SNPs declined within the range of 2.02-10.04 Mbp with an average of 4.28 Mbp. The GWAS scan based on the mean performance of the genotypes across the environments identified 44 significant MTAs across the chromosomes. Twenty-six of these MTAs are novel, whereas the remaining 18 were previously reported and confirmed in this study. We also identified candidate genes for the novel loci potentially regulating the traits. Hence, this study highlights the significance of the Ethiopian durum wheat gene pool for improving durum wheat globally. Furthermore, a breeding strategy focusing on accumulating favorable alleles at these loci could improve durum wheat production in the East African highlands and elsewhere.

Diversification and historical demography of Rhampholeon spectrum in West-Central Africa.

Pygmy Chameleons of the genus Rhampholeon represent a moderately diverse, geographically circumscribed radiation, with most species (18 out of 19 extant taxa) limited to East Africa. The one exception is Rhampholeon spectrum, a species restricted to West-Central African rainforests. We set out to characterize the geographic basis of genetic variation in this disjunctly distributed Rhampholeon species using a combination of multilocus Sanger data and genomic sequences to explore population structure and range-wide phylogeographic patterns. We also employed demographic analyses and niche modeling to distinguish between alternate explanations to contextualize the impact of past geological and climatic events on the present-day distribution of intraspecific genetic variation. Phylogenetic analyses suggest that R. spectrum is a complex of five geographically delimited populations grouped into two major clades (montane vs. lowland). We found pronounced population structure suggesting that divergence and, potentially, speciation began between the late Miocene and the Pleistocene. Sea level changes during the Pleistocene climatic oscillations resulted in allopatric divergence associated with dispersal over an ocean channel barrier and colonization of Bioko Island. Demographic inferences and range stability mapping each support diversification models with secondary contact due to population contraction in lowland and montane refugia during the interglacial period. Allopatric divergence, congruent with isolation caused by geologic uplift of the East African rift system, the "descent into the Icehouse," and aridification of sub-Saharan Africa during the Eocene-Oligocene are identified as the key events explaining the population divergence between R. spectrum and its closely related sister clade from the Eastern Arc Mountains. Our results unveil cryptic genetic diversity in R. spectrum, suggesting the possibility of a species complex distributed across the Lower Guinean Forest and the Island of Bioko. We highlight the major element of species diversification that modelled today's diversity and distributions in most West-Central African vertebrates.

Forest cover and geographical distance influence fine-scale genetic structure of leaf-toed geckos in the tropical dry forests of western Mexico

Abstract The biodiversity within tropical dry forests (TDFs) is astounding and yet poorly catalogued due to inadequate sampling and the presence of cryptic species. In the Mexican TDF, endemic species are common, and the landscape has been continually altered by geological and anthropogenic changes. To understand how landscape and environmental variables have shaped the population structure of endemic species, we studied the recently described species of leaf-toed gecko, Phyllodactylus benedettii, in coastal western Mexico. Using double-digest restriction site-associated DNA sequencing data, we first explore population structure and estimate the number of ancestral populations. The results indicate a high degree of genetic structure with little admixture, and patterns corresponding to both latitudinal and altitudinal gradients. We find that genetic structure cannot be explained purely by geographical distance, and that ecological corridors may facilitate dispersal and gene flow. We then model the spatial distribution of P. benedettii in the TDF through time and find that the coastline has been climatically suitable for the species since the Last Glacial Maximum. Landscape genetic analyses suggest that the combination of isolation by distance (IBD) and isolation by resistance (IBR; forest cover) has influenced the spatial genetic structure of the species. Overall, our genomic data demonstrate fine-scale population structure in TDF habitat, a complex colonization history, and spatial patterns consistent with both IBD and other ecological factors. These results further highlight the Mexican TDF as a diversity hotspot and suggest that continued anthropogenic changes are likely to affect native fauna.

Population structure of Caspian Kutum (Rutilus frisii, Nordmann, 1840) in the southern coast of Caspian Sea using Genome-wide single nucleotide polymorphism markers

The Caspian kutum (Rutilus frisii, Nordmann, 1840) is one of the most important semi-anadromous endemic fish species in the southern part of the Caspian Sea. The fisheries from wild populations of R. frisii are declined drastically due to variability of sea levels, overfishing, and illegal fishing. For this regard, genetic studies is very noticed to maintain of remained this endangered species. Here, we used restriction site-associated DNA sequencing (RAD-seq) approach to explore population structure, genetic variability, and single nucleotide polymorphisms (SNPs) from 146 individuals from eight rivers in Iran. In total 86,274 SNP loci were determined, and after filtering, 16,193 high-quality SNP markers were applied in statistical analyses. Twenty-eight percent of selected SNPs were located in genic regions and about 65% of these SNPs were in the coding sequence. Of the 102 SNP loci, 42 loci were found to be polymorphic and bi-allelic. The minor allele frequency (MAF) of these SNPs varied from 0.156 to 0.471. The observed heterozygosity (Ho) per locus ranged from 0.119 to 0.864, while the expected heterozygosity (He) per locus was from 0.142 to 0.569. Among these SNP loci, 22 SNP loci were significantly deviated from the Hardy–Weinberg equilibrium (P < 0.01) after Bonferroni correction. Both Principal coordinate analysis (PCA) and Bayesian clustering satisfyingly showed population genomic structure in Caspian kutum populations, indicating the description of distinct population segments. The first development of SNP markers in this study provides the better monitoring of population genetic variation in Caspian kutum at a more detailed level and in the scarce of genetic information about this species.

A multiple phenotype imputation method for genetic diversity and core collection in Taiwanese vegetable soybean.

Establishment of vegetable soybean (edamame) [Glycine max (L.) Merr.] germplasms has been highly valued in Asia and the United States owing to the increasing market demand for edamame. The idea of core collection (CC) is to shorten the breeding program so as to improve the availability of germplasm resources. However, multidimensional phenotypes typically are highly correlated and have different levels of missing rate, often failing to capture the underlying pattern of germplasms and select CC precisely. These are commonly observed on correlated samples. To overcome such scenario, we introduced the “multiple imputation” (MI) method to iteratively impute missing phenotypes for 46 morphological traits and jointly analyzed high-dimensional imputed missing phenotypes (ECimpu) to explore population structure and relatedness among 200 Taiwanese vegetable soybean accessions. An advanced maximization strategy with a heuristic algorithm and PowerCore was used to evaluate the morphological diversity among the ECimpu. In total, 36 accessions (denoted as CCimpu) were efficiently selected representing high diversity and the entire coverage of the ECimpu. Only 4 (8.7%) traits showed slightly significant differences between the CCimpu and ECimpu. Compared to the ECimpu, 96% traits retained all characteristics or had a slight diversity loss in the CCimpu. The CCimpu exhibited a small percentage of significant mean difference (4.51%), and large coincidence rate (98.1%), variable rate (138.76%), and coverage (close to 100%), indicating the representativeness of the ECimpu. We noted that the CCimpu outperformed the CCraw in evaluation properties, suggesting that the multiple phenotype imputation method has the potential to deal with missing phenotypes in correlated samples efficiently and reliably without re-phenotyping accessions. Our results illustrated a significant role of imputed missing phenotypes in support of the MI-based framework for plant-breeding programs.

Incomplete lineage sorting and ancient admixture, and speciation without morphological change in ghost-worm cryptic species.

Morphologically similar species, that is cryptic species, may be similar or quasi-similar owing to the deceleration of morphological evolution and stasis. While the factors underlying the deceleration of morphological evolution or stasis in cryptic species remain unknown, decades of research in the field of paleontology on punctuated equilibrium have originated clear hypotheses. Species are expected to remain morphologically identical in scenarios of shared genetic variation, such as hybridization and incomplete lineage sorting, or in scenarios where bottlenecks reduce genetic variation and constrain the evolution of morphology. Here, focusing on three morphologically similar Stygocapitella species, we employ a whole-genome amplification method (WGA) coupled with double-digestion restriction-site associated DNA sequencing (ddRAD) to reconstruct the evolutionary history of the species complex. We explore population structure, use population-level statistics to determine the degree of connectivity between populations and species, and determine the most likely demographic scenarios which generally reject for recent hybridization. We find that the combination of WGA and ddRAD allowed us to obtain genomic-level data from microscopic eukaryotes (∼1 millimetre) opening up opportunities for those working with population genomics and phylogenomics in such taxa. The three species share genetic variance, likely from incomplete lineage sorting and ancient admixture. We speculate that the degree of shared variation might underlie morphological similarity in the Atlantic species complex.

Population epigenetic divergence exceeds genetic divergence in the Eastern oyster Crassostrea virginica in the Northern Gulf of Mexico.

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome‐wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site‐associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene‐wide averages of F ST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low F ST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in F ST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population‐specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene‐body mean F ST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

Exploring Population Structure with Admixture Models and Principal Component Analysis.

Population structure is a commonplace feature of genetic variation data, and it has importance in numerous application areas, including evolutionary genetics, conservation genetics, and human genetics. Understanding the structure in a sample is necessary before more sophisticated analyses are undertaken. Here we provide a protocol for running principal component analysis (PCA) and admixture proportion inference-two of the most commonly used approaches in describing population structure. Along with hands-on examples with CEPH-Human Genome Diversity Panel and pragmatic caveats, readers will learn to analyze and visualize population structure on their own data.

Genome-wide SNP discovery through genotyping by sequencing, population structure, and linkage disequilibrium in Brazilian peach breeding germplasm

Genotyping by sequencing (GBS) is a flexible and cost-effective strategy for genome-wide SNP discovery and high-throughput genotyping. Here, we employed the GBS approach to explore population structure, genetic variability, and patterns of linkage disequilibrium (LD) among 220 peach genotypes (Prunus persica), representative of the Brazilian breeding germplasm. This selected panel is mainly composed by locally adapted peaches developed in mild winter and high relative humidity conditions, and represents a worldwide reference of peach germplasm for low-chilling areas. A total of 93,353 SNP markers were discovered, and after filtering, 18,373 high-quality SNPs were used in analyses. Thirty-four percent of selected SNPs were located in genic regions and about 70% of these in the coding sequence. Principal coordinate analysis (PCoA) and Bayesian clustering (fastSTRUCTURE) successfully detected population genetic structure in our germplasm panel, supporting the identification of three main distinct subpopulations. The distribution of the genotypes within subpopulations reflected their fruit-related traits: melting and non-melting flesh cultivars. LD patterns revealed a medium LD level in peach, with the extent of LD highly dependent on the subpopulation and genome regions. The use of GBS-SNPs to perform genome-wide association studies (GWAS) was also exploited by using phenotypic information from five fruit quality traits. The observed significant SNP-trait associations were always in the regions previously predicted by linkage and association mapping analysis, providing valuable insights to incorporate the use of marker-assisted selection in the breeding program.

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus).

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype‐by‐sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.

Variations in body shape of mountain habitat specialist Carabus croaticus and its sister species Carabus caelatus (Coleoptera: Carabidae) populations across Dinaric Alps

There are abundant phenotypic variations of Carabus croaticus Dejean 1826, an endemic species of the Dinaric Alps, and its sister species Carabus caelatus Fabricius 1801, resulting in taxonomic inflation at intraspecific level, synonymy and inconsistency between relevant catalogues and checklists. The main aims of this research were to explore population structure based on morphology and geographic patterns of phenotypic variability, and compare the results with taxa described by nonmetric visual comparisons of morphological traits. Our study included 224 specimens of C. croaticus and 192 specimens of C. caelatus, covering most of their distributional range. Shapes of the pronotum and head were analysed using geometric morphometrics (GMM). Principal component and canonical variate analyses were used to characterize the main features of shape variation between populations and mountain ranges. GMM delimitated interspecific morphological variations but at the intraspecific level it showed many overlaps within populations for both species. Conducting the morphological analyses for the first time on most of the described phenotypic variants of studied species, we wanted to provide a new evidence for the possible solution of the taxonomic relations within these two endemic species by measuring body shape variability, and thus to enable a better understanding of the evolutionary processes in the Dinaric Mountains.

In the footsteps of Wallace: population structure in the breadfruit fruit fly, Bactrocera umbrosa (F.) (Diptera: Tephritidae), suggests disjunction across the Indo‐Australian Archipelago

AbstractUnderstanding the interplay between plant host and insect herbivore diversification underpins many areas of pure and applied research. The tephritid fruit fly Bactrocera umbrosa is a primary pest of a small number of Artocarpus species throughout Southeast Asia and the West Pacific. Recent molecular evidence supports a pattern of eastward migration and species diversification in Artocarpus. Here, we aimed to test whether population structure in B. umbrosa was associated with historical biogeographical barriers such as Wallace's Line and discuss observed patterns in the context of Artocarpus diversification. We used an integrative approach to explore population structure within B. umbrosa based on morphological (wing shape and aedeagus length) and molecular (mitochondrial COI and COII) data. Overall, aedeagi and wing centroid sizes were generally larger and exhibited greater variation in the West Pacific than Southeast Asia. Molecular data agreed with this trend, and COI also showed a subtle but clear disjunction between regions associated with Weber/Lydekker's Lines. Taken together, the West Pacific was supported as the putative origin of B. umbrosa, whereas movement westward into Southeast Asia occurred more recently, likely via a single colonisation event followed by highly restricted gene flow. Population structure in B. umbrosa does not reflect an ancient history of tracking Artocarpus diversification eastward out of Southeast Asia.

Characterization of microsatellite markers for Moricandia moricandioides (Brassicaceae) and related species.

Premise of the StudyPolymorphic microsatellite markers were developed to study population structure and mating patterns of the monocarpic herb Moricandia moricandioides (Brassicaceae).Methods and ResultsIllumina MiSeq sequencing was used to develop a panel of 15 polymorphic microsatellite markers that were tested across 77 individuals from three populations on the Iberian Peninsula. All markers were polymorphic in at least two studied populations, and the number of alleles ranged from one to 11 per locus. The levels of observed and expected heterozygosity ranged from 0.000 to 1.000 and from 0.153 to 0.865, respectively. Nine and 11 loci were successfully amplified in the congeneric species M. arvensis and M. foetida, respectively.ConclusionsThe 15 microsatellite markers will be useful for population genetic studies of the genus Moricandia. These markers will serve as a useful tool for exploring population structure and mating patterns of M. moricandioides.

Adaptation in temporally variable environments: stickleback armor in periodically breaching bar-built estuaries.

The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar-built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re-form within a few weeks or months, thereby re-isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within-population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure and then analysed how stickleback armor traits, the associated genes Eda and Pitx1 and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among-population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within-population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype and the presence of predators, thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within-population variation while reducing, but not eliminating, among-population variation driven by spatial environmental variation.

Dynamism in the upstream invasion edge of a freshwater fish exposes range boundary constraints.

Studying the dynamics of species' borders can provide insight into the mechanisms limiting or promoting range expansion in response to environmental change. In the John Day River, Oregon (USA), rising stream temperatures are facilitating the upstream expansion of invasive smallmouth bass Micropterus dolomieu. Here, where smallmouth bass occupy the upstream limit of its thermal tolerance, we explore population structure and seasonal movement patterns to elucidate the environmental conditions and individual traits that define front edge (where individuals reside year-round) and leading edge (where individuals colonize, but may not establish) limits to its upstream distribution. Reporting on a multi-year, spatially extensive riverscape survey, our results show dramatic ebbs and flows of seasonal occupancies due to individual movement with an overall trend of upstream expansion. We revealed distinct front and leading edge invasion extents, each constrained by different ecological conditions. The front edge is largely constrained by the ability for juveniles to survive an overwinter starvation period, whereas the leading edge is associated with adult growth potential and seasonal hydrological conditions. We also found key morphological traits associated with more mobile individuals. By providing mechanistic insight into the factors that promote or limit range expansion of an invasive riverine species, our study enhances the ability to predict future range shifts and provides critical information to managers tasked with restricting further expansion.

Genetic variability and population structure of grey wolf (Canis lupus) in Serbia

Results of previous morphometric and genetic analyses of grey wolf (Canis lupus L.) population from Serbia indicated different patterns of population subdivision. In order to explore population structure, level of genetic variability, genetic drift, inbreeding and signals of bottleneck for grey wolves from Serbia, we applied highly polymorphic genetic markers (microsatellites). Obtained data are valuable in determination of conservation units and creation of appropriate management plans. We have amplified 18 highly polymorphic microsatellites, in a total sample of 75 grey wolves, from different localities across Serbia and multilocus genotypes were analyzed using appropriate software. Observed values of the basic genetic parameters (H O = 0.69; H E = 0.75) indicated moderate level of genetic variability, similar to genetic variability in other populations belonging to the Dinaric-Balkan population of grey wolf. In STRUCTURE analysis, although ΔK was estimated to be at first peak K = 2, and second peak K = 4, CLUMPAK analyses showed that there’s no structuring for any of assumed K, and therefore the population of grey wolf from Serbia may be considered as one continuous population and treated as one conservation unit in future management plans. Signals of bottleneck haven’t been observed (Wilcoxon test two phase mutation model p = 0.247; and stepwise mutation model p = 0.815).

Genotyping by sequencing resolves shallow population structure to inform conservation of Chinook salmon (Oncorhynchus tshawytscha)

Recent advances in population genomics have made it possible to detect previously unidentified structure, obtain more accurate estimates of demographic parameters, and explore adaptive divergence, potentially revolutionizing the way genetic data are used to manage wild populations. Here, we identified 10 944 single-nucleotide polymorphisms using restriction-site-associated DNA (RAD) sequencing to explore population structure, demography, and adaptive divergence in five populations of Chinook salmon (Oncorhynchus tshawytscha) from western Alaska. Patterns of population structure were similar to those of past studies, but our ability to assign individuals back to their region of origin was greatly improved (>90% accuracy for all populations). We also calculated effective size with and without removing physically linked loci identified from a linkage map, a novel method for nonmodel organisms. Estimates of effective size were generally above 1000 and were biased downward when physically linked loci were not removed. Outlier tests based on genetic differentiation identified 733 loci and three genomic regions under putative selection. These markers and genomic regions are excellent candidates for future research and can be used to create high-resolution panels for genetic monitoring and population assignment. This work demonstrates the utility of genomic data to inform conservation in highly exploited species with shallow population structure.

An algorithm for deciding the number of clusters and validation using simulated data with application to exploring crop population structure

A first step in exploring population structure in crop plants and other organisms is to define the number of subpopulations that exist for a given data set. The genetic marker data sets being generated have become increasingly large over time and commonly are of the high-dimension, low sample size (HDLSS) situation. An algorithm for deciding the number of clusters is proposed, and is validated on simulated data sets varying in both the level of structure and the number of clusters covering the range of variation observed empirically. The algorithm was then tested on six empirical data sets across three small grain species. The algorithm uses bootstrapping, three methods of clustering, and defines the optimum number of clusters based on a common criterion, the Hubert's gamma statistic. Validation on simulated sets coupled with testing on empirical sets suggests that the algorithm can be used for a wide variety of genetic data sets.