Rixun Fang

Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications.

Nonrecombining Y-chromosomal microsatellites (Y-STRs) are widely used to infer population histories, discover genealogical relationships, and identify males for criminal justice purposes. Although a key requirement for their application is reliable mutability knowledge, empirical data are only available for a small number of Y-STRs thus far. To rectify this, we analyzed a large number of 186 Y-STR markers in nearly 2000 DNA-confirmed father-son pairs, covering an overall number of 352,999 meiotic transfers. Following confirmation by DNA sequence analysis, the retrieved mutation data were modeled via a Bayesian approach, resulting in mutation rates from 3.78 × 10(-4) (95% credible interval [CI], 1.38 × 10(-5) - 2.02 × 10(-3)) to 7.44 × 10(-2) (95% CI, 6.51 × 10(-2) - 9.09 × 10(-2)) per marker per generation. With the 924 mutations at 120 Y-STR markers, a nonsignificant excess of repeat losses versus gains (1.16:1), as well as a strong and significant excess of single-repeat versus multirepeat changes (25.23:1), was observed. Although the total repeat number influenced Y-STR locus mutability most strongly, repeat complexity, the length in base pairs of the repeated motif, and the fathers age also contributed to Y-STR mutability. To exemplify how to practically utilize this knowledge, we analyzed the 13 most mutable Y-STRs in an independent sample set and empirically proved their suitability for distinguishing close and distantly related males. This finding is expected to revolutionize Y-chromosomal applications in forensic biology, from previous male lineage differentiation toward future male individual identification.

Evaluation of the Genplex SNP typing system and a 49plex forensic marker panel.

Using a 52 SNP marker set previously developed for forensic analysis, a novel 49plex assay has been developed based on the Genplex typing system, a modification of SNPlex chemistry (both Applied Biosystems) using oligo-ligation of pre-amplified DNA and dye-labeled, mobility modified detection probes. This gives highly predictable electrophoretic mobility of the allelic products generated from the assay to allow detection with standard capillary electrophoresis analyzers. The loci chosen comprise the 48 most informative autosomal SNPs from the SNPforID core discrimination set supplemented with the amelogenin gender marker. These SNPs are evenly distributed across all 22 autosomes, exhibit balanced polymorphisms in three major population groups and have been previously shown to be effective markers for forensic analysis. We tested the accuracy and reproducibility of the Genplex system in three SNPforID laboratories, each using a different Applied Biosystems Genetic Analyzer. Genotyping concordance was measured using replicates of 44 standardized DNA controls and by comparing genotypes for the same samples generated by the TaqMan, SNaPshot and Sequenom iPLEX SNP typing systems. The degree of informativeness of the 48 SNPs for forensic analysis was measured using previously estimated allele frequencies to derive the cumulative match probability and in paternity analysis using 24 trios previously typed with 18 STRs together with three CEPH families with extensive sibships typed with the 15 STRs in the Identifiler kit.

Progress toward an efficient panel of SNPs for ancestry inference.

Many panels of ancestry informative single nucleotide polymorphisms have been proposed in recent years for various purposes including detecting stratification in biomedical studies and determining an individuals ancestry in a forensic context. All of the panels have limitations in their generality and efficiency for routine forensic work. Some panels have used only a few populations to validate them. Some panels are based on very large numbers of SNPs thereby limiting the ability of others to test different populations. We have been working toward an efficient and globally useful panel of ancestry informative markers that is comprised of a small number of highly informative SNPs. We have developed a panel of 55 SNPs analyzed on 73 populations from around the world. We present the details of the panel and discuss its strengths and limitations.

A Whole-Genome Single Nucleotide Polymorphism-Based Approach To Trace and Identify Outbreaks Linked to a Common Salmonella enterica subsp. enterica Serovar Montevideo Pulsed-Field Gel Electrophoresis Type

ABSTRACT In this study, we report a whole-genome single nucleotide polymorphism (SNP)-based evolutionary approach to study the epidemiology of a multistate outbreak of Salmonella enterica subsp. enterica serovar Montevideo. This outbreak included 272 cases that occurred in 44 states between July 2009 and April 2010. A case-control study linked the consumption of salami made with contaminated black and red pepper to the outbreak. We sequenced, on the SOLiD System, 47 isolates with XbaI PFGE pattern JIXX01.0011, a common pulsed-field gel electrophoresis (PFGE) pattern associated with isolates from the outbreak. These isolates represented 20 isolates collected from human sources during the period of the outbreak and 27 control isolates collected from human, food, animal, and environmental sources before the outbreak. Based on 253 high-confidence SNPs, we were able to reconstruct a tip-dated molecular clock phylogeny of the isolates and to assign four human isolates to the actual outbreak. We developed an SNP typing assay to rapidly discriminate between outbreak-related cases and non-outbreak-related cases and tested this assay on an extended panel of 112 isolates. These results suggest that only a very small percentage of the human isolates with the outbreak PFGE pattern and obtained during the outbreak period could be attributed to the actual pepper-related outbreak (20%), while the majority (80%) of the putative cases represented background cases. This study demonstrates that next-generation-based SNP typing provides the resolution and accuracy needed for outbreak investigations of food-borne pathogens that cannot be distinguished by currently used subtyping methods.

Developmental Validation of the Quantifiler® Duo DNA Quantification Kit for Simultaneous Quantification of Total Human and Human Male DNA and Detection of PCR Inhibitors in Biological Samples*

Abstract:  The Quantifiler® Duo DNA Quantification kit enables simultaneous quantification of human DNA and human male DNA as well as detection of inhibitors of PCR in a single real‐time PCR well. Pooled human male genomic DNA is used to generate standard curves for both human (ribonuclease P RNA component H1) and human male (sex determining region Y) specific targets. A shift in the cycle threshold (CT) values for the internal positive control monitors the presence of PCR inhibitors in a sample. The assay is human specific and exhibits a high dynamic range from 0.023 to 50 ng/μL. In addition, the multiplex assay can detect as little as 25 pg/μL of human male DNA in the presence of a 1000‐fold excess of human female DNA. The multiplex assay provides assessment of the DNA extract and guidance for the selection of the appropriate AmpFℓSTR® Amplification Kit to obtain interpretable short tandem repeat profiles.

Escherichia coli serotype O55:H7 diversity supports parallel acquisition of bacteriophage at Shiga toxin phage insertion sites during evolution of the O157:H7 lineage

Enteropathogenic Escherichia coli (EPEC) continues to be a leading cause of mortality and morbidity in children around the world. Two EPEC genomes have been fully sequenced: those of EPEC O127:H6 strain E2348/69 (United Kingdom, 1969) and EPEC O55:H7 strain CB9615 (Germany, 2003). The O55:H7 serotype is a recent precursor to the virulent enterohemorrhagic E. coli O157:H7. To explore the diversity of O55:H7 and better understand the clonal evolution of O157:H7, we fully sequenced EPEC O55:H7 strain RM12579 (California, 1974), which was collected 1 year before the first U.S. isolate of O157:H7 was identified in California. Phage-related sequences accounted for nearly all differences between the two O55:H7 strains. Additionally, O55:H7 and O157:H7 strains were tested for the presence and insertion sites of Shiga toxin gene (stx)-containing bacteriophages. Analysis of non-phage-associated genes supported core elements of previous O157:H7 stepwise evolutionary models, whereas phage composition and insertion analyses suggested a key refinement. Specifically, the placement and presence of lambda-like bacteriophages (including those containing stx) should not be considered stable evolutionary markers or be required in placing O55:H7 and O157:H7 strains within the stepwise evolutionary models. Additionally, we suggest that a 10.9-kb region (block 172) previously believed unique to O55:H7 strains can be used to identify early O157:H7 strains. Finally, we defined two subsets of O55:H7 strains that share an as-yet-unobserved or extinct common ancestor with O157:H7 strains. Exploration of O55:H7 diversity improved our understanding of the evolution of E. coli O157:H7 and suggested a key revision to accommodate existing and future configurations of stx-containing bacteriophages into current models.

Accurate, rapid and high-throughput detection of strain-specific polymorphisms in Bacillus anthracis and Yersinia pestis by next-generation sequencing

BackgroundIn the event of biocrimes or infectious disease outbreaks, high-resolution genetic characterization for identifying the agent and attributing it to a specific source can be crucial for an effective response. Until recently, in-depth genetic characterization required expensive and time-consuming Sanger sequencing of a few strains, followed by genotyping of a small number of marker loci in a panel of isolates at or by gel-based approaches such as pulsed field gel electrophoresis, which by necessity ignores most of the genome. Next-generation, massively parallel sequencing (MPS) technology (specifically the Applied Biosystems sequencing by oligonucleotide ligation and detection (SOLiD™) system) is a powerful investigative tool for rapid, cost-effective and parallel microbial whole-genome characterization.ResultsTo demonstrate the utility of MPS for whole-genome typing of monomorphic pathogens, four Bacillus anthracis and four Yersinia pestis strains were sequenced in parallel. Reads were aligned to complete reference genomes, and genomic variations were identified. Resequencing of the B. anthracis Ames ancestor strain detected no false-positive single-nucleotide polymorphisms (SNPs), and mapping of reads to the Sterne strain correctly identified 98% of the 133 SNPs that are not clustered or associated with repeats. Three geographically distinct B. anthracis strains from the A branch lineage were found to have between 352 and 471 SNPs each, relative to the Ames genome, and one strain harbored a genomic amplification. Sequencing of four Y. pestis strains from the Orientalis lineage identified between 20 and 54 SNPs per strain relative to the CO92 genome, with the single Bolivian isolate having approximately twice as many SNPs as the three more closely related North American strains. Coverage plotting also revealed a common deletion in two strains and an amplification in the Bolivian strain that appear to be due to insertion element-mediated recombination events. Most private SNPs (that is, a, variant found in only one strain in this set) selected for validation by Sanger sequencing were confirmed, although rare false-positive SNPs were associated with variable nucleotide tandem repeats.ConclusionsThe high-throughput, multiplexing capability, and accuracy of this system make it suitable for rapid whole-genome typing of microbial pathogens during a forensic or epidemiological investigation. By interrogating nearly every base of the genome, rare polymorphisms can be reliably discovered, thus facilitating high-resolution strain tracking and strengthening forensic attribution.

Expanding data and resources for forensic use of SNPs in individual identification

The potential value of SNPs for individual identification has been recognized by many researchers and different panels have been proposed. Here we present a new interface in the ALFRED database to access compendia of allele frequencies for several published panels of markers for forensic uses. One of those is our panel of individual identification SNPs (IISNPs) based on samples of 44 populations originating from many parts of the world. Here we also present additional data and additional statistical analyses that continue to support the value of our panel of IISNPs as a universal panel. We also describe initial developments of multiplex methods and various robustness analyses for our 45 marker IISNP panel.

Mini-haplotypes as lineage informative SNPs and ancestry inference SNPs

We propose that haplotyped loci with high heterozygosity can be useful in human identification, especially within families, if recombination is very low among the sites. Three or more SNPs extending over small molecular intervals (<10 KB) can be identified in the human genome to define miniature haplotypes with moderate levels of linkage disequilibrium. Properly selected, these mini-haplotypes (or minihaps) consist of multiple haplotype lineages (alleles) that have evolved from the ancestral human haplotype but show no evidence of recurring recombination, allowing each distinct haplotype to be equated with an allele, all copies of which are essentially identical by descent. Historic recombinants, representing rare events that have drifted to common frequencies over many generations, can be identified in some cases, they do not equate to frequently recurring recombination. We have identified examples in our data collected on various projects and present eight such mini-haplotypes comprised of informative SNPs. We also discuss the ideal characteristics and advantages of minihaps for human familial identification and ancestry inference, and compare them to other types of forensic markers in use and/or that have been proposed. We expect that it is possible to carry out a systematic search and identify a useful panel of mini-haplotypes, with even better properties than the examples presented here.

Autosomal SNP typing of forensic samples with the GenPlex TM HID System: Results of a collaborative study

The GenPlex™ HID System (Applied Biosystems - AB) offers typing of 48 of the 52 SNPforID SNPs and amelogenin. Previous studies have shown a high reproducibility of the GenPlex™ HID System using 250-500pg DNA of good quality. An international exercise was performed by 14 laboratories (9 in Europe and 5 in the US) in order to test the robustness and reliability of the GenPlex™ HID System on forensic samples. Three samples with partly degraded DNA and 10 samples with low amounts of DNA were analyzed in duplicates using various amounts of DNA. In order to compare the performance of the GenPlex™ HID System with the most commonly used STR kits, 500pg of partly degraded DNA from three samples was typed by the laboratories using one or more STR kits. The median SNP typing success rate was 92.3% with 500pg of partly degraded DNA. Three of the fourteen laboratories counted for more than two thirds of the locus dropouts. The median percentage of discrepant results was 0.2% with 500pg degraded DNA. An increasing percentage of locus dropouts and discrepant results were observed when lower amounts of DNA were used. Different success rates were observed for the various SNPs. The rs763869 SNP was the least successful. With the exception of the MiniFiler™ kit (AB), GenPlex™ HID performed better than five other tested STR kits. When partly degraded DNA was analyzed, GenPlex™ HID showed a very low mean mach probability, while all STR kits except MiniFiler™ had very limited discriminatory power.