Genetic Substructure Research Articles

Genomic insight into the population history of Central Han Chinese

Background In recent decades, considerable attention has been paid to exploring the population genetic characteristics of Han Chinese, mainly documenting a north-south genetic substructure. However, the central Han Chinese have been largely underrepresented in previous studies. Aim To infer a comprehensive understanding of the homogenisation process and population history of Han Chinese. Subjects and methods We collected samples from 122 Han Chinese from seven counties of Hubei province in central China and genotyped 534,000 genome-wide SNPs. We compared Hubei Han with both ancient and present-day Eurasian populations using Principal Component Analysis, ADMIXTURE, f statistics, qpWave and qpAdm. Results We observed Hubei Han Chinese are at a genetically intermediate position on the north-south Han Chinese cline. We have not detected any significant genetic substructure in the studied groups from seven different counties. Hubei Han show significant evidence of genetic admixture deriving about 63% of ancestry from Tai-Kadai or Austronesian-speaking southern indigenous groups and 37% from Tungusic or Mongolic related northern populations. Conclusions The formation of Han Chinese has involved extensive admixture with Tai-Kadai or Austronesian-speaking populations in the south and Tungusic or Mongolic speaking populations in the north. The convenient transportation and central location of Hubei make it the key region for the homogenisation of Han Chinese.

Genetic assessment reveals no population substructure and divergent regional and sex-specific histories in the Chachapoyas from northeast Peru.

Many native populations in South America have been severely impacted by two relatively recent historical events, the Inca and the Spanish conquest. However decisive these disruptive events may have been, the populations and their gene pools have been shaped markedly also by the history prior to the conquests. This study focuses mainly on the Chachapoya peoples that inhabit the montane forests on the eastern slopes of the northern Peruvian Andes, but also includes three distinct neighboring populations (the Jívaro, the Huancas and the Cajamarca). By assessing mitochondrial, Y-chromosomal and autosomal diversity in the region, we explore questions that have emerged from archaeological and historical studies of the regional culture (s). These studies have shown, among others, that Chachapoyas was a crossroads for Coast-Andes-Amazon interactions since very early times. In this study, we examine the following questions: 1) was there pre-Hispanic genetic population substructure in the Chachapoyas sample? 2) did the Spanish conquest cause a more severe population decline on Chachapoyan males than on females? 3) can we detect different patterns of European gene flow in the Chachapoyas region? and, 4) did the demographic history in the Chachapoyas resemble the one from the Andean area? Despite cultural differences within the Chachapoyas region as shown by archaeological and ethnohistorical research, genetic markers show no significant evidence for past or current population substructure, although an Amazonian gene flow dynamic in the northern part of this territory is suggested. The data also indicates a bottleneck c. 25 generations ago that was more severe among males than females, as well as divergent population histories for populations in the Andean and Amazonian regions. In line with previous studies, we observe high genetic diversity in the Chachapoyas, despite the documented dramatic population declines. The diverse topography and great biodiversity of the northeastern Peruvian montane forests are potential contributing agents in shaping and maintaining the high genetic diversity in the Chachapoyas region.

Genetic substructure and admixture of Mongolians and Kazakhs inferred from genome-wide array genotyping

Background Mongolian populations are widely distributed geographically, showing abundant ethnic diversity with geographic and tribal differences. Aim To infer the genetic substructure, admixture and ancient genetic sources of Mongolians together with Kazakhs Subjects and methods We genotyped more than 690,000 genome-wide SNPs from 33 Mongolian and Chinese Kazakh individuals and compared these with both ancient and present-day Eurasian populations using Principal Component Analysis (PCA), ADMIXTURE, Refine-IBD, f statistics, qpWave and qpAdm. Results We found genetic substructures within Mongolians corresponding to Ölöd, Chahar, and Inner Mongolian clusters, which was consistent with tribe classifications. Mongolian and Kazakh groups derived about 6–40% of West Eurasian related ancestry, most likely from Bronze Age Steppe populations. The East Asian related ancestry in Mongolian and Kazakh groups was well represented by the Neolithic DevilsCave related nomadic lineage, comprising 42–64% of studied groups. We also detected 10–51% of Han Chinese related ancestry in Mongolian and Kazakh groups, especially in Inner Mongolians. The average admixture times for Inner Mongolian, Mongolian_Chahar, Mongolian_Ölöd and Chinese Kazakh were about 1381, 626, 635 and 632 years ago, respectively, with Han and French as the sources. Conclusion The DevilsCave related ancestry was once widespread westwards covering a wide geographical range from Far East Russia to the Mongolia Plateau. The formation of present-day Mongolic and Turkic-speaking populations has also received genetic influence from agricultural expansion.

A Coordinated Suite of Wild-Introgression Lines in Indica and Japonica Elite Backgrounds.

Rice, Oryza sativa L., is a cultivated, inbreeding species that serves as the staple food for the largest number of people on earth. It has two strongly diverged varietal groups, Indica and Japonica, which result from a combination of natural and human selection. The genetic divergence of these groups reflects the underlying population structure of their wild ancestors, and suggests that a pre-breeding strategy designed to take advantage of existing genetic, geographic and ecological substructure may provide a rational approach to the utilization of crop wild ancestors in plant improvement. Here we describe the coordinated development of six introgression libraries (n = 63 to 81 lines per library) in both Indica (cv. IR64) and Japonica (cv. Cybonnet) backgrounds using three bio-geographically diverse wild donors representing the Oryza rufipogon Species Complex from China, Laos and Indonesia. The final libraries were genotyped using an Infinium 7K rice SNP array (C7AIR) and analyzed under greenhouse conditions for several simply inherited (Mendelian) traits. These six interspecific populations can be used as individual Chromosome Segment Substitution Line libraries and, when considered together, serve as a powerful genetic resource for systematic genetic dissection of agronomic, physiological and developmental traits in rice.

Genome Wide Association Analysis of Iron Overload in the Trans-Omics for Precision Medicine (TOPMed) Sickle Cell Disease Cohorts

Introduction: Transfusional iron (Fe) overload is a significant problem among patients with chronic, transfusion-dependent anemias. Iron overload is an important problem in pediatric sickle cell disease (SCD) patients on chronic transfusion regimens predominantly for primary and secondary prevention of stroke. Coexistent hereditary iron overload conditions contribute to the iron overload phenotype in SCD. For example, the Q248H mutation (rs11568350) in SLC40A1, which encodes ferroportin (FPN), is associated with a mild tendency to increase serum ferritin in the general population and with increased ferritin levels in SCD patients. Nevertheless, the molecular mechanisms underlying the progression to iron overload in SCD patients are poorly understood and more sensitive markers for outcome prediction that can be applied at early clinical stages are lacking. We hypothesize that genetic variation modifies the risk for iron overload in SCD patients and seek to validate previously identified mutation and identify novel genetic markers of iron overload among participants from TOPMed SCD cohorts by performing whole genome sequencing (WGS) association analyses. Methods: The WGS was performed by several national sequencing centers sponsored by NHLBI's TOPMed program at an average depth of 30× using DNA from SCD patient blood samples. Variant calling was performed jointly across TOPMed studies for all samples using the GotCloud pipeline by the TOPMed Informatics Research Center. The TOPMed data Coordinating Center performed quality control for sample identity. The data across the following studies were shared through the database of Genotype and Phenotype (dbGaP) exchange area: Howard PUSH SCD (N=370), OMG SCD (N=636), Walk PHaSST SCD (N=381) and REDS-III Brazil SCD (N=2620) with a total sample size of 4007. The study was approved by the appropriate institutional review boards (IRB) and informed consent was obtained from all participants. Genome Wide Association Analysis of iron overload was carried out using the University of Michigan ENCORE server. We performed single variant tests to test the association of log-transformed serum ferritin levels with single nucleotide variants (SNVs) while adjusting for sex, age, self-reported race, numbers of lifetime red blood transfusions and genetic substructure (PC's 1-10). We used a significance threshold of p<5×10-8 to report anassociation as genome-wide significant for common and rare genetic variants. Results: We first included PUSH SCD, OMG SCD and Walk PHaSST SCD cohorts with 840 serum ferritin samples in the WGS association analyses, which revealed at the genome-wide level a new rare variant rs137929759 (chr7:49538810 (GRCh38.p12), MAF=0.0043532, p=2.25×10−8). A few variants such as rs80097634 in gene AL163195.3 and RNASE11 (Chr14:20579417. MAF=0.050373, p=3.1×10-7) were close to the genome-wide significance level. We confirmed previously identified associations in SLC40A1 for ferritin (rs11568350, Chr2:189565370, MAF = 0.16853, p = 5.2×10−4). We also found several variants in AC105411.1, TJP1 and DCC that were close to genome-wide significance level. Further analysis will be carried out on the cloud-based platform provided by NHLBI BioData Catalyst using data from all the four cohorts to validate the previous analysis and expand to related phenotype such as transferrin, iron-overload status. Discussion: In this study we identified common and rare variants that associate with serum ferritin concentration. The results from this pilot study point to novel gene variants that might contribute to iron overload in SCD patients and serve as new biomarkers. Future analysis is needed to determine whether the identified variants can also help with therapeutics and outcome prediction for early stages of SCD-associated iron overload. Our findings will be useful for the future treatment of SCD patients and design of novel SCD therapeutics. ACKNOWLEDGMENTS: This work was supported by NIH Research Grants (1P50HL118006, and 1R01HL125005). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Disclosures Gordeuk: Imara: Research Funding; Global Blood Therapeutics: Consultancy, Research Funding; CSL Behring: Consultancy, Research Funding; Ironwood: Research Funding; Novartis: Consultancy. Telen:CSL Behring: Membership on an entity's Board of Directors or advisory committees, Research Funding; Forma Therapeutics: Research Funding; Novartis: Consultancy, Membership on an entity's Board of Directors or advisory committees; Pfizer: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; GlycoMimetics Inc.: Consultancy.

Climatic oscillations during the Mio/Pliocene epochs induced cladogenesis in the terrestrial snail genus Gittenedouardia (Mollusca: Gastropoda: Cerastidae) from South Africa

In South Africa, the terrestrial snail genus Gittenedouardia is the most species-rich member of the Cerastidae, where it is primarily distributed in the highly fragmented Afrotemperate and Indian Ocean coastal belt (IOCB) forest biomes. Phylogenetic relationships and cladogenetic events within the genus remain unstudied. In this respect, we reconstructed a dated phylogeny for eight Gittenedouardia species, and two populations identified to genus level using a combined mitochondrial (16S rRNA and COI) DNA sequencing dataset analysed using Bayesian inference and Maximum Likelihood framework. Furthermore, we investigated the population genetic substructure of the three widely distributed species (Gittenedouardia spadicea, G. natalensis and G. arenicola) for the COI locus, while also subsampling these species using the nuclear DNA ITS-2 locus. Phylogenetic results based on the combined mtDNA dataset supported the monophyly of Gittenedouardia and revealed three major clades and deep genetic structure among the three widely distributed species. Divergence-time estimates suggest that diversification within Gittenedouardia occurred during the middle Miocene/late Pliocene, a period characterised by a decrease in precipitation and the contraction of the Afrotemperate and IOCB forest biomes. We used two species delimitation methods, (PTP and STACEY) to infer putative species in G. spadicea, G. natalensis and G. arenicola. The two methods recovered a large number of evolutionary distinct units, with minimal consensus in the exact number of lineages. Our findings suggest the presence of undescribed diversity, necessitating the need for taxonomic revisionary work on Gittenedouardia. We discuss the climatic factors which may have contributed to the observed cladogenesis and compare our results with other studies of forest dwelling faunal taxa.

Signatures of selection in a recent invasion reveal adaptive divergence in a highly vagile invasive species.

A detailed understanding of population genetics in invasive populations helps us to identify drivers of successful alien introductions. Here, we investigate putative signals of selection in Australian populations of invasive common starlings, Sturnus vulgaris, and seek to understand how these have been influenced by introduction history. We used reduced representation sequencing to determine population structure, and identify Single Nucleotide Polymorphisms (SNPs) that are putatively under selection. We found that since their introduction into Australia, starling populations have become genetically differentiated despite the potential for high levels of dispersal, and that starlings have responded to selective pressures imposed by a wide range of environmental conditions across their geographic range. Isolation by distance appears to have played a strong role in determining genetic substructure across the starling's Australian range. Analyses of candidate SNPs that are putatively under selection indicated that aridity, precipitation and temperature may be important factors driving adaptive variation across the starling's invasive range in Australia. However, we also noted that the historic introduction regime may leave footprints on sites flagged as being under adaptive selection, and encourage critical interpretation of selection analyses in non-native populations.

A whole-genome sequenced control population in northern Sweden reveals subregional genetic differences

The number of national reference populations that are whole-genome sequenced are rapidly increasing. Partly driving this development is the fact that genetic disease studies benefit from knowing the genetic variation typical for the geographical area of interest. A whole-genome sequenced Swedish national reference population (n = 1000) has been recently published but with few samples from northern Sweden. In the present study we have whole-genome sequenced a control population (n = 300) (ACpop) from Västerbotten County, a sparsely populated region in northern Sweden previously shown to be genetically different from southern Sweden. The aggregated variant frequencies within ACpop are publicly available (DOI 10.17044/NBIS/G000005) to function as a basic resource in clinical genetics and for genetic studies. Our analysis of ACpop, representing approximately 0.11% of the population in Västerbotten, indicates the presence of a genetic substructure within the county. Furthermore, a demographic analysis showed that the population from which samples were drawn was to a large extent geographically stationary, a finding that was corroborated in the genetic analysis down to the level of municipalities. Including ACpop in the reference population when imputing unknown variants in a Västerbotten cohort resulted in a strong increase in the number of high-confidence imputed variants (up to 81% for variants with minor allele frequency < 5%). ACpop was initially designed for cancer disease studies, but the genetic structure within the cohort will be of general interest for all genetic disease studies in northern Sweden.

Sequence diversity of cytotoxic T cell antigens and satellite marker analysis of Theileria parva informs the immunization against East Coast fever in Rwanda

BackgroundEast Coast fever (ECF) caused by Theileria parva is endemic in Rwanda. In this study, the antigenic and genetic diversity of T. parva coupled with immunization and field challenge were undertaken to provide evidence for the introduction of ECF immunization in Rwanda.MethodsBlood collected from cattle in the field was screened for T. parva using ELISA and PCR targeting the p104 gene. Tp1 and Tp2 gene sequences were generated from field samples and from Gikongoro and Nyakizu isolates. Furthermore, multilocus genotype data was generated using 5 satellite markers and an immunization challenge trial under field conditions using Muguga cocktail vaccine undertaken.ResultsOut of 120 samples, 44 and 20 were positive on ELISA and PCR, respectively. Antigenic diversity of the Tp1 and Tp2 gene sequences revealed an abundance of Muguga, Kiambu and Serengeti epitopes in the samples. A further three clusters were observed on both Tp1 and Tp2 phylogenetic trees; two clusters comprising of field samples and vaccine isolates and the third cluster comprising exclusively of Rwanda samples. Both antigens exhibited purifying selection with no positive selection sites. In addition, satellite marker analysis revealed that field samples possessed both shared alleles with Muguga cocktail on all loci and also a higher proportion of unique alleles. The Muguga cocktail (Muguga, Kiambu and Serengeti) genotype compared to other vaccine isolates, was the most represented in the field samples. Further low genetic sub-structuring (FST = 0.037) coupled with linkage disequilibrium between Muguga cocktail and the field samples was observed. Using the above data to guide a field immunization challenge trial comprising 41 immunized and 40 control animals resulted in 85% seroconversion in the immunized animals and an efficacy of vaccination of 81.7%, implying high protection against ECF.ConclusionsAntigenic and genetic diversity analysis of T. parva facilitated the use of Muguga cocktail vaccine in field conditions. A protection level of 81.7% was achieved, demonstrating the importance of combining molecular tools with field trials to establish the suitability of implementation of immunization campaigns. Based on the information in this study, Muguga cocktail immunization in Rwanda has a potential to produce desirable results.

Contrasting the ancestry patterns of three distinct population groups from the northernmost region of South America.

Colombia, located in the north of the South American subcontinent is a country of great interest for population genetic studies given its high ethnic and cultural diversity represented by the admixed population, 102 indigenous peoples and African descent populations. In this study, an analysis of the genetic structure and ancestry was performed based on 46 ancestry informative INDEL markers (AIM-INDELs) and considering the genealogical and demographic variables of 451 unrelated individuals belonging to nine Native American, two African American, and four multiple ancestry populations. Measures of genetic diversity, ancestry components, and genetic substructure were analyzed to build a population model typical of the northernmost part of the South American continent. The model suggests three types of populations: Native American, African American, and multiple ancestry. The results support hypotheses posed by other authors about issues like the peopling of South America and the existence of two types of Native American ancestry. This last finding could be crucial for future research on the peopling of Colombia and South America in that a single origin of all indigenous communities should not be assumed. It then would be necessary to consider other events that could explain their genetic variability and complexity throughout the continent.

Genetic diversity analysis of a flax (Linum usitatissimum L.) global collection

BackgroundA sustainable breeding program requires a minimum level of germplasm diversity to provide varied options for the selection of new breeding lines. To maximize genetic gain of the North Dakota State University (NDSU) flax breeding program, we aimed to increase the genetic diversity of its parental stocks by incorporating diverse genotypes. For this purpose, we analyzed the genetic diversity, linkage disequilibrium, and population sub-structure of 350 globally-distributed flax genotypes with 6200 SNP markers.ResultsAll the genotypes tested clustered into seven sub-populations (P1 to P7) based on the admixture model and the output of neighbor-joining (NJ) tree analysis and principal coordinate analysis were in line with that of structure analysis. The largest sub-population separation arose from a cluster of NDSU/American genotypes with Turkish and Asian genotypes. All sub-populations showed moderate genetic diversity (average H = 0.22 and I = 0.34). The pairwise Fst comparison revealed a great degree of divergence (Fst > 0.25) between most of the combinations. A whole collection mantel test showed significant positive correlation (r = 0.30 and p < 0.01) between genetic and geographic distances, whereas it was non-significant for all sub-populations except P4 and P5 (r = 0.251, 0.349 respectively and p < 0.05). In the entire collection, the mean linkage disequilibrium was 0.03 and it decayed to its half maximum within < 21 kb distance.ConclusionsTo maximize genetic gain, hybridization between NDSU stock (P5) and Asian individuals (P6) are potentially the best option as genetic differentiation between them is highest (Fst > 0.50). In contrast, low genetic differentiation between P5 and P2 may enhance the accumulation of favorable alleles for oil and fiber upon crossing to develop dual purpose varieties. As each sub-population consists of many genotypes, a Neighbor-Joining tree and kinship matrix assist to identify distantly related genotypes. These results also inform genotyping decisions for future association mapping studies to ensure the identification of a sufficient number of molecular markers to tag all linkage blocks.

Accounting for cryptic population substructure enhances detection of inbreeding depression with genomic inbreeding coefficients: an example from a critically endangered marsupial.

Characterizing inbreeding depression in wildlife populations can be critical to their conservation. Coefficients of individual inbreeding can be estimated from genome-wide marker data. The degree to which sensitivity of inbreeding coefficients to population genetic substructure alters estimates of inbreeding depression in wild populations is not well understood. Using generalized linear models, we tested the power of two frequently used inbreeding coefficients that are calculated from genome-wide SNP markers, FH and F^III , to predict four fitness traits estimated over two decades in an isolated population of the critically endangered Leadbeater's possum. FH estimates inbreeding as excess observed homozygotes relative to equilibrium expectations, whereas F^III quantifies allelic similarity between the gametes that formed an individual, and upweights rare homozygotes. We estimated FH and F^III from 1,575 genome-wide SNP loci in individuals with fitness trait data (N=179-237 per trait), and computed revised coefficients, FHby group and F^IIIby group , adjusted for population genetic substructure by calculating them separately within two different genetic groups of individuals identified in the population. Using FH or F^III in the models, inbreeding depression was detected for survival to sexual maturity, longevity and whether individuals bred during their lifetime. F^IIIby group (but not FHby group ) additionally revealed significant inbreeding depression for lifetime reproductive output (total offspring assigned to each individual). Estimates of numbers of lethal equivalents indicated substantial inbreeding load, but differing between inbreeding estimators. Inbreeding depression, declining population size, and low and declining genetic diversity suggest that genetic rescue may assist in preventing extinction of this unique Leadbeater's possum population.

The immunogenetic diversity of the HLA system in Mexico correlates with underlying population genetic structure

We studied HLA class I (HLA-A, -B) and class II (HLA-DRB1, -DQB1) allele groups and alleles by PCR-SSP based typing in a total of 15,318 mixed ancestry Mexicans from all the states of the country divided into 78 sample sets, providing information regarding allelic and haplotypic frequencies and their linkage disequilibrium, as well as admixture estimates and genetic substructure. We identified the presence of 4268 unique HLA extended haplotypes across Mexico and find that the ten most frequent (HF > 1%) HLA haplotypes with significant linkage disequilibrium (Δ’≥0.1) in Mexico (accounting for 20% of the haplotypic diversity of the country) are of primarily Native American ancestry (A*02~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*08~DQB1*04, A*68~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*14~DQB1*03:01, A*24~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*04~DQB1*03:02, A*02~B*40:02~DRB1*04~DQB1*03:02, A*68~B*35~DRB1*04~DQB1*03:02, A*02~B*15:01~DRB1*04~DQB1*03:02). Admixture estimates obtained by a maximum likelihood method using HLA-A/-B/-DRB1 as genetic estimators revealed that the main genetic components in Mexico as a whole are Native American (ranging from 37.8% in the northern part of the country to 81.5% in the southeastern region) and European (ranging from 11.5% in the southeast to 62.6% in northern Mexico). African admixture ranged from 0.0 to 12.7% not following any specific pattern. We were able to detect three major immunogenetic clusters correlating with genetic diversity and differential admixture within Mexico: North, Central and Southeast, which is in accordance with previous reports using genome-wide data. Our findings provide insights into the population immunogenetic substructure of the whole country and add to the knowledge of mixed ancestry Latin American population genetics, important for disease association studies, detection of demographic signatures on population variation and improved allocation of public health resources.

Genetic Characterization of the "Chusca Lojana", a Creole Goat Reared in Ecuador, and Its Relationship with Other Goat Breeds.

Simple SummaryAn individual from a population presents a series of characteristics that differ from the rest and that increase as the kinship relationships are lower; this leads to the fact that two populations that stop exchanging genetic material through mating eventually come to present characters common to all members of each, but different between the two. This was what happened with the animal populations brought to America more than 500 years ago from the Iberian Peninsula in the colonization period, resulting in Creole populations that inhabit the harshest environments of our immense geography. The Creole goat “Chusca Lojana” has adapted to live in the dry forest region of Southern Ecuador, where environmental conditions are warm-dry, with sparse vegetation and a rather irregular topography. In the present study, the intra-breed genetic diversity of this goat is analyzed as well as its genetic relationships with other breeds. Significant FIS and intra-breed structure show that there is some heterogeneity and structure within the breed. However, inter-breed structure results underline that this breed is differentiated from other Creole breeds, because crossbreeding with other breeds was not detected; therefore, we must take advantage of this valuable genetic resource, and ensure its conservation and selection.The largest population of goats (62%) in Ecuador is in the dry forest region in the south of the country. A Creole goat, named “Chusca Lojana”, has adapted to the dry forest region where environmental conditions are warm-dry, with sparse vegetation. Knowledge of the genetic information of the Creole goat is important to determine intra-racial diversity, the degree of genetic distance among other breeds of goats, and the possible substructure of the population, which is valuable for the conservation of such a species’ genetic resources. A total of 145 samples of the Creole goat was taken from the four biotypes previously identified. Genetic analyses were performed using 38 microsatellites recommended for studies of goat genetic diversity (FAO-ISAG). The results of within-breed genetic diversity showed a mean number of alleles per locus (MNA) of 8, an effective number of alleles (Ae) of 4.3, an expected heterozygosity (He) of 0.71, an observed heterozygosity (Ho) of 0.63, polymorphic information content (PIC) of 0.67, and an FIS value of 0.11. Between-breed genetic diversity among 43 goat populations (native of Spain, American Creole, Europeans, and Africans) showed the following values: FIS = 0.087, FIT = 0.176, and FST = 0.098. Regarding the analysis of the population structure, the results showed that the Creole Chusca Lojana goat population is homogeneous and no genetic separation was observed between the different biotypes (FST = 0.0073). In conclusion, the Chusca Lojana goat has a high genetic diversity, without exhibiting a genetic substructure. Therefore, it should be considered as a distinct population because crossbreeding with other breeds was not detected.

Ancient genome-wide DNA from France highlights the complexity of interactions between Mesolithic hunter-gatherers and Neolithic farmers.

Starting from 12,000 years ago in the Middle East, the Neolithic lifestyle spread across Europe via separate continental and Mediterranean routes. Genomes from early European farmers have shown a clear Near Eastern/Anatolian genetic affinity with limited contribution from hunter-gatherers. However, no genomic data are available from modern-day France, where both routes converged, as evidenced by a mosaic cultural pattern. Here, we present genome-wide data from 101 individuals from 12 sites covering today's France and Germany from the Mesolithic (N = 3) to the Neolithic (N = 98) (7000-3000 BCE). Using the genetic substructure observed in European hunter-gatherers, we characterize diverse patterns of admixture in different regions, consistent with both routes of expansion. Early western European farmers show a higher proportion of distinctly western hunter-gatherer ancestry compared to central/southeastern farmers. Our data highlight the complexity of the biological interactions during the Neolithic expansion by revealing major regional variations.

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.

Genomic studies on Strongyloides stercoralis in northern and western Thailand

BackgroundStrongyloidiasis is a soil borne helminthiasis, which in most cases is caused by Strongyloides stercoralis. Human infections with S. fuelleborni fuelleborni and S. fuelleborni kellyi also occur. Although up to 370 million people are currently estimated to be infected with S. stercoralis, this parasite is frequently overlooked. Strongyloides stercoralis is prevalent among humans in Thailand; however, S. fuelleborni fuelleborni has also been reported. Three recent genomic studies of individual S. stercoralis worms found genetically diverse populations of S. stercoralis, with comparably low heterozygosity in Cambodia and Myanmar, and less diverse populations with high heterozygosity in Japan and southern China that presumably reproduce asexually.MethodsWe isolated individual Strongyloides spp. from different localities in northern and western Thailand and determined their nuclear small ribosomal subunit rDNA (18S rDNA, SSU), in particular the hypervariable regions I and IV (HVR-I and HVR-IV), mitochondrial cytochrome c oxidase subunit 1 (cox1) and for a subset whole genome sequences. These sequences were then compared with each other and with published sequences from different geographical locations.ResultsAll 237 worms isolated from 16 different human hosts were S. stercoralis, no S. fuelleborni was found. All worms had the common S. stercoralis SSU HVR IV haplotype A. Two different SSU HVR I haplotypes (I and II), both previously described in S. stercoralis, were found. No animal heterozygous for the two haplotypes was identified. Among the twelve cox1 haplotypes found, five had not been previously described. Based upon the mitochondrial cox1 and the nuclear whole genome sequences, S. stercoralis in Thailand was phylogenetically intermixed with the samples from other Southeast Asian countries and did not form its own branch. The genomic heterozygosity was even slightly lower than in the samples from the neighboring countries.ConclusionsIn our sample from humans, all Strongyloides spp. were S. stercoralis. The S. stercoralis from northern and western Thailand appear to be part of a diverse, intermixing continental Southeast Asian population. No obvious indication for genetic sub-structuring of S. stercoralis within Thailand or within the Southeast Asian peninsula was detected.

Reconstructing early routes of invasion of the bronze bug Thaumastocoris peregrinus (Hemiptera: Thaumastocoridae): cities as bridgeheads for global pest invasions

Insect pest invasions pose a major threat to agriculture, forestry and many natural ecosystems. Thaumastocoris peregrinus is an invasive sap-sucking pest of significant economic importance to Eucalyptus forestry that has recently invaded several countries worldwide. In this study we identify the origin and retrace the invasion history of T. peregrinus. We analysed samples from six locations in Africa, South America and Australia using microsatellites markers and a combination of clustering methods and scenario testing using Approximate Bayesian Clustering. We detected clear genetic substructure differentiating African and South American samples, with representatives of both present in Australia. The Australian population from New South Wales showed substantially higher genetic diversity than the Queensland source, which could indicate that this region could be part of the core range and evolutionary origin of the species. Africa and South America were colonised by independent introductions that occurred more or less concurrently. The study illustrates the impact of the bridgehead effect on global invasions following an outbreak or ‘invasion’ within a city in the native range of the insect.

Pliocene-Early Pleistocene Geological Events Structure Pacific Martens (Martes caurina).

The complex topography, climate, and geological history of Western North America have shaped contemporary patterns of biodiversity and species distributions in the region. Pacific martens (Martes caurina) are distributed along the northern Pacific Coast of North America with disjunct populations found throughout the Northwestern Forested Mountains and Marine West Coast Forest ecoregions of the West Coast. Martes in this region have been classified into subspecies; however, the subspecific designation has been extensively debated. In this study, we use genomic data to delineate conservation units of Pacific marten in the Sierra-Cascade-Coastal montane belt in the western United States. We analyzed the mitochondrial genome for 94 individuals to evaluate the spatial distribution and divergence times of major lineages. We further genotyped 401 individuals at 13 microsatellite loci to investigate major patterns of population structure. Both nuclear and mitochondrial DNA suggest substantial genetic substructure concordant with historical subspecies designations. Our results revealed that the region contains 2 distinct mitochondrial lineages: a Cascades/Sierra lineage that diverged from the Cascades/coastal lineage 2.23 (1.48-3.14 mya), consistent with orogeny of the Cascade Mountain chain. Interestingly, Pacific Martes share phylogeographic patterns similar with other sympatric taxa, suggesting that the complex geological history has shaped the biota of this region. The information is critical for conservation and management efforts, and further investigation of adaptive diversity is warranted following appropriate revision of conservation management designations.

Population genetic analysis of Shaanxi male Han Chinese population reveals genetic differentiation and homogenization of East Asians.

BackgroundShaanxi province, located in the upper Yellow River, has been evidenced as the geographic origin of Chinese civilization, Sino‐Tibetan‐speaking language, and foxtail or broomcorn millet farmers via the linguistic phylogenetic spectrum, archeological documents, and genetic evidence. Nowadays, Han Chinese is the dominant population in this area. The formation process of modern Shaanxi Han population reconstructed via the ancient DNA is on the way, however, the patterns of genetic relationships of modern Shaanxi Han, allele frequency distributions of high mutated short tandem repeats (STRs) and corresponding forensic parameters are remained to be explored.MethodsHere, we successfully genotyped 23 autosomal STRs in 630 unrelated Shaanxi male Han individuals using the recently updated Huaxia Platinum PCR amplification system. Forensic allele frequency and parameters of all autosomal STRs were assessed. And comprehensive population genetic structure was explored via various typical statistical technologies.ResultsPopulation genetic analysis based on the raw‐genotype dataset among 15,803 Eurasian individuals and frequency datasets among 56 populations generally illustrated that linguistic stratification is significantly associated with the genetic substructure of the East Asian population. Principal component analysis, multidimensional scaling plots and phylogenetic tree further demonstrated that Shaanxi Han has a close genetic relationship with geographically close Shanxi Han, and showed that Han Chinese is a homogeneous population during the historic and recent admixture from the STR variations. Except for Sinitic‐speaking populations, Shaanxi Han harbored more alleles sharing with Tibeto‐Burman‐speaking populations than with other reference populations. Focused on the allele frequency correlation and forensic parameters, all loci are in accordance with the minimum requirements of HWE and LD. The observed combined probability of discrimination of 8.2201E‐28 and the cumulative power of exclusion of 0.9999999995 in Shaanxi Han demonstrated that the studied STR loci are informative and polymorphic, and this system can be used as a powerful routine forensic tool in personal identification and parentage testing.ConclusionBoth the geographical and linguistic divisions have shaped the genetic structure of modern East Asian. And more forensic reference data should be obtained for ethnically, culturally, geographically and linguistically different populations for better routine forensic practice and population genetic studies.