Polymerase slippage in relation to the uniformity of tetrameric repeat stretches

Effect existence of aging on stutter ratio evaluated via Bayesian inference

The stochastic behavior of the stutter ratio (SR) in capillary electrophoresis-based DNA typing is currently described and predicted using statistical models in forensic genetics. Clarifying this behavior can help obtain more objective and robust evidence to the court in terms of mixture interpretation. This study aimed to investigate the effect existence of aging on SR via a Bayesian framework. Nail scrapings and clippings were collected from 68 healthy individuals with informed consent. Samples were classified by age-class: young group (0-16 years; n = 36) and older-adult group (>61 years; n = 32). Then, they were compared in terms of their SRs for each simple repeat locus included in GlobalFiler Kit. Bayesian modeling was performed with lognormal distribution model, which implemented multiple linear regression, allele and age-class as explanatory variables. For all simple repeat loci, the median of the posterior distribution of the age-class parameter was a positive value. For CSF1PO and D7S820, the 95% credible interval of the posterior distribution did not include 0. Our data suggested that aging slightly increases the SR. These findings might help elucidate the stochastic behavior of SR.

Modeling the minus two base pair stutter ratio of the D1S1656 locus: A sequence-based mixture distribution model

In this study, we propose a stutter ratio for a minus two base pair stutter (-2bpSR) model of the D1S1656 locus in capillary electrophoresis (CE)-based short tandem repeat (STR) typing. DNA from a total of 108 Japanese individuals was analyzed via massively parallel sequencing to investigate the length of the longest uninterrupted stretch of two base repeat motif (2bpLUS value) within repetitive structures involving the flanking region. Additionally, -2bpSR data was collected using the GlobalFiler Kit on a 3500xL Genetic Analyzer. As a result of sequencing analysis, all alleles were classified into two types by their 2bpLUS values. The -2bpSR differed significantly between the types. Then, we modeled the -2bpSR with a mixture log-normal distribution using the classification of alleles based on the 2bpLUS values. Furthermore, probabilities of the sequence type within each repeat number in the mixture log-normal distribution model were estimated using logistic regression for each of the five major detected populations. This study is expected to enable interpretation of STR typing while considering minus two base pair stutter at the D1S1656 locus.

A reference allelic ladder for Western Capercaillie (Tetrao urogallus) and Black Grouse (Tetrao tetrix) enables linking grouse genetic data across studies

Rapid anthropogenic climate change and progressing habitat degradation are considered top threats to biodiversity. The employment of demanding umbrella species as indicators for ecosystem health is a popular and cost-effective strategy that facilitates continuous monitoring and evaluation within a long-term conservation management scheme. The Western Capercaillie (Tetrao urogallus) and the Black Grouse (Tetrao tetrix) are both considered viable candidates due to their extensive habitat requirements, the possibility for conservative, non-invasive sampling, and their broad popular appeal. Regional population surveys based on genetic data from Short Sequence Repeat (SSR) analysis are being conducted throughout the Palearctic. However, to ensure reliable comparability among laboratories, standardization is required. Here, we report a catalogue of fifty fully characterized reference alleles from twelve SSR loci and the construction of a customizable allelic ladder for genotyping and individualization in Western Capercaillie and Black Grouse. This methodological improvement will help to cost-efficiently generate and collate supraregional data from different grouse surveys and thereby contribute to conservation management. Reference alleles and ladders can be obtained on demand.

Developmental validation of the monSTR identity panel, a forensic STR multiplex assay for massively parallel sequencing

The 21-plex STR panel monSTR was designed for high-fidelity forensic genotyping on the Illumina MiSeq platform. In this study, the panel's performance was validated according to the recommended validation guidelines of the Scientific Working Group for DNA Analysis Methods (SWGDAM). Concordance, repeatability and reproducibility, sensitivity of detection, mixture analysis, species-specificity, and the ability to analyze mock samples were assessed. Sequence data was analyzed using the genotyping software toaSTR. The assay performance was evaluated by measuring the read on-target ratio, the genotype accuracy, the inter-locus balance, the heterozygosity balance, and the signal-to-noise ratio. Results showed that profiles of NIST reference DNA samples as well as GEDNAP proficiency samples were fully concordant with CE-based methods. In addition, inter-run and intra-run variation experiments indicated high precision. Furthermore, full profiles could be obtained using 62.5 pg of DNA input amount with proper inter-locus balance and read on-target ratio; 76.4% of alleles were correctly called with 7.8 pg DNA input amount. It was demonstrated that 94.4% of minor contributor alleles were resolved accurately in a 1:49 mixture. Results suggested that the minor contribution could be precisely calculated based on the minor component allele frequency. Validation results described here demonstrate that the monSTR forensic identity panel is a valid tool for forensic STR genotyping using massively parallel sequencing.

Human Leukocyte Antigen alleles as an aid to STR in complex forensic DNA samples

Human biological samples with multiple contributors remain one of the most challenging aspects of DNA typing within a forensic science context. With the increasing sensitivity of commercially available kits allowing detection of low template DNA, complex mixtures are now a standard component of forensic DNA evidence. Over the years, various methods and techniques have been developed to try to resolve the issue of mixed profiles. However, forensic DNA analysis has relied on the same markers to generate DNA profiles for the past 30 years causing considerable challenges in the deconvolution of complex mixed samples. The future of resolving complicated DNA mixtures may rely on utilising markers that have been previously applied to gene typing of non-forensic relevance. With Massively Parallel Sequencing (MPS), techniques becoming more popular and accessible even epigenetic markers have become a source of interest for forensic scientists. The aim of this review is to consider the potential of alleles from the Human Leukocyte Antigen (HLA) complex as effective forensic markers. While Massively Parallel Sequencing of HLA is routinely used in clinical laboratories in fields such as transplantation, pharmacology or population studies, there have not been any studies testing its suitability for forensic casework samples.

Distribution of Cryptosporidium parvum gp60 subtypes in calf herds of Saxony, Germany

Cryptosporidiosis is a common protozoan parasitic infection that causes diarrhoea in neonatal calves. The high shedding of environmentally resistant oocysts facilitates outbreaks of cryptosporidiosis in humans. In total, 58 farms (512 calves) in Germany (Saxony and Brandenburg) were visited three times each. Faecal samples of pre-weaned calves were microscopically examined for oocysts of Cryptosporidium spp. using Heine staining and were scored with regard to their consistency. Overall, 88.9% of calves tested microscopically positive for Cryptosporidium spp. in at least one sample, and the excretion of oocysts was significantly (P < 0.01) associated with a higher faecal score (diarrhoea). After DNA extraction from pooled farm isolates, 47 samples were successfully subtyped by sequence analysis of the 60 kDa glycoprotein gene (gp60). All isolates belonged to subtype family IIa. IIaA15G2R1 was the most common subtype (present on 66% of the farms), followed by IIaA16G3R1 (13%). Subtypes IIaA14G1R1, IIaA14G2R1, IIaA1612R1, IIaA16G2R1, IIaA17G1R1, IIaA17G2R1, IIaA17G4R1 and IIaA19G2R1 were found sporadically. This is the first description of gp60 subtype IIaA17G4R1 in cattle in Germany.

toaSTR: A web application for forensic STR genotyping by massively parallel sequencing

Massively parallel sequencing (MPS) is emerging within the forensic community as a promising technique for high-resolution short tandem repeat (STR) genotyping, discovering both length and sequence polymorphisms. Conversely, the application of MPS to routine casework poses new challenges to the DNA analyst in view of the complex sequence data that is generated with this technology. We developed the web application toaSTR to help forensic experts work with MPS data simply and efficiently. An intuitive graphical user interface guides through the STR genotyping workflow. This versatile software handles data from various popular MPS platforms and supports both commercial and in-house multiplex PCR kits. Users can define locus-specific stutter thresholds and create custom sets of STR markers to be analyzed. toaSTR's innovative sequence-based stutter model predicts and identifies common stutter variants. The algorithm automatically differentiates biological (iso-)alleles from stutter and other artefacts to assist the interpretation of mixed samples. toaSTR features a comprehensive data visualization with interactive diagrams and a dynamic tabular overview of sequence observations. The software provides an interface to biostatistics tools and enables PDF result export in compliance with the sequence notation recommended by the International Society for Forensic Genetics (ISFG). An initial compatibility and concordance study confirmed the software's independent functionality and precise allele calling with data of different MPS platforms, STR amplification kits, and library preparation chemistries. Discussion of genotyping results for single source and mixed samples demonstrates toaSTR's advantages and includes suggestions for future MPS software development. The beta version of toaSTR is freely accessible at www.toastr.online.

Flanking Variation Influences Rates of Stutter in Simple Repeats

It has been posited that the longest uninterrupted stretch (LUS) of tandem repeats, as defined by the number of exactly matching repeating motif units, is a better predictor of rates of stutter than the parental allele length (PAL). While there are cases where this hypothesis is likely correct, such as the 9.3 allele in the TH01 locus, there can be situations where it may not apply as well. For example, the PAL may capture flanking indel variations while remaining insensitive to polymorphisms in the repeat, and these haplotypic changes may impact the stutter rate. To address this, rates of stutter were contrasted against the LUS as well as the PAL on different flanking haplotypic backgrounds. This study shows that rates of stutter can vary substantially depending on the flanking haplotype, and while there are cases where the LUS is a better predictor of stutter than the PAL, examples to the contrary are apparent in commonly assayed forensic markers. Further, flanking variation that is 7 bp from the repeat region can impact rates of stutter. These findings suggest that non-proximal effects, such as DNA secondary structure, may be impacting the rates of stutter in common forensic short tandem repeat markers.

A forensic multiplex of nine novel pentameric-repeat STRs

Pentameric-repeat short tandem repeats (STRs), consisting of loci with repeat units of five base-pairs, have the advantage of reduced stutter products compared to their tetrameric-repeat STR counterparts. This characteristic potentially helps the interpretation of mixed DNA profiles when minor component alleles may coincide with stutter peaks from the major components. To develop a simple but informative forensic multiplex with the capability to aid mixture interpretation, we designed an 11-plex assay of nine pentameric STRs new to forensic analysis plus two male- specific markers: DYS391 and the Y-Indel rs2032678 used in GlobalFiler™ (Life Technologies). East Asian-specific variation in the recently adopted Y-Indel rs2032678 is reported in this study for the first time in its forensic use as a sex marker. We estimated the levels of variation observed in the nine pentameric STRs in three of the major population groups sampled in the HGDP-CEPH human genome diversity panel: African, European and East Asian (combining individual populations as their sample sizes were too small for STR allele frequency estimations); and we include genotype data from a population sample of Northwest Spain. From this data, forensic informativeness metrics were estimated when applying the nine novel STRs in identification or kinship analyses. The assay was assessed for forensic sensitivity and ability to successfully genotype highly degraded DNA. In the profiles from the 11-plex assay we observed an average 2.15% stutter ratio in all the pentameric loci compared to 7.32% across equivalently-sized tetrameric STRs in the Promega Powerplex^® ESX-17 kit.

Genotyping and interpretation of STR-DNA: Low-template, mixtures and database matches-Twenty years of research and development

The introduction of Short Tandem Repeat (STR) DNA was a revolution within a revolution that transformed forensic DNA profiling into a tool that could be used, for the first time, to create National DNA databases. This transformation would not have been possible without the concurrent development of fluorescent automated sequencers, combined with the ability to multiplex several loci together. Use of the polymerase chain reaction (PCR) increased the sensitivity of the method to enable the analysis of a handful of cells. The first multiplexes were simple: 'the quad', introduced by the defunct UK Forensic Science Service (FSS) in 1994, rapidly followed by a more discriminating 'six-plex' (Second Generation Multiplex) in 1995 that was used to create the world's first national DNA database. The success of the database rapidly outgrew the functionality of the original system - by the year 2000 a new multiplex of ten-loci was introduced to reduce the chance of adventitious matches. The technology was adopted world-wide, albeit with different loci. The political requirement to introduce pan-European databases encouraged standardisation - the development of European Standard Set (ESS) of markers comprising twelve-loci is the latest iteration. Although development has been impressive, the methods used to interpret evidence have lagged behind. For example, the theory to interpret complex DNA profiles (low-level mixtures), had been developed fifteen years ago, but only in the past year or so, are the concepts starting to be widely adopted. A plethora of different models (some commercial and others non-commercial) have appeared. This has led to a confusing 'debate' about the 'best' to use. The different models available are described along with their advantages and disadvantages. A section discusses the development of national DNA databases, along with details of an associated controversy to estimate the strength of evidence of matches. Current methodology is limited to searches of complete profiles - another example where the interpretation of matches has not kept pace with development of theory. STRs have also transformed the area of Disaster Victim Identification (DVI) which frequently requires kinship analysis. However, genotyping efficiency is complicated by complex, degraded DNA profiles. Finally, there is now a detailed understanding of the causes of stochastic effects that cause DNA profiles to exhibit the phenomena of drop-out and drop-in, along with artefacts such as stutters. The phenomena discussed include: heterozygote balance; stutter; degradation; the effect of decreasing quantities of DNA; the dilution effect.

Correlated occurrence and bypass of frame-shifting insertion-deletions (InDels) to give functional proteins

Short insertions and deletions (InDels) comprise an important part of the natural mutational repertoire. InDels are, however, highly deleterious, primarily because two-thirds result in frame-shifts. Bypass through slippage over homonucleotide repeats by transcriptional and/or translational infidelity is known to occur sporadically. However, the overall frequency of bypass and its relation to sequence composition remain unclear. Intriguingly, the occurrence of InDels and the bypass of frame-shifts are mechanistically related - occurring through slippage over repeats by DNA or RNA polymerases, or by the ribosome, respectively. Here, we show that the frequency of frame-shifting InDels, and the frequency by which they are bypassed to give full-length, functional proteins, are indeed highly correlated. Using a laboratory genetic drift, we have exhaustively mapped all InDels that occurred within a single gene. We thus compared the naive InDel repertoire that results from DNA polymerase slippage to the frame-shifting InDels tolerated following selection to maintain protein function. We found that InDels repeatedly occurred, and were bypassed, within homonucleotide repeats of 3–8 bases. The longer the repeat, the higher was the frequency of InDels formation, and the more frequent was their bypass. Besides an expected 8A repeat, other types of repeats, including short ones, and G and C repeats, were bypassed. Although obtained in vitro, our results indicate a direct link between the genetic occurrence of InDels and their phenotypic rescue, thus suggesting a potential role for frame-shifting InDels as bridging evolutionary intermediates.

Homonucleotide repeats are exceptionally prone to insertions and/or deletions of bases (InDels). However, unless they occur in a multiplicity of 3 bases, InDels disrupt the reading frame and are thus expected to be purged from coding regions. Homonucleotide repeats, however, are also vulnerable to slippage by RNA polymerases and the ribosome. Using laboratory evolution techniques, we systematically mapped the occurrence of InDels within a given gene, before and after selection. Our data indicate that frame-shifting InDels were frequently bypassed to give functional proteins at surprisingly high frequencies. Further, we found a strict correlation between the repeat length, the frequency of occurrence of InDels at the DNA level, and the likelihood of bypass by transcriptional/translational slippage. Our results suggest that frame-shifting InDels might comprise functional evolutionary intermediates, and an effective mean of sequence divergence (e.g. when an adjacent InDel restores the frame, resulting in altered sequence and, potentially, in an altered protein structure).

A 27-locus STR assay to meet all United States and European law enforcement agency standards

Different national and international agencies have selected specific STR sets for forensic database use. To enhance database comparison across national and international borders, a 27-locus multiplex system was developed comprising all 15 STR loci of the European standard set, the current 13 STR loci of the CODIS core, the proposed 22 STR loci of the expanded CODIS core, 4 additional commonly used STR loci, and the amelogenin locus. Development required iterative primer design to resolve primer-related artifacts, amplicon sizing, and locus-to-locus balance issues. The 19.5-min assay incorporated newly developed six-dye chemistry analyzed using a novel microfluidic electrophoresis instrument capable of simultaneous detection and discrimination of 8 or more fluorescent dyes. The 27-locus multiplex offers the potential for a new international STR standard permitting laboratories in any jurisdiction to use a single reaction to determine profiles for loci they typically generate plus an expanded common STR profiling set of global interest.

Comparison of southern Chinese Han and Brazilian Caucasian mutation rates at autosomal short tandem repeat loci used in human forensic genetics

The short tandem repeat (STR) loci used in human genetic studies are characterized by having relatively high mutation rates. In particular, to ensure an appropriate evaluation of genetic evidence in parentage and forensic analyses, it is essential to have accurate estimates of the mutation rates associated with the commonly used autosomal and sex chromosome STR loci. Differences in STR mutation rates between different ethnic groups should also be determined. Mutation data from two laboratories working with different ethnic groups were extracted from many meiotic transmissions ascertained for 15 autosomal STR loci currently used in forensic routine. Forty-five thousand and eighty-five trios were checked for the biological consistency of maternity and paternity through the analysis of a minimum of 15 loci. Mutations were scored as paternal, maternal, or ambiguous according to the most parsimonious explanation for the inconsistency, using always the least requiring hypothesis in terms of number of repeat differences. The main findings are: (a) the overall mutation rate across the 15 loci was 9.78 × 10(-4) per gamete per generation (95% CI = 9.30 × 10(-4)-1.03 × 10(-3)), and with just 48 (out of 1,587) exceptions, all of the mutations were single-step; (b) repeat gains were more frequent than losses; (c) longer alleles were found to be more mutable; and (d) the mutation rates differ at some loci between the two ethnic groups. Large worldwide meiotic transmission datasets are still needed to measure allele-specific mutation rates at the STR loci consensually used in forensic genetics.

Characterising stutter in forensic STR multiplexes

Stutter is an artefact seen when amplifying short tandem repeats and typically occurs at one repeat unit shorter in length than the parent allele. In forensic analysis, stutter complicates the analysis of DNA profiles from multiple contributors, known as mixed profiles, a common profile type. Consequently it is important to both understand and predict stutter behaviour in order to improve our understanding of the resolution and interpretation of these profiles. Whilst stutter is well recognised and documented, little information is available that identifies and quantifies what influences the formation of stutter. In this work we use a novel approach to examine this. We have used synthetic oligonucleotides comprising multiple repeat units to test; the influence of repeat number, the influence of repeat sequence and the impact of interruptions to the repeat sequence length. Using multiple replicates allows detailed statistical analysis. We have confirmed a linear relationship between stutter ratio and repeat number. We have shown that increased A-T content increases stutter ratio and that interruptions in repeating sequences decreased stutter ratios to levels similar to the longest uninterrupted repeat stretch. We also found that there was no relationship between stutter ratio and repeat number for a repeat unit with an A-T content of 1/4 and that half of the interrupted repeat sequences stuttered significantly less than their longest uninterrupted repeat stretches. We have applied the knowledge gained to examine specific features of the loci present in the AmpFlSTR(®) SGM Plus(®) multiplex kit used in our laboratory.

Insight into the sialome of the Black Fly, Simulium vittatum

Adaptation to vertebrate blood feeding includes development of a salivary "magic potion" that can disarm host hemostasis and inflammatory reactions. Within the lower Diptera, a vertebrate blood-sucking mode evolved in the Psychodidae (sand flies), Culicidae (mosquitoes), Ceratopogonidae (biting midges), Simuliidae (black flies), and the frog-feeding Corethrellidae. Sialotranscriptome analyses from several species of mosquitoes and sand flies and from one biting midge indicate divergence in the evolution of the blood-sucking salivary potion, manifested in the finding of many unique proteins within each insect family, and even genus. Gene duplication and divergence events are highly prevalent, possibly driven by vertebrate host immune pressure. Within this framework, we describe the sialome (from Greek sialo, saliva) of the black fly Simulium vittatum and discuss the findings within the context of the protein families found in other blood-sucking Diptera. Sequences and results of Blast searches against several protein family databases are given in Supplemental Tables S1 and S2, which can be obtained from http://exon.niaid.nih.gov/transcriptome/S_vittatum/T1/SV-tb1.zip and http://exon.niaid.nih.gov/transcriptome/S_vittatum/T2/SV-tb2.zip .

Variation in the ratio of nucleotide substitution and indel rates across genomes in mammals and bacteria

Rates of nucleotide substitution and insertion/deletion (indel) are known to vary across the functional components of a genome. Little attention has been paid, however, to the quantitative relationship between the two. Here we investigate the ratio of nucleotide substitutions to indels (S/I) in different regions of 4 primates, 70 bacteria, and 8 other genomes. We find that the ratio differs at 5.4-times between coding and noncoding, 3.3-times between conserved and less conserved coding sequences, and 1.46-times between nonrepeat and repeat regions. The S/I ratio is also positively correlated with the level of divergence between the genomes compared. Our results suggest that the S/I ratio may reflect differences in the efficacy of selection against indels. Due to the sensitivity of indel density in different regions, this ratio varies over a much larger range. With the recent discovery suggesting that indels act as local enhancers of mutation in surrounding sequences, nucleotide substitution rates are expected to be accelerated in regions of low constraint, where indels tend to accumulate, but will otherwise be modulated in proportion to the level of a sequence's functional constraint. Indels, therefore, may play a nontrivial role in controlling differences in genetic variation and divergence across functional regions of a genome.

Genotyping of simple sequence repeats--factors implicated in shadow band generation revisited

PCR amplification of microsatellite sequences generates, besides the main product corresponding to allele size, also additional, undesired products usually shorter by multiples of the repeated unit. These extra products known as shadow bands or stutter products may complicate genotyping. The mechanism by which these artifacts are formed is not well understood and so no effective remedy has been found to cope with these spurious products. In this study, using the DNA templates containing the CAG/CTG repeats flanked by gene-specific sequences and universal priming sites, we analyzed the effects of many PCR variables on the shadow band generation. The most important result was that at the decreased temperature of the denaturation step during PCR cycling the shadow bands were either not formed or were strongly suppressed. Several possible sources of this effect are discussed.

Mosaic maternal uniparental isodisomy for chromosome 7q21-qter

Uniparental disomy (UPD) for several human chromosomes is associated with clinical abnormalities. We report the case of a 2-year-old boy with severe intrauterine and post-natal growth retardation (IUGR/PNGR) and highly variable sweat chloride concentrations. The patient was identified as heterozygous for the F508del mutation of the CFTR (cystic fibrosis transmembrane conductance regulator) gene. Unexpectedly, the signal corresponding to the maternally inherited F508del allele appeared much more intense than the paternally derived wild allele. Molecular analysis including polymorphic marker studies, microsatellites and single-nucleotide polymorphisms subsequently showed that the boy was a carrier of a de novo mosaic maternal isodisomy of a chromosome 7 segment while there was a biparental inheritance of the rest of the chromosome. This is the first report of a mosaic partial UPD7. The matUPD7 segment at 7q21-qter extends for 72.7 Mb. The karyotype (550 bands) of our patient was normal, and fluorescence in situ hybridization with probes mapping around the CFTR gene allowed us to rule out a partial duplication. The detection of this chromosomal rearrangement confirms the hypothesis that the 7q31-qter segment is a candidate for the localization of human imprinted genes involved in the control of IUGR and PNGR. It also emphasizes the importance of searching for UPD7 in severe, isolated and unexplained IUGR and PNGR.

Exact tandem repeats analyzer (E-TRA): a new program for DNA sequence mining

Exact Tandem Repeats Analyzer 1.0 (E-TRA) combines sequence motif searches with keywords such as 'organs', 'tissues', 'cell lines' and 'development stages' for finding simple exact tandem repeats as well as non-simple repeats. E-TRA has several advanced repeat search parameters/options compared to other repeat finder programs as it not only accepts GenBank, FASTA and expressed sequence tags (EST) sequence files, but also does analysis of multiple files with multiple sequences. The minimum and maximum tandem repeat motif lengths that E-TRA finds vary from one to one thousand. Advanced user defined parameters/options let the researchers use different minimum motif repeats search criteria for varying motif lengths simultaneously. One of the most interesting features of genomes is the presence of relatively short tandem repeats (TRs). These repeated DNA sequences are found in both prokaryotes and eukaryotes, distributed almost at random throughout the genome. Some of the tandem repeats play important roles in the regulation of gene expression whereas others do not have any known biological function as yet. Nevertheless, they have proven to be very beneficial in DNA profiling and genetic linkage analysis studies. To demonstrate the use of E-TRA, we used 5,465,605 human EST sequences derived from 18,814,550 GenBank EST sequences. Our results indicated that 12.44% (679,800) of the human EST sequences contained simple and non-simple repeat string patterns varying from one to 126 nucleotides in length. The results also revealed that human organs, tissues, cell lines and different developmental stages differed in number of repeats as well as repeat composition, indicating that the distribution of expressed tandem repeats among tissues or organs are not random, thus differing from the un-transcribed repeats found in genomes.