Feeding up to 91% concentrate to Holstein and Jersey dairy cows: Effects on enteric methane emission, rumen fermentation and bacterial community, digestibility, production, and feeding behavior.

Due to climate change, periods of drought might be longer and occur more frequently, which challenges roughage production and requires changed feeding of dairy cattle by increasing the grain content of the diet. This study investigated the effect of diets with concentrate proportions up to 91% of dry matter on dry matter intake (DMI), milk production, enteric methane emission, rumen fermentation, rumen bacterial community structure, nutrient digestibility, and feeding behavior of Holstein and Jersey dairy cows. Twelve Danish Holstein and 12 Danish Jersey cows were fed ad libitum with one of 3 total mixed rations differing in concentrate proportion in a continuous design with staggered approach over 19 to 29 d. Dietary concentrate proportions were 49% (C49), 70% (C70), and 91% (C91) on dry matter basis, and were based on increasing proportions of chopped barley straw, dried beet pulp, barley, NaOH-treated wheat, dried distillers grain, and rapeseed cake at the expense of grass/clover silage, corn silage and soybean meal. Cows were adapted to the diets over a 12- to 19-d period, before rumination activity was measured over 3 d. Subsequently, spot samples of feces were collected for digestibility determination over 2 d, and gas exchange was measured on the last 3 d of the experimental period. Shortly after chamber stay, rumen liquid was collected using an oro-ruminal device. Dry matter intake was higher for Holstein than Jersey. Methane emissions (all expressions) were affected by the interaction between breed and diet. Methane per kilogram of DMI was lowered by 18 and 48% for Holstein fed C70 and C91, respectively, compared with C49, whereas this was 17 and 22% respectively for Jersey. Rumen propionate molar proportion increased more, rumen bacterial community was less diverse, and rumination time and rumination chews relative to DMI reduced less for Holstein than for Jersey cows with increasing concentrate level. In conclusion, Holstein dairy cows responded stronger to increased dietary concentrate level regarding methane mitigation, changes in rumen VFA profile, and effect on the rumen bacterial community structure than Jersey cows, whereas Jersey cows responded stronger with regard to rumination time and rumination chews (per kilogram of DMI and per kilogram of neutral detergent fiber intake) than Holstein cows. Thus, diets high in concentrates are a less effective methane mitigation strategy for Jersey than for Holstein.

Evaluation of a custom QIAseq targeted DNA panel with 164 ancestry informative markers sequenced with the Illumina MiSeq.

Introduction of new methods requires meticulous evaluation before they can be applied to forensic genetic case work. Here, a custom QIAseq Targeted DNA panel with 164 ancestry informative markers was assessed using the MiSeq sequencing platform. Concordance, sensitivity, and the capability for analysis of mixtures were tested. The assay gave reproducible and nearly concordant results with an input of 10 and 2 ng DNA. Lower DNA input led to an increase in both locus and allele drop-outs, and a higher variation in heterozygote balance. Locus or allele drop-outs in the samples with less than 2 ng DNA input were not necessarily associated with the overall performance of a locus. Thus, the QIAseq assay will be difficult to implement in a forensic genetic setting where the sample material is often scarce and of poor quality. With equal or near equal mixture ratios, the mixture DNA profiles were easily identified by an increased number of imbalanced heterozygotes. For more skewed mixture ratios, the mixture DNA profiles were identified by an increased noise level. Lastly, individuals from Great Britain and the Middle East were investigated. The Middle Eastern individuals showed a greater affinity with South European populations compared to North European populations.

Assessment of the effectiveness of the EUROFORGEN NAME and Precision ID Ancestry panel markers for ancestry investigations.

The EUROFORGEN NAME panel is a regional ancestry panel designed to differentiate individuals from the Middle East, North Africa, and Europe. The first version of the panel was developed for the MassARRAY system and included 111 SNPs. Here, a custom AmpliSeq EUROFORGEN NAME panel with 102 of the original 111 loci was used to sequence 1098 individuals from 14 populations from Europe, the Middle East, North Africa, North-East Africa, and South-Central Asia. These samples were also sequenced with a global ancestry panel, the Precision ID Ancestry Panel. The GenoGeographer software was used to assign the AIM profiles to reference populations and calculate the weight of the evidence as likelihood ratios. The combination of the EUROFORGEN NAME and Precision ID Ancestry panels led to fewer ambiguous assignments, especially for individuals from the Middle East and South-Central Asia. The likelihood ratios showed that North African individuals could be separated from European and Middle Eastern individuals using the Precision ID Ancestry Panel. The separation improved with the addition of the EUROFORGEN NAME panel. The analyses also showed that the separation of Middle Eastern populations from European and South-Central Asian populations was challenging even when both panels were applied.

Enteric methane emission and digestion in dairy cows fed wheat or molasses.

The aim of this experiment was to measure enteric methane (CH4) emission and its relation with rumen digestion in dairy cows fed diets rich in 1 of the 2 carbohydrate sources, starch or sugar. The rations were based on late first-cut grass-clover silage supplemented with wheat (Wh), NaOH-treated wheat (Wh+NaOH), sugar beet molasses (Mo), or sugar beet molasses with addition of sodium bicarbonate (Mo+Bic). Wheat and molasses made up 35% of dry matter in the 2 diets with molasses and wheat, respectively. Four cows fitted with ruminal, duodenal, and ileal canulae were used in a 4 × 4 Latin square design. Nutrient digestibility was measured using chromium oxide and titanium oxide as flow markers, and emissions of CH4 and hydrogen were measured via open-circuit indirect calorimetry on 4 consecutive days. Data were analyzed using PROC MIXED of SAS (version 9.4; SAS Institute Inc., Cary, NC) with treatment and period as fixed effects and cow as random effect. Furthermore, orthogonal contrasts were calculated. The cows produced 32.5, 33.6, 36.2, and 35.1 L of CH4/kg of dry matter intake (DMI) on diets Wh, Wh+NaOH, Mo, and Mo+Bic, respectively. The emission of CH4 per day, per kilogram of DMI, and per kilogram of energy-corrected milk as well as daily hydrogen emission were higher on the Mo diet compared with the Wh diet. With the present inclusion of wheat and molasses in the diet, no effects of NaOH treatment of wheat or of sodium bicarbonate supplementation to the Mo diet could be demonstrated on CH4 emission expressed per kilogram of DMI or per kilogram of energy-corrected milk. The duodenal flow of starch was higher when wheat was treated with NaOH. Under the conditions in the present experiment, ruminal NDF digestibility was not affected by carbohydrate source, NaOH treatment of wheat, or bicarbonate supplementation. Total volatile fatty acid concentration in the rumen and the proportions of acetate and propionate were not affected by carbohydrate source, NaOH treatment of wheat, or bicarbonate supplementation. Likewise, we could not show any influence of diet on microbial protein synthesis or efficiency of microbial protein synthesis expressed as grams of microbial protein synthesis per kilogram of true rumen-digested organic matter. We concluded that CH4 emission was increased when wheat was replaced by molasses, whereas no effect of manipulating rumen fermentation by NaOH treatment of wheat or addition of bicarbonate to molasses could be found with a level of approximately 25% of dry matter from starch and sugar, respectively.

Analysis of Uyghur and Kazakh populations using the Precision ID Ancestry Panel.

Autosomal ancestry informative markers (AIMs) are important markers for inferring ancestry of humans. In the present study, we typed 105 Uyghurs and 94 Kazakhs with the Precision ID Ancestry Panel that amplifies 165 autosomal AIMs. No statistically significant deviation from Hardy-Weinberg equilibrium and no linkage disequilibrium between loci was observed after Bonferroni correction. STRUCTURE and PCA analyses showed that Uyghurs and Kazakhs appeared as admixed individuals of primarily European and East Asian ancestry and were clearly differentiated from Europeans, Middle Easterners, South/Central Asians, and East Asians. However, it was not possible to differentiate the two populations from each other and they were also difficult to differentiate from Greenlanders, a population with European/Inuit admixture. GenoGeographer was used to evaluate the weight of the evidence. Initially, the results showed that the majority of AIM profiles from Uyghur and Kazakh individuals were not represented by any of the 36 reference populations of the GenoGeographer database. Consequently, it was not reasonable to infer the ancestry of these individuals. A randomly selected subset of the studied populations (75 Uyghur and 75 Kazakh individuals) was used to construct two new reference populations for GenoGeographer, and ancestry prediction was performed on the remaining test individuals. A total of 42 out of 49 test individuals were represented by at least one population after the introduction of Uyghur and Kazakh reference populations. Likelihood ratios ≥106 were obtained when the alternative hypothesis was that the individual belonged to the South/Central Asian, East Asian, Middle Eastern, European, or the admixed Greenlandic population.

Development and validation of the EUROFORGEN NAME (North African and Middle Eastern) ancestry panel.

Inference of biogeographic origin is an important factor in clinical, population and forensic genetics. The information provided by AIMs (Ancestry Informative Markers) can allow the differentiation of major continental population groups, and several AIM panels have been developed for this purpose. However, from these major population groups, Eurasia covers a wide area between two continents that is difficult to differentiate genetically. These populations display a gradual genetic cline from West Europe to South Asia in terms of allele frequency distribution. Although differences have been reported between Europe and South Asia, Middle East populations continue to be a target of further investigations due to the lack of genetic variability, therefore hampering their genetic differentiation from neighboring populations. In the present study, a custom-built ancestry panel was developed to analyze North African and Middle Eastern populations, designated the 'NAME' panel. The NAME panel contains 111 SNPs that have patterns of allele frequency differentiation that can distinguish individuals originating in North Africa and the Middle East when combined with a previous set of 126 Global AIM-SNPs.

The Danish STR sequence database: duplicate typing of 363 Danes with the ForenSeq™ DNA Signature Prep Kit.

Some STR loci have internal sequence variations, which are not revealed by the standard STR typing methods used in forensic genetics (PCR and fragment length analysis by capillary electrophoresis (CE)). Typing of STRs with next-generation sequencing (NGS) uncovers the sequence variation in the repeat region and in the flanking regions. In this study, 363 Danish individuals were typed for 56 STRs (26 autosomal STRs, 24 Y-STRs, and 6 X-STRs) using the ForenSeq™ DNA Signature Prep Kit to establish a Danish STR sequence database. Increased allelic diversity was observed in 34 STRs by the PCR-NGS assay. The largest increases were found in DYS389II and D12S391, where the numbers of sequenced alleles were around four times larger than the numbers of alleles determined by repeat length alone. Thirteen SNPs and one InDel were identified in the flanking regions of 12 STRs. Furthermore, 36 single positions and five longer stretches in the STR flanking regions were found to have dubious genotyping quality. The combined match probability of the 26 autosomal STRs was 10,000 times larger using the PCR-NGS assay than by using PCR-CE. The typical paternity indices for trios and duos were 500 and 100 times larger, respectively, than those obtained with PCR-CE. The assay also amplified 94 SNPs selected for human identification. Eleven of these loci were not in Hardy-Weinberg equilibrium in the Danish population, most likely because the minimum threshold for allele calling (30 reads) in the ForenSeq™ Universal Analysis Software was too low and frequent allele dropouts were not detected.

Forensic SNP genotyping with SNaPshot: Technical considerations for the development and optimization of multiplexed SNP assays.

This review explores the key factors that influence the optimization, routine use, and profile interpretation of the SNaPshot single-base extension (SBE) system applied to forensic single-nucleotide polymorphism (SNP) genotyping. Despite being a mainly complimentary DNA genotyping technique to routine STR profiling, use of SNaPshot is an important part of the development of SNP sets for a wide range of forensic applications with these markers, from genotyping highly degraded DNA with very short amplicons to the introduction of SNPs to ascertain the ancestry and physical characteristics of an unidentified contact trace donor. However, this technology, as resourceful as it is, displays several features that depart from the usual STR genotyping far enough to demand a certain degree of expertise from the forensic analyst before tackling the complex casework on which SNaPshot application provides an advantage. In order to provide the basis for developing such expertise, we cover in this paper the most challenging aspects of the SNaPshot technology, focusing on the steps taken to design primer sets, optimize the PCR and single-base extension chemistries, and the important features of the peak patterns observed in typical forensic SNP profiles using SNaPshot. With that purpose in mind, we provide guidelines and troubleshooting for multiplex-SNaPshot-oriented primer design and the resulting capillary electrophoresis (CE) profile interpretation (covering the most commonly observed artifacts and expected departures from the ideal conditions).

A collaborative EDNAP exercise on SNaPshot™-based mtDNA control region typing.

A collaborative European DNA Profiling (EDNAP) Group exercise was undertaken to assess the performance of an earlier described SNaPshot™-based screening assay (denoted mini-mtSNaPshot) (Weiler et al., 2016) [1] that targets 18 single nucleotide polymorphism (SNP) positions in the mitochondrial (mt) DNA control region and allows for discrimination of major European mtDNA haplogroups. Besides the organising laboratory, 14 forensic genetics laboratories were involved in the analysis of 13 samples, which were centrally prepared and thoroughly tested prior to shipment. The samples had a variable complexity and comprised straightforward single-source samples, samples with dropout or altered peak sizing, a point heteroplasmy and two-component mixtures resulting in one to five bi-allelic calls. The overall success rate in obtaining useful results was high (97.6%) given that some of the participating laboratories had no previous experience with the typing technology and/or mtDNA analysis. The majority of the participants proceeded to haplotype inference to assess the feasibility of assigning a haplogroup and checking phylogenetic consistency when only 18 SNPs are typed. To mimic casework procedures, the participants compared the SNP typing data of all 13 samples to a set of eight mtDNA reference profiles that were described according to standard nomenclature (Parson et al., 2014) [2], and indicated whether these references matched each sample or not. Incorrect scorings were obtained for 2% of the comparisons and derived from a subset of the participants, indicating a need for training and guidelines regarding mini-mtSNaPshot data interpretation.

Thirty autosomal insertion-deletion polymorphisms analyzed using the Investigator® DIPplex Kit in populations from Iraq, Lithuania, Slovenia, and Turkey.

Thirty autosomal insertion-deletion (InDel) polymorphisms were analyzed in four populations from Iraq, Lithuania, Slovenia, and Turkey using the commercial kit Investigator® DIPplex. Genotyping issues were encountered for five of the 30 InDels. They were most probably caused by polymorphisms located in the primer binding sites. Population and forensic parameters were calculated. No significant deviations from Hardy-Weinberg equilibrium or significant linkage disequilibrium were detected. The observed heterozygosities ranged from 33% to 61% depending on the marker and the population. The combined probability of exclusion for the 30 markers was 99.7% in all four populations and the matching probabilities were 1 in 3-4×1012 individuals. The multidimensional scaling plot drawn from FST distances showed a good concordance between the relative position of the 15 populations included in the plot and their geographic locations.

Inter-laboratory evaluation of the EUROFORGEN Global ancestry-informative SNP panel by massively parallel sequencing using the Ion PGM™.

The EUROFORGEN Global ancestry-informative SNP (AIM-SNPs) panel is a forensic multiplex of 128 markers designed to differentiate an individual's ancestry from amongst the five continental population groups of Africa, Europe, East Asia, Native America, and Oceania. A custom multiplex of AmpliSeq™ PCR primers was designed for the Global AIM-SNPs to perform massively parallel sequencing using the Ion PGM™ system. This study assessed individual SNP genotyping precision using the Ion PGM™, the forensic sensitivity of the multiplex using dilution series, degraded DNA plus simple mixtures, and the ancestry differentiation power of the final panel design, which required substitution of three original ancestry-informative SNPs with alternatives. Fourteen populations that had not been previously analyzed were genotyped using the custom multiplex and these studies allowed assessment of genotyping performance by comparison of data across five laboratories. Results indicate a low level of genotyping error can still occur from sequence misalignment caused by homopolymeric tracts close to the target SNP, despite careful scrutiny of candidate SNPs at the design stage. Such sequence misalignment required the exclusion of component SNP rs2080161 from the Global AIM-SNPs panel. However, the overall genotyping precision and sensitivity of this custom multiplex indicates the Ion PGM™ assay for the Global AIM-SNPs is highly suitable for forensic ancestry analysis with massively parallel sequencing.

Forensic ancestry analysis with two capillary electrophoresis ancestry informative marker (AIM) panels: Results of a collaborative EDNAP exercise.

There is increasing interest in forensic ancestry tests, which are part of a growing number of DNA analyses that can enhance routine profiling by obtaining additional genetic information about unidentified DNA donors. Nearly all ancestry tests use single nucleotide polymorphisms (SNPs), but these currently rely on SNaPshot single base extension chemistry that can fail to detect mixed DNA. Insertion-deletion polymorphism (Indel) tests have been developed using dye-labeled primers that allow direct capillary electrophoresis detection of PCR products (PCR-to-CE). PCR-to-CE maintains the direct relationship between input DNA and signal strength as each marker is detected with a single dye, so mixed DNA is more reliably detected. We report the results of a collaborative inter-laboratory exercise of 19 participants (15 from the EDNAP European DNA Profiling group) that assessed a 34-plex SNP test using SNaPshot and a 46-plex Indel test using PCR-to-CE. Laboratories were asked to type five samples with different ancestries and detect an additional mixed DNA sample. Statistical inference of ancestry was made by participants using the Snipper online Bayes analysis portal plus an optional PCA module that analyzes the genotype data alongside calculation of Bayes likelihood ratios. Exercise results indicated consistent genotyping performance from both tests, reaching a particularly high level of reliability for the Indel test. SNP genotyping gave 93.5% concordance (compared to the organizing laboratory's data) that rose to 97.3% excluding one laboratory with a large number of miscalled genotypes. Indel genotyping gave a higher concordance rate of 99.8% and a reduced no-call rate compared to SNP analysis. All participants detected the mixture from their Indel peak height data and successfully assigned the correct ancestry to the other samples using Snipper, with the exception of one laboratory with SNP miscalls that incorrectly assigned ancestry of two samples and did not obtain informative likelihood ratios for a third. Therefore, successful ancestry assignments were achieved by participants in 92 of 95 Snipper analyses. This exercise demonstrates that ancestry inference tests based on binary marker sets can be readily adopted by laboratories that already have well-established CE regimes in place. The Indel test proved to be easy to use and allowed all exercise participants to detect the DNA mixture as well as achieving complete and concordant profiles in nearly all cases. Lastly, two participants successfully ran parallel next-generation sequencing analyses (each using different systems) and achieved high levels of genotyping concordance using the exercise PCR primer mixes unmodified.

Population and forensic data for three sets of forensic genetic markers in four ethnic groups from Iran: Persians, Lurs, Kurds and Azeris.

A total of 255 individuals (Persians, Lurs, Kurds and Azeris) from Iran were typed for three sets of forensic genetic markers with the NGM SElect™, DIPplex(®) and Argus X-12 kits. Statistically significant deviations (P≤0.002) from Hardy-Weinberg expectations were observed for the insertion-deletion markers HLD97 and HLD93 after Holm-Sidák correction. Statistically significant (P<0.05) levels of linkage disequilibrium were observed between markers within two of the four studied X-chromosomal linkage groups. AMOVA analyses of the three sets of markers did not show population structure when the individuals were grouped according to their ethnic group. The Iranian population grouped closely to populations living geographically near to Iran based on pairwise FST distances. The matching probabilities ranged from 1 in 3.2×10(7) males by using haplotype frequencies of four X-chromosomal haplogroups to 1 in 3.4×10(21) individuals for the 16 autosomal STRs.

Results for five sets of forensic genetic markers studied in a Greek population sample.

A population sample of 223 Greek individuals was typed for five sets of forensic genetic markers with the kits NGM SElect™, SNPforID 49plex, DIPplex®, Argus X-12 and PowerPlex® Y23. No significant deviation from Hardy-Weinberg expectations was observed for any of the studied markers after Holm-Sidák correction. Statistically significant (P<0.05) levels of linkage disequilibrium were observed between markers within two of the studied X-chromosome linkage groups. AMOVA analyses of the five sets of markers did not show population structure when the individuals were grouped according to their geographic origin. The Greek population grouped closely to the other European populations measured by F(ST)(*) distances. The match probability ranged from a value of 1 in 2×10(7) males by using haplotype frequencies of four X-chromosome haplogroups in males to 1 in 1.73×10(21) individuals for 16 autosomal STRs.

Building a forensic ancestry panel from the ground up: The EUROFORGEN Global AIM-SNP set.

Emerging next-generation sequencing technologies will enable DNA analyses to add pigmentation predictive and ancestry informative (AIM) SNPs to the range of markers detectable from a single PCR test. This prompted us to re-appraise current forensic and genomics AIM-SNPs and from the best sets, to identify the most divergent markers for a five population group differentiation of Africans, Europeans, East Asians, Native Americans and Oceanians by using our own online genome variation browsers. We prioritized careful balancing of population differentiation across the five group comparisons in order to minimize bias when estimating co-ancestry proportions in individuals with admixed ancestries. The differentiation of European from Middle East or South Asian ancestries was not chosen as a characteristic in order to concentrate on introducing Oceanian differentiation for the first time in a forensic AIM set. We describe a complete set of 128 AIM-SNPs that have near identical population-specific divergence across five continentally defined population groups. The full set can be systematically reduced in size, while preserving the most informative markers and the balance of population-specific divergence in at least four groups. We describe subsets of 88, 55, 28, 20 and 12 AIMs, enabling both new and existing SNP genotyping technologies to exploit the best markers identified for forensic ancestry analysis.

Non-uniform phenotyping of D12S391 resolved by second generation sequencing.

Non-uniform phenotyping of five case work samples were observed in the D12S391 locus. The samples were typed at least twice with the AmpFℓSTR NGM SElect PCR Amplification Kit and different alleles were called with GeneMapper ID-X in the different experiments. Detailed analyses of the electropherograms suggested that the individuals were heterozygous with two alleles that differed in size by one nucleotide. This was confirmed by amplifying the samples with the PowerPlex ESX 17 system. D12S391 is a complex STR with variable numbers of AGAT and AGAC repeats. Second generation sequencing revealed that separation of two alleles differing by one nucleotide in length was poor if the number of AGAT repeats in the short allele was higher than in the long allele. A total of 45 individuals with microvariants or off-ladder alleles in D12S391 were sequenced. Thirty different alleles were detected and sixteen of these were not previously reported.

Evaluation of the iPLEX® Sample ID Plus Panel designed for the Sequenom MassARRAY® system. A SNP typing assay developed for human identification and sample tracking based on the SNPforID panel.

Sequenom launched the first commercial SNP typing kit for human identification, named the iPLEX(®) Sample ID Plus Panel. The kit amplifies 47 of the 52 SNPs in the SNPforID panel, amelogenin and two Y-chromosome SNPs in one multiplex PCR. The SNPs were analyzed by single base extension (SBE) and Matrix Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF MS). In this study, we evaluated the accuracy and sensitivity of the iPLEX(®) Sample ID Plus Panel by comparing the typing results of the iPLEX(®) Sample ID Plus Panel with those obtained with our ISO 17025 accredited SNPforID assay. The average call rate for duplicate typing of any one SNPs in the panel was 90.0% when the mass spectra were analyzed automatically with the MassARRAY(®) TYPER 4.0 genotyping software in real time. Two reproducible inconsistencies were observed (error rate: 0.05%) at two different SNP loci. In addition, four inconsistencies were observed once. The optimal amount of template DNA in the PCR was ≥10ng. There was a relatively high risk of allele and locus drop-outs when ≤1ng template DNA was used. We developed an R script with a stringent set of "forensic analysis parameters" based on the peak height and the signal to noise data exported from the TYPER 4.0 software. With the forensic analysis parameters, all inconsistencies were eliminated in reactions with ≥10ng DNA. However, the average call rate decreased to 69.9%. The iPLEX(®) Sample ID Plus Panel was tested on 10 degraded samples from forensic case-work. Two samples could not be typed, presumably because the samples contained PCR and SBE inhibitors. The average call rate was generally lower for degraded DNA samples and the number of inconsistencies higher than for pristine DNA. However, none of the inconsistencies were reproduced and the highest match probability for the degraded samples typed with the panel was 1.7E-9 using the stringent forensic analysis parameters. Although the relatively low sensitivity of the iPLEX(®) Sample ID Plus Panel makes it inappropriate for typing of trace samples from crime scenes, the panel may be interesting for relationship testing and for identification of e.g. samples in biobanks because of the low reagent costs, the limited hands-on time of the iPLEX(®) assay and the automatic analysis of the mass spectra.

Analysis of 49 autosomal SNPs in three ethnic groups from Iran: Persians, Lurs and Kurds.

A total number of 149 individuals from Iran (Persians, Lurs and Kurds) were analyzed for 49 autosomal SNPs using PCR, SBE and capillary electrophoresis. No deviation from Hardy-Weinberg expectations was observed. One SNP pair (rs1015250-rs251934) showed significant linkage disequilibrium in Kurds. However, this was most likely due to chance. High intrapopulation variability and no significant population structure were observed among the three ethnic groups from Iran. Pairwise FST values obtained from the mean numbers of pairwise differences between SNP profiles were calculated for Persians, Lurs, Kurds and eighteen other worldwide populations. For each of the three Iranian ethnic groups, the lowest FST values calculated between an Iranian and non-Iranian populations were observed between Iranians and populations in Iraq and Turkey. The three Iranian ethnic groups grouped together with other West Asian populations in the MDS plot drawn from the FST values. Statistical parameters of forensic interest calculated for the Iranian ethnic groups showed values of the same order of magnitudes as those obtained for Asians. The mean match probability calculated for the 49 SNPs ranged from 1.7x10(-18) for Kurds to 1.3x10(-19) for Persians. Despite the low level of genetic structure observed among Persians, Lurs and Kurds, a single autosomal SNP database should be used with care when extending its forensic application to other Iranian ethnic groups.