Combining markers into haplotypes can improve population structure inference

Are rabid raccoons (Procyon lotor) ready for the rapture? Determining the geographic origin of rabies virus-infected raccoons using RADcapture and microhaplotypes

North America is recognized for the exceptional richness of rabies virus (RV) wildlife reservoir species. Management of RV is accomplished through vaccination targeting mesocarnivore reservoir populations, such as the raccoon (Procyon lotor) in Eastern North America. Raccoons are a common generalist species, and populations may reach high densities in developed areas, which can result in contact with humans and pets with potential exposures to the raccoon variant of RV throughout the eastern United States. Understanding the spatial movement of RV by raccoon populations is important for monitoring and refining strategies supporting the landscape-level control and local elimination of this lethal zoonosis. We developed a high-throughput genotyping panel for raccoons based on hundreds of microhaplotypes to identify population structure and genetic diversity relevant to rabies management programs. Throughout the eastern United States, we identified hierarchical population genetic structure with clusters that were connected through isolation-by-distance. We also illustrate that this genotyping approach can be used to support real-time management priorities by identifying the geographic origin of a rabid raccoon that was collected in an area of the United States that had been raccoon RV-free for 8 years. The results from this study and the utility of the microhaplotype panel and genotyping method will provide managers with information on raccoon ecology that can be incorporated into future management decisions.

Non-random mating within an Island rookery of Hawaiian hawksbill turtles: demographic discontinuity at a small coastline scale

Hawksbill sea turtles (Eretmochelys imbricata) from the Hawaiian archipelago form a small and genetically isolated population, consisting of only a few tens of individuals breeding annually. Most females nest on the island of Hawai'i, but little is known about the demographics of this rookery. This study used genetic relatedness, inferred from 135 microhaplotype markers, to determine breeding sex-ratios, estimate female nesting frequency and assess relationships between individuals nesting on different beaches. Samples were collected during the 2017 nesting season and final data included 13 nesting females and 1002 unhatched embryos, salvaged from 41 nests, of which 13 had no observed mother. Results show that most females used a single nesting beach laying 1-5 nests each. From female and offspring alleles, the paternal genotypes of 12 breeding males were reconstructed and many showed high relatedness to their mates. Pairwise relatedness of offspring revealed one instance of polygyny but otherwise suggested a 1 : 1 breeding-sex ratio. Relatedness analysis and spatial-autocorrelation of genotypes indicate that turtles from different nesting areas do not regularly interbreed, suggesting that strong natal homing tendencies in both sexes result in non-random mating across the study area. Complexes of nearby nesting beaches also showed unique patterns of inbreeding across loci, further indicating that Hawaiian hawksbill turtles have demographically discontinuous nesting populations separated by only tens of km.

Genetic evidence of transoceanic migration of the small brown planthopper between China and Japan

BACKGROUND

The small brown planthopper, Laodelphax striatellus (Fallén), is an important pest of rice. It is suspected of migrating over the sea from China to Japan. However, where in China it comes from and how it affects Japanese populations remain unclear.

RESULTS

Here, we studied the genetic structure of 15 L. striatellus populations sampled from Japan and China using single nucleotide polymorphisms generated by the double digest restriction site-associated DNA sequencing technique. We found weak genetic differentiation between the Chinese and Japanese populations. Our data revealed migration signals of L. striatellus from China to southern and northern Japan. However, the source regions of the immigrants remain unclear due to the low genetic differentiation between populations. Our results also pointed to the possibility of backward gene flow from Japanese to Chinese populations. We suspect that the south-eastern wind associated with the East Asian summer monsoon may facilitate the reverse migration of L. striatellus from Japan to China. Interestingly, we found that the X chromosome displayed relatively higher genetic differentiation among populations and suffered more intensive selection pressure than autosomes.

CONCLUSION

We provide genetic evidence of transoceanic migration of L. striatellus from China to Japan and found that the X chromosome can aid the deciphering of the migration trajectories of species with low genetic differentiation. These findings have implications for forecasting the outbreak of this pest and also provide insights into how to improve the tracking of the migration routes of small insects via population genomics. © 2022 Society of Chemical Industry.

Variants in linkage status at D5S818 detected by multiple STR kits comparison and Sanger sequencing

Background

D5S818 discrepancies have been reported in forensic parental testing due to null alleles. However, more cases may be ignored since proportional null alleles were missed without detection of heredity discrepancy between parents and offspring.

Results

In this study, null allele 12 at D5S818 was detected by the PowerPlex^® 21 System with a higher occurrence rate on the basis of review on 2824 samples from the 1282 routine cases in Chinese Han population. Sequencing results revealed novel variant of guanine (G) into adenine (A) in the 7th [AGAT] repeats in the core repeat region accompanied by rs1187948322 in the samples with null allele 12.

Conclusions

Forensic STR typing may benefit from this discovery: (1) primer design of CE profiling system could be improved for sensitive population and (2) polymorphic information could be enriched for the accuracy and precision of NGS genotyping system. Peak area of D5S818 was also analyzed through different commercial STR kits. It is suggested that more attention should be paid on observed homozygosity with reduced peak area, especially for the samples from Chinese Han population.

Haplotype structure in commercial maize breeding programs in relation to key founder lines

High-density haplotype analysis revealed significant haplotype sharing between ex-PVPs registered from 1976 to 1992 and key maize founders, and uncovered similarities and differences in haplotype sharing patterns by company and heterotic group. Proprietary inbreds developed by the private seed industry have been the major source for driving genetic gain in successful North American maize hybrids for decades. Much of the history of industry germplasm can be traced back to key founder lines, some of which were pivotal in the development of prominent heterotic groups. Previous studies have summarized pedigree-based relationships, genetic diversity and population structure among commercial inbreds with expired Plant Variety Protection (ex-PVP). However, less is known about the extent of haplotype sharing between historical founders and ex-PVPs. A better understanding of the relationships between founders and ex-PVPs provides insight into the haplotype and heterotic group structure among industry germplasm. We performed high-density haplotype analysis with 11.3 million SNPs on 212 maize inbreds, which included 157 ex-PVPs registered 1976-1992 and 55 public lines relevant to PVPs. Among these lines were 12 key founders identified in literature review: 207, A632, B14, B37, B73, LH123HT, LH82, Mo17, Oh43, OH7, PHG39 and Wf9. Our results revealed that, on average, 81.6% of an ex-PVP's genome is shared with at least 1 of these 12 founder lines and more than half when limited to B73, Mo17 and 207. Quantifiable similarities and contrasts among heterotic groups and major US seed industry companies were also observed. The results from this study provide high-resolution haplotype data on ex-PVP germplasm, confirm founder relationship trends observed in previous studies, uncover region-specific haplotype structure differences and demonstrate how haplotype sharing analysis can be used as a tool to explore germplasm diversity.

Resolving the cause of recurrent Plasmodium vivax malaria probabilistically

Relapses arising from dormant liver-stage Plasmodium vivax parasites (hypnozoites) are a major cause of vivax malaria. However, in endemic areas, a recurrent blood-stage infection following treatment can be hypnozoite-derived (relapse), a blood-stage treatment failure (recrudescence), or a newly acquired infection (reinfection). Each of these requires a different prevention strategy, but it was not previously possible to distinguish between them reliably. We show that individual vivax malaria recurrences can be characterised probabilistically by combined modelling of time-to-event and genetic data within a framework incorporating identity-by-descent. Analysis of pooled patient data on 1441 recurrent P. vivax infections in 1299 patients on the Thailand-Myanmar border observed over 1000 patient follow-up years shows that, without primaquine radical curative treatment, 3 in 4 patients relapse. In contrast, after supervised high-dose primaquine only 1 in 40 relapse. In this region of frequent relapsing P. vivax, failure rates after supervised high-dose primaquine are significantly lower (∼3%) than estimated previously.

The Relationship Between Haplotype-Based FST and Haplotype Length

The population-genetic statistic [Formula: see text] is used widely to describe allele frequency distributions in subdivided populations. The increasing availability of DNA sequence data has recently enabled computations of [Formula: see text] from sequence-based "haplotype loci." At the same time, theoretical work has revealed that [Formula: see text] has a strong dependence on the underlying genetic diversity of a locus from which it is computed, with high diversity constraining values of [Formula: see text] to be low. In the case of haplotype loci, for which two haplotypes that are distinct over a specified length along a chromosome are treated as distinct alleles, genetic diversity is influenced by haplotype length: longer haplotype loci have the potential for greater genetic diversity. Here, we study the dependence of [Formula: see text] on haplotype length. Using a model in which a haplotype locus is sequentially incremented by one biallelic locus at a time, we show that increasing the length of the haplotype locus can either increase or decrease the value of [Formula: see text], and usually decreases it. We compute [Formula: see text] on haplotype loci in human populations, finding a close correspondence between the observed values and our theoretical predictions. We conclude that effects of haplotype length are valuable to consider when interpreting [Formula: see text] calculated on haplotypic data.

Estimating Relatedness Between Malaria Parasites

Understanding the relatedness of individuals within or between populations is a common goal in biology. Increasingly, relatedness features in genetic epidemiology studies of pathogens. These studies are relatively new compared to those in humans and other organisms, but are important for designing interventions and understanding pathogen transmission. Only recently have researchers begun to routinely apply relatedness to apicomplexan eukaryotic malaria parasites, and to date have used a range of different approaches on an ad hoc basis. Therefore, it remains unclear how to compare different studies and which measures to use. Here, we systematically compare measures based on identity-by-state (IBS) and identity-by-descent (IBD) using a globally diverse data set of malaria parasites, Plasmodium falciparum and P. vivax, and provide marker requirements for estimates based on IBD. We formally show that the informativeness of polyallelic markers for relatedness inference is maximized when alleles are equifrequent. Estimates based on IBS are sensitive to allele frequencies, which vary across populations and by experimental design. For portability across studies, we thus recommend estimates based on IBD. To generate estimates with errors below an arbitrary threshold of 0.1, we recommend ∼100 polyallelic or 200 biallelic markers. Marker requirements are immediately applicable to haploid malaria parasites and other haploid eukaryotes. C.I.s facilitate comparison when different marker sets are used. This is the first attempt to provide rigorous analysis of the reliability of, and requirements for, relatedness inference in malaria genetic epidemiology. We hope it will provide a basis for statistically informed prospective study design and surveillance strategies.

Investigation of RFLP Haplotypes β-Globin Gene Cluster in Beta-Thalassemia Patients in Central Iran

Introduction: Beta-thalassemia is one of the most prevalent inherited blood diseases among Iranians. The aim of this study was to elucidate the chromosomal background of beta-thalassemia mutations in Esfahan province, Iran. Materials and Methods: In this study, we investigated three frequent mutations (c.315+1G>A, c.93-21G>A and c.92+5G>C in β-globin gene, the frequency of RFLP haplotypes, and LD between markers at β-globin gene cluster) in 150 beta-thalassemia patients and 50 healthy individuals. The molecular and population genetic investigations were performed on RFLP markers HindIII in the c.315+1G>A of Gγ (HindIIIG) and Aγ (HindIIIA) genes, AvaII in the c.315+1G>A of β-globin gene and BamHI 3' to the β-globin gene. All statistical analyses were performed using Power Marker software and SISA server. Results: Fifty percent of beta-thalasemia patients were associated with these mutations. Haplotype I was the most prevalent haplotype among beta-thalassemia patients (39.33%) and normal individuals (46%). The commonest c.315+1G>A mutation in our population was tightly linked with haplotype III (43.75%) and haplotype I (31.25%). The second prevalent mutation, c.92+5G>C, was 90%, 6.66%, and 3.33% in linkage disequilibrium with haplotypes I, VII, and III, respectively. The c.93-21G>A mutation indicated a strong association with haplotype I (80%). Conclusion: Our study participants like beta-thalassemia patients from Kermanshah province was found to possess a similar haplotype background for common mutations. The emergence of most prevalent mutations on chromosomes with different haplotypes can be explained by gene conversion and recombination. High linkage of a mutation with specific haplotype is consistent with the hypothesis that chromosomes carrying beta-thalassemia mutations experienced positive selection pressure, probably because of the protection against malaria experienced by beta-thalassemia carriers.

Anadromy Redux? Genetic Analysis to Inform Development of an Indigenous American River Steelhead Broodstock

The construction of dams and water diversions has severely limited access to spawning habitat for anadromous fishes. To mitigate for these impacts, hatchery programs rear and release millions of juvenile salmonids, including steelhead, the anadromous ecotype of the species Oncorhynchus mykiss. These programs sometimes use nonindigenous broodstock sources that may have negative effects on wild populations. In California, however, only one anadromous fish hatchery program currently uses nonnative broodstock: the steelhead program at Nimbus Fish Hatchery on the American River, a tributary of the Sacramento River in the California Central Valley. The goal of this study was to determine if potentially appropriate sources to replace the broodstock for the Nimbus Hatchery steelhead program exist in the Upper American River, above Nimbus and Folsom dams. We show that all Upper American River O. mykiss sampled share ancestry with other populations in the Central Valley steelhead distinct population segment, with limited introgression from out-of-basin sources in some areas. Furthermore, some Upper American River populations retain adaptive genomic variation associated with a migratory life history, supporting the hypothesis that these populations display adfluvial migratory behavior. Together, these results provide insights into the evolution of trout populations above barrier dams. We conclude that some Upper American River O. mykiss populations represent genetically appropriate sources from which fisheries managers could potentially develop a new broodstock for the Nimbus Hatchery steelhead program to reestablish a native anadromous population in the Lower American River and contribute to recovery of the threatened Central Valley steelhead distinct population segment.

New perspectives on multilocus ancestry informativeness

We present an axiomatic approach for multilocus informativeness measures for determining the amount of information that a set of polymorphic genetic markers provides about individual ancestry. We then reveal several surprising properties of a decision-theoretic based measure that is consistent with the set of proposed criteria for multilocus informativeness. In particular, these properties highlight the interplay between information originating from population priors and the information extractable from the population genetic variants. This analysis then reveals a certain deficiency of mutual information based multilocus informativeness measures when such population priors are incorporated. Finally, we analyse and quantify the inevitable inherent decrease in informativeness due to learning from finite population samples.

Inferring biogeographic ancestry with compound markers of slow and fast evolving polymorphisms

Bio-geographic ancestry is an area of considerable interest in the medical genetics, anthropology and forensics. Although genome-wide panels are ideal as they provide dense genotyping data, small sets of ancestry informative marker provide a cost-effective way to investigate genetic ancestry and population structure. Here, we investigate the performance of a reduced marker set that combine different types of autosomal markers through haplotype analysis. In particular, recently described DIP-STR markers should offer the advantage of comprising both, low mutation rate Indels (DIPs), to study human history over longer time scale; and high mutation rate STRs, to trace relatively recent demographic events. In this study, we assessed the ability of an initial set of 23 DIP-STRs to distinguish major population groups using the HGDP-CEPH reference samples. The results obtained applying the STRUCTURE algorithm show that the discrimination capacity of the DIP-STRs is comparable to currently used small-scale ancestry informative markers by approaching seven major demographic groups. Yet, the DIP-STRs show an improved success rate in assigning individuals to populations of Europe and Middle East. These data show a remarkable ability of a preliminary set of 23 DIP-STR markers to infer major biogeographic origins. A novel set of DIP-STRs preselected to contain ancestry information should lead to further improvements.

Human cytochrome P450 2B6 genetic variability in Botswana: a case of haplotype diversity and convergent phenotypes

Identification of inter-individual variability for drug metabolism through cytochrome P450 2B6 (CYP2B6) enzyme is important for understanding the differences in clinical responses to malaria and HIV. This study evaluates the distribution of CYP2B6 alleles, haplotypes and inferred metabolic phenotypes among subjects with different ethnicity in Botswana. A total of 570 subjects were analyzed for CYP2B6 polymorphisms at position 516 G > T (rs3745274), 785 A > G (rs2279343) and 983 T > C (rs28399499). Samples were collected in three districts of Botswana where the population belongs to Bantu (Serowe/Palapye and Chobe) and San-related (Ghanzi) ethnicity. The three districts showed different haplotype composition according to the ethnic background but similar metabolic inferred phenotypes, with 59.12%, 34.56%, 2.10% and 4.21% of the subjects having, respectively, an extensive, intermediate, slow and rapid metabolic profile. The results hint at the possibility of a convergent adaptation of detoxifying metabolic phenotypes despite a different haplotype structure due to the different genetic background. The main implication is that, while there is substantial homogeneity of metabolic inferred phenotypes among the country, the response to drugs metabolized via CYP2B6 could be individually associated to an increased risk of treatment failure and toxicity. These are important facts since Botswana is facing malaria elimination and a very high HIV prevalence.

The fine-scale genetic structure and evolution of the Japanese population

The contemporary Japanese populations largely consist of three genetically distinct groups—Hondo, Ryukyu and Ainu. By principal-component analysis, while the three groups can be clearly separated, the Hondo people, comprising 99% of the Japanese, form one almost indistinguishable cluster. To understand fine-scale genetic structure, we applied powerful haplotype-based statistical methods to genome-wide single nucleotide polymorphism data from 1600 Japanese individuals, sampled from eight distinct regions in Japan. We then combined the Japanese data with 26 other Asian populations data to analyze the shared ancestry and genetic differentiation. We found that the Japanese could be separated into nine genetic clusters in our dataset, showing a marked concordance with geography; and that major components of ancestry profile of Japanese were from the Korean and Han Chinese clusters. We also detected and dated admixture in the Japanese. While genetic differentiation between Ryukyu and Hondo was suggested to be caused in part by positive selection, genetic differentiation among the Hondo clusters appeared to result principally from genetic drift. Notably, in Asians, we found the possibility that positive selection accentuated genetic differentiation among distant populations but attenuated genetic differentiation among close populations. These findings are significant for studies of human evolution and medical genetics.

Fine-Scale Genetic Structure in Finland

Coupling dense genotype data with new computational methods offers unprecedented opportunities for individual-level ancestry estimation once geographically precisely defined reference data sets become available. We study such a reference data set for Finland containing 2376 such individuals from the FINRISK Study survey of 1997 both of whose parents were born close to each other. This sampling strategy focuses on the population structure present in Finland before the 1950s. By using the recent haplotype-based methods ChromoPainter (CP) and FineSTRUCTURE (FS) we reveal a highly geographically clustered genetic structure in Finland and report its connections to the settlement history as well as to the current dialectal regions of the Finnish language. The main genetic division within Finland shows striking concordance with the 1323 borderline of the treaty of Nöteborg. In general, we detect genetic substructure throughout the country, which reflects stronger regional genetic differences in Finland compared to, for example, the UK, which in a similar analysis was dominated by a single unstructured population. We expect that similar population genetic reference data sets will become available for many more populations in the near future with important applications, for example, in forensic genetics and in genetic association studies. With this in mind, we report those extensions of the CP + FS approach that we found most useful in our analyses of the Finnish data.

Recent advances in the study of fine-scale population structure in humans

Empowered by modern genotyping and large samples, population structure can be accurately described and quantified even when it only explains a fraction of a percent of total genetic variance. This is especially relevant and interesting for humans, where fine-scale population structure can both confound disease-mapping studies and reveal the history of migration and divergence that shaped our species' diversity. Here we review notable recent advances in the detection, use, and understanding of population structure. Our work addresses multiple areas where substantial progress is being made: improved statistics and models for better capturing differentiation, admixture, and the spatial distribution of variation; computational speed-ups that allow methods to scale to modern data; and advances in haplotypic modeling that have wide ranging consequences for the analysis of population structure. We conclude by outlining four important open challenges: the limitations of discrete population models, uncertainty in individual origins, the incorporation of both fine-scale structure and ancient DNA in parametric models, and the development of efficient computational tools, particularly for haplotype-based methods.

HaploPOP: a software that improves population assignment by combining markers into haplotypes

Background

In ecology and forensics, some population assignment techniques use molecular markers to assign individuals to known groups. However, assigning individuals to known populations can be difficult if the level of genetic differentiation among populations is small. Most assignment studies handle independent markers, often by pruning markers in Linkage Disequilibrium (LD), ignoring the information contained in the correlation among markers due to LD.

Results

To improve the accuracy of population assignment, we present an algorithm, implemented in the HaploPOP software, that combines markers into haplotypes, without requiring independence. The algorithm is based on the Gain of Informativeness for Assignment that provides a measure to decide if a pair of markers should be combined into haplotypes, or not, in order to improve assignment. Because complete exploration of all possible solutions for constructing haplotypes is computationally prohibitive, our approach uses a greedy algorithm based on windows of fixed sizes. We evaluate the performance of HaploPOP to assign individuals to populations using a split-validation approach. We investigate both simulated SNPs data and dense genotype data from individuals from Spain and Portugal.

Conclusions

Our results show that constructing haplotypes with HaploPOP can substantially reduce assignment error. The HaploPOP software is freely available as a command-line software at www.ieg.uu.se/Jakobsson/software/HaploPOP/.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

Nonstationary patterns of isolation-by-distance: inferring measures of local genetic differentiation with Bayesian kriging

Patterns of isolation-by-distance (IBD) arise when population differentiation increases with increasing geographic distances. Patterns of IBD are usually caused by local spatial dispersal, which explains why differences of allele frequencies between populations accumulate with distance. However, spatial variations of demographic parameters such as migration rate or population density can generate nonstationary patterns of IBD where the rate at which genetic differentiation accumulates varies across space. To characterize nonstationary patterns of IBD, we infer local genetic differentiation based on Bayesian kriging. Local genetic differentiation for a sampled population is defined as the average genetic differentiation between the sampled population and fictive neighboring populations. To avoid defining populations in advance, the method can also be applied at the scale of individuals making it relevant for landscape genetics. Inference of local genetic differentiation relies on a matrix of pairwise similarity or dissimilarity between populations or individuals such as matrices of F_ST between pairs of populations. Simulation studies show that maps of local genetic differentiation can reveal barriers to gene flow but also other patterns such as continuous variations of gene flow across habitat. The potential of the method is illustrated with two datasets: single nucleotide polymorphisms from human Swedish populations and dominant markers for alpine plant species.

Resequencing data indicate a modest effect of domestication on diversity in barley: a cultigen with multiple origins

The levels of diversity and extent of linkage disequilibrium in cultivated species are largely determined by diversity in their wild progenitors. We report a comparison of nucleotide sequence diversity in wild and cultivated barley (Hordeum vulgare ssp. spontaneum and ssp. vulgare) at 7 nuclear loci totaling 9296bp, using sequence from Hordeum bulbosum to infer the ancestral state of mutations. The sample includes 36 accessions of cultivated barley, including 23 landraces (cultivated forms not subject to modern breeding) and 13 cultivated lines and genetic stocks compared to either 25 or 45 accessions of wild barley for the same loci. Estimates of nucleotide sequence diversity indicate that landraces retain >80% of the diversity in wild barley. The primary population structure in wild barley, which divides the species into eastern and western populations, is reflected in significant differentiation at all loci in wild accessions and at 3 of 7 loci in landraces. "Oriental" landraces have slightly higher diversity than "Occidental" landraces. Genetic assignment suggests more admixture from Occidental landraces into Oriental landraces than the converse, which may explain this difference. Based on θπ for silent sites, modern western cultivars have ~73% of the diversity found in landraces and ~71% of the diversity in wild barley.

Tracing the geographic origins of weedy Ipomoea purpurea in the southeastern United States

Ipomoea purpurea (common morning glory) is an annual vine native to Mexico that is well known for its large, showy flowers. Humans have spread morning glories worldwide, owing to the horticultural appeal of morning glory flowers. Ipomoea purpurea is an opportunistic colonizer of disturbed habitats including roadside and agricultural settings, and it is now regarded as a noxious weed in the Southeastern US. Naturalized populations in the Southeastern United States are highly polymorphic for a number of flower color morphs, unlike native Mexican populations that are typically monomorphic for the purple color morph. Although I. purpurea was introduced into the United States from Mexico, little is known about the specific geographic origins of US populations relative to the Mexican source. We use resequencing data from 11 loci and 30 I. purpurea accessions collected from the native range of the species in Central and Southern Mexico and 8 accessions from the Southeastern United States to infer likely geographic origins in Mexico. Based on genetic assignment analysis, haplotype composition, and the degree of shared polymorphism, I. purpurea samples from the Southeastern United States are genetically most similar to samples from the Valley of Mexico and Veracruz State. This supports earlier speculation that I. purpurea in the Southeastern United States was likely to have been introduced by European colonists from sources in Central Mexico.

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over-exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded FST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between-ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17-SNP panel was developed in Atlantic BFT by cross-species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (FST  = 0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.

Genome-wide diversity in the levant reveals recent structuring by culture

The Levant is a region in the Near East with an impressive record of continuous human existence and major cultural developments since the Paleolithic period. Genetic and archeological studies present solid evidence placing the Middle East and the Arabian Peninsula as the first stepping-stone outside Africa. There is, however, little understanding of demographic changes in the Middle East, particularly the Levant, after the first Out-of-Africa expansion and how the Levantine peoples relate genetically to each other and to their neighbors. In this study we analyze more than 500,000 genome-wide SNPs in 1,341 new samples from the Levant and compare them to samples from 48 populations worldwide. Our results show recent genetic stratifications in the Levant are driven by the religious affiliations of the populations within the region. Cultural changes within the last two millennia appear to have facilitated/maintained admixture between culturally similar populations from the Levant, Arabian Peninsula, and Africa. The same cultural changes seem to have resulted in genetic isolation of other groups by limiting admixture with culturally different neighboring populations. Consequently, Levant populations today fall into two main groups: one sharing more genetic characteristics with modern-day Europeans and Central Asians, and the other with closer genetic affinities to other Middle Easterners and Africans. Finally, we identify a putative Levantine ancestral component that diverged from other Middle Easterners ∼23,700–15,500 years ago during the last glacial period, and diverged from Europeans ∼15,900–9,100 years ago between the last glacial warming and the start of the Neolithic.

Population stratification caused by nonrandom mating between groups of the same species is often due to geographical distances leading to physical separation followed by genetic drift of allele frequencies in each group. In humans, population structures are also often driven by geographical barriers or distances; however, humans might also be structured by abstract factors such as culture, a consequence of their reasoning and self-awareness. Religion in particular, is one of the unusual conceptual factors that can drive human population structures. This study explores the Levant, a region flanked by the Middle East and Europe, where individual and population relationships are still strongly influenced by religion. We show that religious affiliation had a strong impact on the genomes of the Levantines. In particular, conversion of the region's populations to Islam appears to have introduced major rearrangements in populations' relations through admixture with culturally similar but geographically remote populations, leading to genetic similarities between remarkably distant populations like Jordanians, Moroccans, and Yemenis. Conversely, other populations, like Christians and Druze, became genetically isolated in the new cultural environment. We reconstructed the genetic structure of the Levantines and found that a pre-Islamic expansion Levant was more genetically similar to Europeans than to Middle Easterners.

Inference of population structure using dense haplotype data

The advent of genome-wide dense variation data provides an opportunity to investigate ancestry in unprecedented detail, but presents new statistical challenges. We propose a novel inference framework that aims to efficiently capture information on population structure provided by patterns of haplotype similarity. Each individual in a sample is considered in turn as a recipient, whose chromosomes are reconstructed using chunks of DNA donated by the other individuals. Results of this "chromosome painting" can be summarized as a "coancestry matrix," which directly reveals key information about ancestral relationships among individuals. If markers are viewed as independent, we show that this matrix almost completely captures the information used by both standard Principal Components Analysis (PCA) and model-based approaches such as STRUCTURE in a unified manner. Furthermore, when markers are in linkage disequilibrium, the matrix combines information across successive markers to increase the ability to discern fine-scale population structure using PCA. In parallel, we have developed an efficient model-based approach to identify discrete populations using this matrix, which offers advantages over PCA in terms of interpretability and over existing clustering algorithms in terms of speed, number of separable populations, and sensitivity to subtle population structure. We analyse Human Genome Diversity Panel data for 938 individuals and 641,000 markers, and we identify 226 populations reflecting differences on continental, regional, local, and family scales. We present multiple lines of evidence that, while many methods capture similar information among strongly differentiated groups, more subtle population structure in human populations is consistently present at a much finer level than currently available geographic labels and is only captured by the haplotype-based approach. The software used for this article, ChromoPainter and fineSTRUCTURE, is available from http://www.paintmychromosomes.com/.