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.

Massively parallel sequencing-enabled mixture analysis of mitochondrial DNA samples.

The mitochondrial genome has a number of characteristics that provide useful information to forensic investigations. Massively parallel sequencing (MPS) technologies offer improvements to the quantitative analysis of the mitochondrial genome, specifically the interpretation of mixed mitochondrial samples. Two-person mixtures with nuclear DNA ratios of 1:1, 5:1, 10:1, and 20:1 of individuals from different and similar phylogenetic backgrounds and three-person mixtures with nuclear DNA ratios of 1:1:1 and 5:1:1 were prepared using the Precision ID mtDNA Whole Genome Panel and Ion Chef, and sequenced on the Ion PGM or Ion S5 sequencer (Thermo Fisher Scientific, Waltham, MA, USA). These data were used to evaluate whether and to what degree MPS mixtures could be deconvolved. Analysis was effective in identifying the major contributor in each instance, while SNPs from the minor contributor's haplotype only were identified in the 1:1, 5:1, and 10:1 two-person mixtures. While the major contributor was identified from the 5:1:1 mixture, analysis of the three-person mixtures was more complex, and the mixed haplotypes could not be completely parsed. These results indicate that mixed mitochondrial DNA samples may be interpreted with the use of MPS technologies.

Identification and analysis of mtDNA genomes attributed to Finns reveal long-stagnant demographic trends obscured in the total diversity.

In Europe, modern mitochondrial diversity is relatively homogeneous and suggests an ubiquitous rapid population growth since the Neolithic revolution. Similar patterns also have been observed in mitochondrial control region data in Finland, which contrasts with the distinctive autosomal and Y-chromosomal diversity among Finns. A different picture emerges from the 843 whole mitochondrial genomes from modern Finns analyzed here. Up to one third of the subhaplogroups can be considered as Finn-characteristic, i.e. rather common in Finland but virtually absent or rare elsewhere in Europe. Bayesian phylogenetic analyses suggest that most of these attributed Finnish lineages date back to around 3,000-5,000 years, coinciding with the arrival of Corded Ware culture and agriculture into Finland. Bayesian estimation of past effective population sizes reveals two differing demographic histories: 1) the 'local' Finnish mtDNA haplotypes yielding small and dwindling size estimates for most of the past; and 2) the 'immigrant' haplotypes showing growth typical of most European populations. The results based on the local diversity are more in line with that known about Finns from other studies, e.g., Y-chromosome analyses and archaeology findings. The mitochondrial gene pool thus may contain signals of local population history that cannot be readily deduced from the total diversity.

More comprehensive forensic genetic marker analyses for accurate human remains identification using massively parallel DNA sequencing.

BACKGROUND: Although the primary objective of forensic DNA analyses of unidentified human remains is positive identification, cases involving historical or archaeological skeletal remains often lack reference samples for comparison. Massively parallel sequencing (MPS) offers an opportunity to provide biometric data in such cases, and these cases provide valuable data on the feasibility of applying MPS for characterization of modern forensic casework samples. In this study, MPS was used to characterize 140-year-old human skeletal remains discovered at a historical site in Deadwood, South Dakota, United States. The remains were in an unmarked grave and there were no records or other metadata available regarding the identity of the individual. Due to the high throughput of MPS, a variety of biometric markers could be typed using a single sample. RESULTS: Using MPS and suitable forensic genetic markers, more relevant information could be obtained from a limited quantity and quality sample. Results were obtained for 25/26 Y-STRs, 34/34 Y SNPs, 166/166 ancestry-informative SNPs, 24/24 phenotype-informative SNPs, 102/102 human identity SNPs, 27/29 autosomal STRs (plus amelogenin), and 4/8 X-STRs (as well as ten regions of mtDNA). The Y-chromosome (Y-STR, Y-SNP) and mtDNA profiles of the unidentified skeletal remains are consistent with the R1b and H1 haplogroups, respectively. Both of these haplogroups are the most common haplogroups in Western Europe. Ancestry-informative SNP analysis also supported European ancestry. The genetic results are consistent with anthropological findings that the remains belong to a male of European ancestry (Caucasian). Phenotype-informative SNP data provided strong support that the individual had light red hair and brown eyes. CONCLUSIONS: This study is among the first to genetically characterize historical human remains with forensic genetic marker kits specifically designed for MPS. The outcome demonstrates that substantially more genetic information can be obtained from the same initial quantities of DNA as that of current CE-based analyses.

Whole mitochondrial genome genetic diversity in an Estonian population sample.

Mitochondrial DNA is a useful marker for population studies, human identification, and forensic analysis. Commonly used hypervariable regions I and II (HVI/HVII) were reported to contain as little as 25% of mitochondrial DNA variants and therefore the majority of power of discrimination of mitochondrial DNA resides in the coding region. Massively parallel sequencing technology enables entire mitochondrial genome sequencing. In this study, buccal swabs were collected from 114 unrelated Estonians and whole mitochondrial genome sequences were generated using the Illumina MiSeq system. The results are concordant with previous mtDNA control region reports of high haplogroup HV and U frequencies (47.4 and 23.7% in this study, respectively) in the Estonian population. One sample with the Northern Asian haplogroup D was detected. The genetic diversity of the Estonian population sample was estimated to be 99.67 and 95.85%, for mtGenome and HVI/HVII data, respectively. The random match probability for mtGenome data was 1.20 versus 4.99% for HVI/HVII. The nucleotide mean pairwise difference was 27 ± 11 for mtGenome and 7 ± 3 for HVI/HVII data. These data describe the genetic diversity of the Estonian population sample and emphasize the power of discrimination of the entire mitochondrial genome over the hypervariable regions.

High sensitivity multiplex short tandem repeat loci analyses with massively parallel sequencing.

STR typing in forensic genetics has been performed traditionally using capillary electrophoresis (CE). However, CE-based method has some limitations: a small number of STR loci can be used; stutter products, dye artifacts and low level alleles. Massively parallel sequencing (MPS) has been considered a viable technology in recent years allowing high-throughput coverage at a relatively affordable price. Some of the CE-based limitations may be overcome with the application of MPS. In this study, a prototype multiplex STR System (Promega) was amplified and prepared using the TruSeq DNA LT Sample Preparation Kit (Illumina) in 24 samples. Results showed that the MinElute PCR Purification Kit (Qiagen) was a better size selection method compared with recommended diluted bead mixtures. The library input sensitivity study showed that a wide range of amplicon product (6-200ng) could be used for library preparation without apparent differences in the STR profile. PCR sensitivity study indicated that 62pg may be minimum input amount for generating complete profiles. Reliability study results on 24 different individuals showed that high depth of coverage (DoC) and balanced heterozygote allele coverage ratios (ACRs) could be obtained with 250pg of input DNA, and 62pg could generate complete or nearly complete profiles. These studies indicate that this STR multiplex system and the Illumina MiSeq can generate reliable STR profiles at a sensitivity level that competes with current widely used CE-based method.

Evaluation of a novel material, Diomics X-Swab™, for collection of DNA.

Success of DNA typing is related to the amount of target material recovered from an evidentiary item. Generally, the more DNA that is recovered, the better the chance is of obtaining a typing result that will be robust and reliable. One method of collecting stain materials is by swabbing. Recovery of DNA from a number of commercially available swabs is not an efficient process. The X-Swab™ (Diomics Corporation, La Jolla, CA) is a unique bio-specimen collection material with highly absorptive properties and can be dissolved during certain extraction conditions. Therefore, more DNA may be collected from a substrate and be released from the swab matrix than other swabs. The ability to recover DNA from X-Swab material and success in STR typing were compared with the Copan 4N6FLOQSwab™ (Brescia, Italy), a device which utilizes a proprietary flocked-swab technology to maximize DNA collection and elution efficiency. Both types of swabs were impregnated with known amounts of DNA and body fluids and allowed to air dry. In addition, blood was placed onto glass slides, allowed to dry and collected using both types of swabs. DNA recovery was assessed by DNA quantitation and by STR typing. Results suggested that X-Swab material yielded greater DNA recovery, particularly of low quantity samples (defined as diluted neat samples), compared with the 4N6FLOQSwab. Results also indicated that X-Swab material itself enhances yield of PCR products.

A high volume extraction and purification method for recovering DNA from human bone.

DNA recovery, purity and overall extraction efficiency of a protocol employing a novel silica-based column, Hi-Flow(®) (Generon Ltd., Maidenhead, UK), were compared with that of a standard organic DNA extraction methodology. The quantities of DNA recovered by each method were compared by real-time PCR and quality of DNA by STR typing using the PowerPlex(®) ESI 17 Pro System (Promega Corporation, Madison, WI) on DNA from 10 human bone samples. Overall, the Hi-Flow method recovered comparable quantities of DNA ranging from 0.8ng±1 to 900ng±159 of DNA compared with the organic method ranging from 0.5ng±0.9 to 855ng±156 of DNA. Complete profiles (17/17 loci tested) were obtained for at least one of three replicates for 3/10 samples using the Hi-Flow method and from 2/10 samples with the organic method. All remaining bone samples yielded partial profiles for all replicates with both methods. Compared with a standard organic DNA isolation method, the results indicated that the Hi-Flow method provided equal or improved recovery and quality of DNA without the harmful effects of organic extraction. Moreover, larger extraction volumes (up to 20mL) can be employed with the Hi-Flow method which enabled more bone sample to be extracted at one time.

High-quality and high-throughput massively parallel sequencing of the human mitochondrial genome using the Illumina MiSeq.

Mitochondrial DNA typing in forensic genetics has been performed traditionally using Sanger-type sequencing. Consequently sequencing of a relatively-large target such as the mitochondrial genome (mtGenome) is laborious and time consuming. Thus, sequencing typically focuses on the control region due to its high concentration of variation. Massively parallel sequencing (MPS) has become more accessible in recent years allowing for high-throughput processing of large target areas. In this study, Nextera(®) XT DNA Sample Preparation Kit and the Illumina MiSeq™ were utilized to generate quality whole genome mitochondrial haplotypes from 283 individuals in a both cost-effective and rapid manner. Results showed that haplotypes can be generated at a high depth of coverage with limited strand bias. The distribution of variants across the mitochondrial genome was described and demonstrated greater variation within the coding region than the non-coding region. Haplotype and haplogroup diversity were described with respect to whole mtGenome and HVI/HVII. An overall increase in haplotype or genetic diversity and random match probability, as well as better haplogroup assignment demonstrates that MPS of the mtGenome using the Illumina MiSeq system is a viable and reliable methodology.