Hierarchical Population Structure Research Articles

Living in a post‐industrial landscape: repeated patterns of genetic divergence in brown trout (Salmo trutta L.) across the British Isles

AbstractAimThe British Isles have been worked for millennia to extract metal ores to feed industrial development, leaving a legacy of mine water pollution that continues to impact freshwater communities in many regions. Brown trout (Salmo trutta L.) have long been observed to persist in these metal‐impacted systems as apex predators, with previous studies showing a small number of impacted populations to be highly genetically divergent. We sought to understand the scale of genetic diversity across regions and the repeatability of genetic divergence in trout populations affected by metal pollution.LocationWe examined four mine water‐impacted regions across the British Isles: west Wales, northeast England, southwest England and southeast Ireland.MethodsWe employed a panel of 95 SNP loci to screen 1236 individuals from 71 sites representing paired metal‐impacted and clean sites from across the four regions. From these, we obtained diversity statistics, assessed genetic structuring of populations and modelled historical demographic scenarios to understand which factors most credibly explain genetic variation in divergent populations.ResultsWe evidenced hierarchical population structure in the regions studied, in line with expectations from phylogeographic history. However, in a hierarchical analysis of genetic structuring the first level of differentiation was driven by the divergence of the metal‐impacted trout of Cornwall in southwest England. Within regions we observed reduced genetic diversity and repeated patterns of local genetic sub‐structuring between paired samples from metal‐impacted and relatively clean sites. Demographic history analyses suggested the timing of these splits to be relatively recent and to be associated with periods of peak mining activity.Main conclusionsOur findings demonstrate distinct patterns of genetic isolation and reduced diversity arising from legacy pollution in freshwater ecosystems, with impacts being most apparent where both chemical pollution and physical barriers are present. Management should focus on the amelioration of mine water wash‐out and the removal of barriers to fish movement to safeguard genetic diversity in impacted populations.

A Comparison of Entropic Diversity and Variance in the Study of Population Structure.

AMOVA is a widely used approach that focuses on variance within and among strata to study the hierarchical genetic structure of populations. The recently developed Shannon Informational Diversity Translation Analysis (SIDTA) instead tackles exploration of hierarchical genetic structure using entropic allelic diversity. A mix of artificial and natural population data sets (including allopolyploids) is used to compare the performance of SIDTA (a 'q = 1' diversity measure) vs. AMOVA (a 'q = 2' measure) under different conditions. An additive allelic differentiation index based on entropic allelic diversity measuring the mean difference among populations (ΩAP) was developed to facilitate the comparison of SIDTA with AMOVA. These analyses show that the genetic population structure seen by AMOVA is notably different in many ways from that provided by SIDTA, and the extent of this difference is greatly affected by the stability of the markers employed. Negative among group values are lacking with SIDTA but occur with AMOVA, especially with allopolyploids. To provide more focus on measuring allelic differentiation among populations, additional measures were also tested including Bray-Curtis Genetic Differentiation (BCGD) and several expected heterozygosity-based indices (e.g., GST, G″ST, Jost's D, and DEST). Corrections, such as almost unbiased estimators, that were designed to work with heterozygosity-based fixation indices (e.g., FST, GST) are problematic when applied to differentiation indices (eg., DEST, G″ST, G'STH).

Identifying New Clusterons: Application of TBEV Analyzer 3.0.

Early knowledge about novel emerging viruses and rapid determination of their characteristics are crucial for public health. In this context, development of theoretical approaches to model viral evolution are important. The clusteron approach is a recent bioinformatics tool which analyzes genetic patterns of a specific E protein fragment and provides a hierarchical network structure of the viral population at three levels: subtype, lineage, and clusteron. A clusteron is a group of strains with identical amino acid (E protein fragment) signatures; members are phylogenetically closely related and feature a particular territorial distribution. This paper announces TBEV Analyzer 3.0, an analytical platform for rapidly characterizing tick-borne encephalitis virus (TBEV) strains based on the clusteron approach, workflow optimizations, and simplified parameter settings. Compared with earlier versions of TBEV Analyzer, we provide theoretical and practical enhancements to the platform. Regarding the theoretical aspect, the model of the clusteron structure, which is the core of platform analysis, has been updated by analyzing all suitable TBEV strains available in GenBank, while the practical enhancements aim at improving the platform's functionality. Here, in addition to expanding the strain sets of prior clusterons, we introduce eleven novel clusterons through our experimental results, predominantly of the European subtype. The obtained results suggest effective application of the proposed platform as an analytical and exploratory tool in TBEV surveillance.

AdmixPipe v3: facilitating population structure delimitation from SNP data.

Quantifying genetic clusters (=populations) from genotypic data is a fundamental, but non-trivial task for population geneticists that is compounded by: hierarchical population structure, diverse analytical methods, and complex software dependencies. AdmixPipe v3 ameliorates many of these issues in a single bioinformatic pipeline that facilitates all facets of population structure analysis by integrating outputs generated by several popular packages (i.e. CLUMPAK, EvalAdmix). The pipeline interfaces disparate software packages to parse Admixture outputs and conduct EvalAdmix analyses in the context of multimodal population structure results identified by CLUMPAK. We further streamline these tasks by packaging AdmixPipe v3 within a Docker container to create a standardized analytical environment that allows for complex analyses to be replicated by different researchers. This also grants operating system flexibility and mitigates complex software dependencies. Source code, documentation, example files, and usage examples are freely available at https://github.com/stevemussmann/admixturePipeline. Installation is facilitated via Docker container available from https://hub.docker.com/r/mussmann/admixpipe. Usage under Windows operating systems requires the Windows Subsystem for Linux.

Integrative characterization of genetic and phenotypic differentiation in an ant species complex with strong hierarchical population structure and low dispersal abilities.

Low dispersal, occurrence of asexual reproduction and geographic discontinuity increase genetic differentiation between populations, which ultimately can lead to speciation. In this work, we used a multidisciplinary framework to characterize the genetic and phenotypic differentiation between and within two cryptic ant species with restricted dispersal, Cataglyphis cursor and C. piliscapa and used behavioral experiments to test for reproductive isolation. Their distribution is segregated by the Rhône River and they have been traditionally distinguished only by hair numbers, although a statistical assessment is still lacking. We found strong genetic (microsatellites, nuclear and mitochondrial sequences), morphological (number of hairs, tibia length, male genitalia) and chemical (cuticular hydrocarbons) differentiation not only between species but also among localities within species. However, inter-specific differentiation was slightly higher than intra-specific differentiation for most markers. Overall, this pattern could either reflect reproductive isolation or could result from a longer period of geographic isolation between species than among localities within species without necessarily involving reproductive isolation. Interestingly, our behavioral experiments showed an absence of mating between species associated to a higher aggressiveness of workers towards heterospecific males. This suggests that sexual selection may, at least partially, fuel reproductive isolation. We also showed that cuticular hydrocarbons, mtDNA sequences and number of hairs provide reliable criteria allowing species discrimination. Overall, this species complex offers a case study to further investigate varying stages of a speciation continuum by estimating reproductive isolation between pairs of localities varying by their level of genetic differentiation.

Comparison of Tetraploid Alfalfa (Medicago sativa L.) Populations Collected from Turkey and Former Soviet Countries

Tetraploid Medicago sativa L. subspecies; M. sativa subsp. sativa, M. sativa subsp. falcata and M. sativa subsp. varia form the primary gene pool of the alfalfa. The center of diversity for this subspecies is seen as the Caucasus, Northwest Iran and Northeast Turkey, and its natural range is the former Soviet Union and the southern border of North Africa as the northern border. Genetic diversity among the primary gene pool and comparisons between regions of diversity provides a good reference for breeders when utilizing genetic resources. The United States Department of Agriculture Genetic Resources Information Network (USDA-GRIN) System provides reference data from its entire natural range, including Turkey and the Former Soviet region. In this study, seven populations collected from Turkey and thirteen alfalfa populations collected from Former Soviet Countries held in the USDA GRIN System were evaluated using 20 SSR markers. Within the scope of the study, the information between the locations of the subspecies was compared to reveal the hierarchical population structure. The results obtained from the STRUCTURE and PCA analyzes show that the populations are clustered in two main groups for both countries, but there is a high similarity in hybrid genome contents in the subspecies belonging to the Former Soviet countries. In addition, it was determined by AMOVA analysis that the variance within the populations was higher than that between the populations according to the subspecies analyzed from both countries. It is thought that the results will be effective in terms of using alfalfa genetic resources of these countries in breeding programs.

A multi-objective chicken swarm optimization algorithm based on dual external archive with various elites

Multi-objective optimization problems (MOPs) that widely exist in real world concern all optimal solutions compromised among multiple objectives. Chicken swarm optimization algorithm derived from emergent behaviors of organisms provides an effective way for handling MOPs. To speed up convergence and improve uniformity of Pareto-optimal solutions, a multi-objective chicken swarm optimization algorithm based on dual external archives and boundary learning strategy (MOCSO-DABL) is proposed in this paper. Dual external archives are employed to distinguish and choose two types of elite solutions, with the purpose of more effectively guiding individual evolution. A boundary learning strategy guides the chickens to learn from boundary individuals in the later stage of evolution. Moreover, fast non-dominated sorting is adopted to establish the hierarchical social structure of a chicken population, and learning strategies of roosters, hens and chicks are improved to meet the requirements of MOPs. Experimental results on 14 benchmark functions show that the proposed MOCSO-DABL outperforms other five state-of-the-art algorithms significantly.

Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads.

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units.

Independent origins of fetal liver haematopoietic stem and progenitor cells.

Self-renewal and differentiation are tightly controlled to maintain haematopoietic stem cell (HSC) homeostasis in the adult bone marrow1,2. During fetal development, expansion of HSCs (self-renewal) and production of differentiated haematopoietic cells (differentiation) are both required to sustain the haematopoietic system for body growth3,4. However, it remains unclear how these two seemingly opposing tasks are accomplished within the short embryonic period. Here we used in vivo genetic tracing in mice to analyse the formation of HSCs and progenitors from intra-arterial haematopoietic clusters, which contain HSC precursors and express the transcription factor hepatic leukaemia factor (HLF). Through kinetic study, we observed the simultaneous formation of HSCs and defined progenitors-previously regarded as descendants of HSCs5-from the HLF+ precursor population, followed by prompt formation of the hierarchical haematopoietic population structure in the fetal liver in an HSC-independent manner. The transcription factor EVI1 is heterogeneously expressed within the precursor population, with EVI1hi cells being predominantly localized to intra-embryonic arteries and preferentially giving rise to HSCs. By genetically manipulating EVI1 expression, we were able to alter HSC and progenitor output from precursors in vivo. Using fate tracking, we also demonstrated that fetal HSCs are slowly used to produce short-term HSCs at late gestation. These data suggest that fetal HSCs minimally contribute to the generation of progenitors and functional blood cells before birth. Stem cell-independent pathways during development thus offer a rational strategy for the rapid and simultaneous growth of tissues and stem cell pools.

Accounting for population structure in genomic predictions of Eucalyptus globulus.

Genetic groups have been widely adopted in tree breeding to account for provenance effects within pedigree-derived relationship matrices. However, provenances or genetic groups have not yet been incorporated into single-step genomic BLUP ("HBLUP") analyses of tree populations. To quantify the impact of accounting for population structure in Eucalyptus globulus, we used HBLUP to compare breeding value predictions from models excluding base population effects and models including either fixed genetic groups or the marker-derived proxies, also known as metafounders. Full-sib families from 2 separate breeding populations were evaluated across 13 sites in the "Green Triangle" region of Australia. Gamma matrices (Γ) describing similarities among metafounders reflected the geographic distribution of populations and the origins of 2 land races were identified. Diagonal elements of Γ provided population diversity or allelic covariation estimates between 0.24 and 0.56. Genetic group solutions were strongly correlated with metafounder solutions across models and metafounder effects influenced the genetic solutions of base population parents. The accuracy, stability, dispersion, and bias of model solutions were compared using the linear regression method. Addition of genomic information increased accuracy from 0.41 to 0.47 and stability from 0.68 to 0.71, while increasing bias slightly. Dispersion was within 0.10 of the ideal value (1.0) for all models. Although inclusion of metafounders did not strongly affect accuracy or stability and had mixed effects on bias, we nevertheless recommend the incorporation of metafounders in prediction models to represent the hierarchical genetic population structure of recently domesticated populations.

Population structure and genetic diversity of U.S. wheat varieties.

The United States is a major wheat producer with more than a century of wheat (Triticum aestivum L.) research and breeding. Using a panel of 753 historical and modern wheat varieties grown in the United States from the early 1800s to present day, we examined population structure and changes in genetic diversity. We used previously mapped high-quality single-nucleotide polymorphism (SNP) markers from the wheat 90K SNP array for genotyping. The wheat varieties had a slight hierarchical population structure based on growth habit and then by kernel color within spring varieties and by kernel hardness within winter varieties, which corresponds with geographical distribution of the varieties. Classifying varieties by market class, which is a combination of habit, hardness, and color, accounted for the greatest amount of variation (13.3%). We did not find evidence of decreased genetic diversity of either spring or winter varieties after the release of the first semidwarf wheat variety in 1961. On the contrary, northern and Pacific spring varieties, hard red spring (HRS), hard white spring (HWS), and soft white winter (SWW) had increases in both SNP and haplotype genetic diversity after 1961. The soft white spring (SWS) and soft red winter (SRW) market classes already had high genetic diversity in varieties before 1961 and showed some evidence of decreased diversity after 1961. Examination of temporal trends in genetic diversity also did not indicate long-term decline in diversity despite occasional fluctuations.

Population connectivity across a highly fragmented distribution: Phylogeography of the Chalcophaps doves

Chalcophaps is a morphologically conserved genus of ground-walking doves distributed from India to mainland China, south to Australia, and across the western Pacific to Vanuatu. Here, we reconstruct the evolutionary history of this genus using DNA sequence data from two nuclear genes and one mitochondrial gene, sampled from throughout the geographic range of Chalcophaps. We find support for three major evolutionary lineages in our phylogenetic reconstruction, each corresponding to the three currently recognized Chalcophaps species. Despite this general concordance, we identify discordant mitochondrial and nuclear ancestries in the subspecies C. longirostris timorensis, raising further questions about the evolutionary history of this Timor endemic population. Within each of the three species, we find evidence for isolation by distance or hierarchical population structure, indicating an important role for geography in the diversification of this genus. Despite being distributed broadly across a highly fragmented geographic region known as a hotspot for avian diversification, the Chalcophaps doves show modest levels of phenotypic and genetic diversity, a pattern potentially explained by strong population connectivity owing to high overwater dispersal capability.

The population genomic structure of green turtles (Chelonia mydas) suggests a warm-water corridor for tropical marine fauna between the Atlantic and Indian oceans during the last interglacial.

The occasional westward transport of warm water of the Agulhas Current, "Agulhas leakage", around southern Africa has been suggested to facilitate tropical marine connectivity between the Atlantic and Indian oceans, but the "Agulhas leakage" hypothesis does not explain the signatures of eastward gene flow observed in many tropical marine fauna. We investigated an alternative hypothesis: the establishment of a warm-water corridor during comparatively warm interglacial periods. The "warm-water corridor" hypothesis was investigated by studying the population genomic structure of Atlantic and Southwest Indian Ocean green turtles (N = 27) using 12,035 genome-wide single nucleotide polymorphisms (SNPs) obtained via ddRAD sequencing. Model-based and multivariate clustering suggested a hierarchical population structure with two main Atlantic and Southwest Indian Ocean clusters, and a Caribbean and East Atlantic sub-cluster nested within the Atlantic cluster. Coalescent-based model selection supported a model where Southwest Indian Ocean and Caribbean populations diverged from the East Atlantic population during the transition from the last interglacial period (130-115 thousand years ago; kya) to the last glacial period (115-90 kya). The onset of the last glaciation appeared to isolate Atlantic and Southwest Indian Ocean green turtles into three refugia, which subsequently came into secondary contact in the Caribbean and Southwest Indian Ocean when global temperatures increased after the Last Glacial Maximum. Our findings support the establishment of a warm-water corridor facilitating tropical marine connectivity between the Atlantic and Southwest Indian Ocean during warm interglacials.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.

HierCC: a multi-level clustering scheme for population assignments based on core genome MLST.

SummaryMotivationRoutine infectious disease surveillance is increasingly based on large-scale whole-genome sequencing databases. Real-time surveillance would benefit from immediate assignments of each genome assembly to hierarchical population structures. Here we present pHierCC, a pipeline that defines a scalable clustering scheme, HierCC, based on core genome multi-locus typing that allows incremental, static, multi-level cluster assignments of genomes. We also present HCCeval, which identifies optimal thresholds for assigning genomes to cohesive HierCC clusters. HierCC was implemented in EnteroBase in 2018 and has since genotyped >530 000 genomes from Salmonella, Escherichia/Shigella, Streptococcus, Clostridioides, Vibrio and Yersinia. Availability and implementation https://enterobase.warwick.ac.uk/ and Source code and instructions: https://github.com/zheminzhou/pHierCCSupplementary information Supplementary data are available at Bioinformatics online.

An elite-guided hierarchical differential evolution algorithm

Population structure has an impact on the performance of metaheuristic algorithms. To better improve the performance of differential evolution (DE), an elite-guided hierarchical differential evolution algorithm (EHDE) is proposed. First, an elite-guided hierarchical mutation mechanism is presented, which integrates elite elements into the hierarchical population structure. During each generation, the population is divided into three groups according to fitness values, each group playing a unique role in its hierarchy. The best individual on the top layer is used to avoid the local optimal by random reinitialization or Levy flight. The (k − 1) elite individuals on the middle layer focus on the local search around the best individual. The remaining non-elite individuals on the bottom layer pay more attention to a more considerable range search by the guidance of the k elite individuals. Second, to accommodate diverse optimization problems and seek the balance between exploration and exploitation, the adaptive strategy of EHDE control parameters has added the random component and the time-varying component. Finally, for the sake of evaluating the performance of EHDE, sensitivity analysis to the size of elite individuals, efficiency analysis of the control parameters adaptive strategy, and comparisons with nine advanced DE variants and three non-DE algorithms on 29 universal benchmark function in terms of convergence accuracy and convergence speed have been taken out. All the obtained results show that the proposed EHDE has excellent optimization performance.

Complex patterns of genetic and morphological differentiation in the Smallmouth Bass subspecies (Micropterus dolomieu dolomieu and M. d. velox) of the Central Interior Highlands

Due to geologic processes and recent anthropogenic introductions, patterns of genetic and morphological diversity within the Smallmouth Bass (Micropterus dolomieu), which are endemic to the central and eastern United States (USA), are poorly understood. We assessed genetic and morphological differentiation between the widespread Northern Smallmouth Bass (M. d. dolomieu) and the more restricted Neosho Smallmouth Bass (M. d. velox) where their ranges meet in the Central Interior Highlands ecoregion (CIH). Data from 14 microsatellite loci were used to conduct Structure and principal components analyses to evaluate diversity across populations and screen for hybridization with sympatric Spotted Bass (M. punctulatus). We also tested for morphological differences using five morphometric traits and one meristic trait. We found support for three genetic clusters corresponding to previously described taxonomic variation; five clusters largely corresponding to river systems; and nine clusters representing hierarchical population structure within both ranges. We found evidence of a unique genetic cluster in tributaries of the White River within the Northern Smallmouth Bass range and admixture between the subspecies throughout the Neosho range. We also found evidence of morphological differentiation between subspecies; Neosho Smallmouth Bass exhibited larger head length than Northern Smallmouth Bass relative to total length, and there was a significant interaction of subspecies and orbital length, possibly indicating differential growth patterns between subspecies. Our results reveal multiple levels of divergence, suggesting the CIH harbors greater and more complex Smallmouth Bass diversity than previously thought.

UNITED IMMUNOLOGICAL FORUM: CURRENT TRENDS IN THE DEVELOPMENT OF FUNDAMENTAL AND APPLIED ONCOIMMUNOLOGY (NOVOSIBIRSK, 2019)

The work provides information on the results of the Joint Immunological Forum, which was held from June 24 to 29, 2019 in Novosibirsk. The modern directions of the development of fundamental and applied immunology are analyzed. Particular attention is paid to the discussion of the most issues identified in the section Immunopathogenetic bases of tumor growth, which also presented the results of studies conducted at the National Medical Research Centre for Oncology of the Ministry of Health of Russia for the study of isolation and study of the biological properties of tumor stem cells. Noteworthy are the new advances in modern immunology, which clarify the hierarchical structure of lymphocyte populations, with the separation of various minor subpopulations based on the phenotypic, molecular genetic and functional properties of cells, whose role in ensuring the integrity of the body has not been fully studied. In addition to theoretical reports, during this Forum the results of using new methodological approaches to study the structural and functional organization of individual links of innate and adaptive immunity both under model conditions and during the development of various human diseases were presented, the most promising ways to improve analytical and technological platforms were identified. In the crayfish of the Forum, several advanced training programs were implemented for employees of various levels of practical health care and fundamental science.

First steps towards the identification of evolutionarily significant units in Mexican native trout: An assessment of microsatellite variation

In the remote high-elevation streams and rivers in the northwest of Mexico, there is an endemic and highly diverse group of native trout related to the rainbow trout lineage (Oncorhynchus mykiss). Currently, only two species have been formally described from this trout complex, the Mexican golden trout, O. chrysogaster, and the San Pedro Martir rainbow trout, O. m. nelsoni, but only the former occurs in the Sierra Madre Occidental (SMO). South and north of the O. chrysogaster range, there are other undescribed trout that have been recently discovered. Given the various human impacts on aquatic systems occupied by these Mexican native trout, the lack of taxonomic descriptions and the imminent risk of losing populations without their even being documented, we sought to infer a hierarchical population structure of the Mexican native trout, determine the geographic patterns of genetic diversity and gene flow among them, and identify evolutionarily significant units for a prompt conservation strategy. The allelic diversity across eleven microsatellite loci was used to assess genetic divergences of wild trout collected throughout its range in Mexico. Our results reject the old hypothesis that many of the trout in Mexico are introduced or translocated hatchery-reared rainbow trout and strongly support the new hypothesis of distinct genetic populations and lineages in the northwest of Mexico. According to the different statistical analysis of the genetic data, we distinguish between species and subspecies (O. chrysogaster, O. mykiss, and O. m. nelsoni) of the other regional lineages of undescribed trout (northern and southern) inhabiting the SMO. We also found a strong genetic substructure within these large clusters. We propose the presence of seven Evolutionarily Significant Units in the drainage basins of the SMO (two occur in the northern region, four in the southern, and one for the area of distribution of O. chrysogaster) and one in north-western Baja California. Based in microsatellite data a minimum gene flow among populations is evidenced, resulting the presence the highly divergent populations.

Confidently identifying the correct K value using the ΔK method: When does K=2?

Populations delineated based on genetic data are commonly used for wildlife conservation and management. Many studies use the program structure combined with the ΔK method to identify the most probable number of populations (K). We recently found K=2 was identified more often when studies used ΔK compared to studies that did not. We suggested two reasons for this: hierarchical population structure leads to underestimation, or the ΔK method does not evaluate K=1 causing an overestimation. The present contribution aims to develop a better understanding of the limits of the method using one, two and three population simulations across migration scenarios. From these simulations we identified the "best K" using model likelihood and ΔK. Our findings show that mean probability plots and ΔK are unable to resolve the correct number of populations once migration rate exceeds 0.005. We also found a strong bias towards selecting K=2 using the ΔK method. We used these data to identify the range of values where the ΔK statistic identifies a value of K that is not well supported. Finally, using the simulations and a review of empirical data, we found that the magnitude of ΔK corresponds to the level of divergence between populations. Based on our findings, we suggest researchers should use the ΔK method cautiously; they need to report all relevant data, including the magnitude of ΔK, and an estimate of connectivity for the research community to assess whether meaningful genetic structure exists within the context of management and conservation.