Development of a multiplex panel with 31 multi-allelic InDels for forensic DNA typing.

Insertion/Deletion (InDel) polymorphic genetic markers are abundant in human genomes. Diallelic InDel markers have been widely studied for forensic purposes, yet the low polymorphic information content limits their application and current InDel panels remain to be improved. In this study, multi-allelic InDels located out of low complexity sequence regions were selected in the datasets from East Asian populations, and a multiplex amplification system containing 31 multi-allelic InDel markers and the Amelogenin marker (FA-HID32plex) was constructed and optimized. The preliminary study on sensitivity, species specificity, inhibitor tolerance, mixture resolution, and the detection of degraded samples demonstrates that the FA-HID32plex is highly sensitive, specific, and robust for traces and degraded samples. The combined power of discrimination (CPD) of 31 multi-allelic InDel markers was 0.999 999 999 999 999 999 85, and the cumulative probability of exclusion (CPE) was 0.999 920 in a Chinese Han population, which indicates a high discrimination power. Altogether, the FA-HID32plex panel could provide reliable supplements or stand-alone information in individual identification and paternity testing, especially for challenging samples.

Evaluation of one-step RT-PCR multiplex assay for body fluid identification.

The discrimination of body fluid stains provides crucial evidence during the investigation of criminal cases. Previous studies have demonstrated the practical value of mRNA profiling in body fluid identification. Conventional strategy of mRNA profiling entails reverse transcription and PCR amplification in two separate procedures with different buffer systems. In this study, we subjected the one-step multiplex reverse transcription PCR strategy to mRNA profiling with the inclusion of the same 18 tissue-specific biomarkers in the F18plex system targeting peripheral blood, menstrual blood, vaginal secretion, saliva, semen, and urine. The Qiagen OneStep RT-PCR kit and Titanium One-Step RT-PCR kit were applied to multiplex construction, while reproducible profiling results were obtained with both kits. Compared to the F18plex system, similar expression profiles of biomarkers were obtained in targeted tissues, while expected cross-reaction was observed in non-targeted body fluids. However, CYP2B7P1 and SPINK5 were detected in menstrual blood samples, which was not observed using the F18plex system. Full-profiling results were obtained in all samples using 0.1 ng peripheral blood and semen RNA, and 1 ng menstrual blood, vaginal secretion, saliva, and urine RNA. In conclusion, the application of one-step mRNA profiling strategy could be a reliable and economical method for the simplified, specific, and simultaneous analysis of tissue-specific biomarkers for the discrimination of body fluid origin.

Application of ITO Method in Determination of Special Half-Sibling Relationship.

OBJECTIVES: To identify whether the relationship between Zhang A, Zhang B, Zhang C and Zhang X is the half-sibling relationship whose mother is sister (hereinafter referred to as the special half-sibling relationship) or the common first cousin relationship and discuss the application of ITO method in discriminating the special kinship. METHODS: DNA was extracted from blood stain of four identified individuals, PowerPlex® 21 System and AGCU 21+1 STR kit were used to detect autosomal STR genetic markers. Investigator® Argus X-12 QS kit was used to detect the X chromosome STR genetic markers, the special half-sibling index (SHSI) and first cousin index (FCI) and their likelihood ratio (LR) were calculated by ITO method. RESULTS: The LR results of SHSI to FCI, which were calculated based on autosomal STR genotyping and the analysis of X-STR genotyping results suggested that the relationship between Zhang A, Zhang B, Zhang C and Zhang X was inclined to be a special half-sibling relationship. CONCLUSIONS: For the identification of special kinship, it is necessary to comprehensively apply various genetic markers according to the case. After the conclusion that shared alleles cannot be excluded from the analysis, ITO method can be further used to establish discriminant assumptions according to the specific case to obtain objective and reliable identification opinions.

Genetic analysis of tri-allelic patterns at the CODIS STR loci.

In the routine of autosomal STR genotyping for forensic aims, tri-allelic patterns could be occasionally observed at a single locus in phenotypically normal individuals. Two predominant types of tri-allelic variants have been nominated. Uneven intensities of three alleles are normally considered as the Type 1 pattern, and balanced height of three alleles are considered as the Type 2 pattern. In this study, the prevalence of tri-allelic patterns at the CODIS STR loci was investigated in global populations based on previous reports. The frequencies of the Type 1 and Type 2 pattern manifest a correlation with the germline mutation rates at the CODIS STR loci. The irregular high frequencies of the Type 2 pattern at TPOX with low germline mutation rates could attribute to the stable inheritance of genomic rearrangement from ancestral origin. Furthermore, results from genetic pattern analysis show that only a single allele from STRs with the Type 1 pattern could be transmitted from parents to offsprings, while a single allele and a combination of two alleles from STRs with the Type 2 pattern present an equal opportunity of transmission from parents to offsprings. Altogether, these results provide a genetic portrait of STRs with tri-allelic patterns, which will help the genetic interpretation of tri-allelic patterns in forensic practice.

Genetic analysis of type 2 tri-allelic pattern at TPOX locus in the Chinese Han population.

Tri-allelic patterns can occasionally be observed during the profiling of short tandem repeats (STRs) in routine forensic practice. In previous studies, the Type 2 tri-allelic pattern at TPOX has been widely studied in African and Brazilian populations. In this study, we investigated the incidence, rearrangement, and inheritance of the Type 2 tri-allelic pattern at the TPOX locus in a Chinese Han population. The frequency of the Type 2 pattern at TPOX was approximately 0.0189%, and the major extra allele was allele 11 in the Chinese Han population. Two major allelic combinations, 8/11 and 11/12, were observed, which are different from the configuration of that in both African and Brazilian populations. Tight linkage between alleles 11 and 12 was observed in the majority of the Type 2 pattern at TPOX in the Chinese Han population, while the location of the extra copy on chromosome 2 was validated, which shows an identical ancestral origin. The excess allelic combination 8/11 implies a homogeneous origin and tight linkage relationship. However, the rearrangement in the Type 2 pattern with the 8/11 allelic combination remained unknown. Altogether, these results show the configuration of the Type 2 tri-allelic pattern at the TPOX locus in the Chinese Han population, which will assist in the understanding of the Type 2 tri-allelic pattern at the TPOX locus in the global population.

Characterization of the extra copy of TPOX locus with tri-allelic pattern.

BACKGROUND: An STR locus with tri-allelic pattern is occasionally observed in routine forensic casework. The extra copy of TPOX locus with tri-allelic pattern in populations has been assumed to be inserted into an X chromosome, which took place forth before the Bantu expansion in Africa. Nonetheless, the exact location of the duplication and the form of rearrangement in the human genome has not been clarified yet. RESULTS: In this study, we investigated the extra copy of type 2 tri-allelic pattern at TPOX in various populations. While allele 10 is the major third allele in Africa, allele 11 appears more frequent in America and overwhelming in Chinese and Korean populations, which might attribute to the population substructures. Results from the investigation of family cases showed that the transmission of the extra allele had a similar genetic pattern of autosomal genes. Furthermore, a whole-genome sequencing followed by bioinformatics analysis revealed that the intact form of chromosomal duplication and rearrangement occurred ~ 407 kb away from the authentic TPOX locus on chromosome 2 in two cases. The breakpoints of the insertion were further validated in most other tri-allelic subjects, which can imply the identical origin from the ancestral extra copy. Nevertheless, de novo chromosomal duplication and rearrangement at thyroid peroxidase gene occur in populations. CONCLUSIONS: Instead of the extra allele 10 in African populations, the main third allele at TPOX with tri-allelic pattern is allele 11 in Chinese and Korean populations. The insertion of the extra copy into chromosome 2 occurs in most subjects with tri-allelic pattern at TPOX and demonstrates the transmission of the third allele from parents to offspring. The breakpoints of the ancestral extra copy are defined, which shows evidence of its inheritance from African populations. In addition, the simple validation method would help improve tri-allelic pattern calling, distinguish de novo chromosomal rearrangements, and also count the frequencies among different geographic regions. Therefore, the statistical interpretation of tri-allelic pattern at TPOX could be enhanced during forensic practice.

Characterization of tissue-specific biomarkers with the expression of circRNAs in forensically relevant body fluids.

Messenger RNA (mRNA) markers have been extensively investigated for the identification of forensically relevant body fluids and tissues based on their expression profiles among cell types. As products of the backsplicing of pre-mRNAs, circular RNAs (circRNAs) share exonic sequences with their linear counterparts. The inclusion of circRNAs in mRNA profiling is shown to facilitate the detection of biomarkers in the identification of body fluids. In this study, we identified the expression of circRNAs of 14 out of 45 biomarkers from five body fluid types using outward-facing primer sets and revealed the ratio of circular to total transcripts of biomarkers by RNase R treatment. Furthermore, our results of qPCR analysis show that the inclusion of circRNAs in the detection of biomarkers, including HBA and ALAS2 for blood; MMP7 and MMP10 for menstrual blood; HTN3 for saliva; SPINK5, SERPINB3, ESR1, and CYP2B7P1 for vaginal secretions; TGM4, KLK3, and PRM2 for semen; and SLC22A6 and MIOX for urine, does not impair the specificity of these biomarkers. Additionally, a high copy number of targets from linear transcripts could be employed to increase the detection sensitivity of TGM4 and KLK3 with a low expression level of circRNAs in urine samples. Altogether, these results will help with the development of robust multiplex assays for body fluid identification.

Evaluation of the inclusion of circular RNAs in mRNA profiling in forensic body fluid identification.

The use of messenger RNA (mRNA) profiling is considered a promising method in the identification of forensically relevant body fluids which can provide crucial information for reconstructing a potential crime. However, casework samples are usually of limited quantity or have been subjected to degradation, which requires improvement of body fluid identification. Circular RNAs (circRNAs), a class of products from the backsplicing of pre-mRNAs, are shown to have high abundance, remarkable stability, and cell type-specific expression in human cells. In this study, we investigated whether the inclusion of circRNAs in mRNA profiling improve the detection of biomarkers including δ-aminolevulinate synthase 2 (ALAS2) and matrix metallopeptidase 7 (MMP7) in body fluid identification. The major circRNAs of ALAS2 and MMP7 were first identified and primer sets for the simultaneous detection of linear and circular transcripts were developed. The inclusion of circRNAs in mRNA profiling showed improved detection sensitivity and stability of biomarkers revealed by using serial dilutions, mixed samples, and menstrual bloodstains as well as degraded and aged samples. Therefore, the inclusion of circRNAs in mRNA profiling should facilitate the detection of mRNA markers in forensic body fluid identification.

Evaluating the potential of housekeeping genes, rRNAs, snRNAs, microRNAs and circRNAs as reference genes for the estimation of PMI.

The precise estimation of postmortem interval (PMI) is a critical step in death investigation of forensic cases. Detecting the degradation of RNA in tissues by real time quantitative polymerase chain reaction (RT-qPCR) technology provides a new theoretical basis for estimation of PMI. However, most commonly used reference genes degrade over time, while previous studies seldom consider this when selecting suitable reference genes for the estimation of PMI. Studies have shown microRNAs (miRNAs) are very stable and circular RNAs (circRNAs) have recently emerged as a novel class of RNAs with high stability. We aimed to evaluate the stability of the two kinds of RNAs and normal reference genes using geNorm and NormFinder algorithms to identify tissue-specific reference genes for PMI estimation. The content of candidate RNAs from mouse heart, liver and skeletal muscle tissues were dynamically examined in 8 consecutive days after death. Among the 11 candidate genes (ß-actin, Gapdh, Rps18, 5S, 18S, U6, miR-133a, miR-122, circ-AFF1, LC-Ogdh and LC-LRP6), the following genes showed prioritized stability: miR-122, miR-133a and 18S in heart tissues; LC-Ogdh, circ-AFF1 and miR-122 in liver tissues; and miR-133a, circ-AFF1 and LC-LRP6 in skeletal muscle tissues. Our results suggested that miRNAs and circRNAs were more stable as reference genes than other kinds of RNAs regarding PMI estimation. The appropriate internal control genes were not completely the same across tissue types.

Mutation analysis of 19 autosomal short tandem repeats in Chinese Han population from Shanghai.

The mutation of short tandem repeat (STR) loci is affected by several factors, such as sex, age, and DNA architectures. Previous studies have shown a different profile of mutation rates at autosomal STR loci among populations. It is important to provide population data and reveal underlying factors influencing the evaluation of STR mutation rates. In this study, we performed a comprehensive analysis on the mutation of 19 autosomal STR loci through 124,773 parent-child allelic transfers from 5846 paternity testing cases. A total of 197 mutations were observed including 187 single-step mutations. The observed mutation rates ranged from 0.15 × 10-3 (TH01) to 4.57 × 10-3 (FGA), and the average mutation rate across all the 19 loci was 1.58 × 10-3. Furthermore, the average mutation rate of STR loci increases with the paternal conception ages and remains relatively stable in different maternal age groups, which suggest the profile of paternal conception ages as a potential factor influencing the evaluation of STR mutation rates and the ratio of paternal versus maternal mutation rate in populations. Multidimensional scaling analysis (MDS) shows a difference in the profile of mutation rates at 13 CODIS STR loci among ethnical groups. Based on our data, our results support that short alleles are biased towards expansion mutation and longer alleles favor contraction mutation. In conclusion, our results provide useful information for further investigation on STR mutation in forensic genetics and population genetics.

[Research Progress of circRNA and Its Significance in Forensic Science].

RNA has received more attention in the field of forensic medicine and the development of the new biological markers based on RNA shows great significance in the analysis of complex cases. circular RNA (circRNA) is a kind of non-coding RNA which is widely reported recently. Although the regulatory mechanisms of generation and expression are not fully clear, the existing research indicates that circRNA has important biological functions. CircRNA has a cell-type-specific expression with great stability and a high expression level, which makes it meaningful in forensic applications potentially. In this paper, the research progress, the generation and regulation of circRNA as well as its biological characteristics and functions are summarized, which will provide references for related studies and forensic applications.

[Application of the Peak Area Ratio of STR Loci to Amelogenin Locus in the Estimation of DNA Degradation].

OBJECTIVE: To explore the change rules of peak area ratio of STR loci to Amelogenin (AMEL) locus (STR/AMEL), a sex-determining gene in DNA degradation, and to evaluate the application of STR/AMEL value in the estimation of DNA degradation degree. METHODS: DNA was extracted from iliopsoas, and the variations of STR/AMEL value (Penta E/AMEL, Penta D/AMEL, FGA/AMEL) were analyzed after the artificial degradation was made by DNase I, and the changes of these three ratios of the iliopsoas naturally degraded in an outdoor environment were also analyzed. The regression curves were analyzed using the periods of DNA degradation and outside the body as the independent variable (x) and the STR/AMEL value as the dependent variable (y) and three curve equations under two conditions were established. RESULTS: Both under the conditions of artificial and natural degradation, STR/AMEL value had a negative relationship with the degradation time. The relationship between STR/AMEL and degradation time can be well simulated by the cubic function. R2 was over 0.99 under controlled degradation condition and over 0.86 under natural degradation condition. CONCLUSION: The STR/AMEL value (Penta E/AMEL, Penta D/AMEL, FGA/AMEL) is negatively related with the DNA degradation degree, which follows mathematical regression models strictly, and it might be applied to evaluate the DNA degradation degree.

Identification and characterization of the highly polymorphic locus D14S739 in the Han Chinese population.

AIM: To systemically select and evaluate short tandem repeats (STRs) on the chromosome 14 and obtain new STR loci as expanded genotyping markers for forensic application. METHODS: STRs on the chromosome 14 were filtered from Tandem Repeats Database and further selected based on their positions on the chromosome, repeat patterns of the core sequences, sequence homology of the flanking regions, and suitability of flanking regions in primer design. The STR locus with the highest heterozygosity and polymorphism information content (PIC) was selected for further analysis of genetic polymorphism, forensic parameters, and the core sequence. RESULTS: Among 26 STR loci selected as candidates, D14S739 had the highest heterozygosity (0.8691) and PIC (0.8432), and showed no deviation from the Hardy-Weinberg equilibrium. 14 alleles were observed, ranging in size from 21 to 34 tetranucleotide units in the core region of (GATA)9-18 (GACA)7-12 GACG (GACA)2 GATA. Paternity testing showed no mutations. CONCLUSION: D14S739 is a highly informative STR locus and could be a suitable genetic marker for forensic applications in the Han Chinese population.

Validation of phylogenetic informative Y-InDels in Y-chromosomal haplogroup O-M175.

The Y-chromosomal haplogroup tree, which consists of a group of Y-chromosomal loci with phylogenetic information, has been widely applied in anthropology, archaeology and population genetics. With the continuous updating of the phylogenetic structure, Y-chromosomal haplogroup tree provides more information for recalling the biogeographical origin of Y chromosomes. Generally, Y-chromosomal insertion-deletion polymorphisms (Y-InDels) are genetically stable as Y-chromosomal single nucleotide polymorphisms (Y-SNPs), and therefore carry mutations that can accumulate over generations. In this study, potential phylogenetic informative Y-InDels were filtered out in haplogroup O-M175, which is dominant in East Asia, based on population data retrieved from the 1000 Genomes Project. A group of 22 phylogenetic informative Y-InDels were identified and then assigned to their corresponding subclades of haplogroup O-M175, which provided a supplement for the update and application of Y-chromosomal markers. Especially, four Y-InDels were introduced to define subclades determined using a single Y-SNP.

Investigation of Linear Amplification Using Abasic Site-Containing Primers Coupled to Routine STR Typing for LT-DNA Analysis.

Obtaining a full short tandem repeat (STR) profile from a low template DNA (LT-DNA) still presents a challenge for conventional methods due to significant stochastic effects and polymerase slippage. A novel amplification method with a lower cost and higher accuracy is required to improve the DNA amount. Previous studies suggested that DNA polymerases without bypass activity could not perform processive DNA synthesis beyond abasic sites in vitro and our results showed a lack of bypass activity for Phusion, Pfu and KAPA DNA polymerases in this study. Based on this feature, we developed a novel linear amplification method, termed Linear Aamplification for double-stranded DNA using primers with abasic sites near 3' end (abLAFD), to limit the replication error. The amplification efficiency was evaluated by qPCR analysis with a result of approximately a 130-fold increase in target DNA. In a LT-DNA analysis, the abLAFD method can be employed as a pre-PCR. Similar to nested PCRs, primer sets used for the abLAFD method were designed as external primers suitable for commercial multiplex STR amplification assays. The practical performance of the abLAFD method was evaluated by coupling it to a routine PP21 STR analysis using 50 pg and 25 pg DNA. Compared to reference profiles, all abLAFD profiles showed significantly recovered alleles, increased average peak height and heterozygote balance with a comparable stutter ratio. Altogether, our results support the theory that the abLAFD method is a promising strategy coupled to STR typing for forensic LT-DNA analysis.

Identification and assessment of a subset of Y-SNPs with recurrent mutation for forensic purpose.

Y chromosome has an important role in the forensic practice due to its unique paternal inheritance pattern. Y-chromosomal single nucleotide polymorphisms (Y-SNPs) could provide supplementary information while the application of Y-chromosomal STR (Y-STR) haplotypes encounter their limitations. Y-SNPs with recurrent mutation can be seen in different Y-chromosomal haplogroups, which might help discriminate different paternal pedigrees. In this study, a host of candidate Y-SNPs with recurrent mutation were obtained based on population data from 1000 Genome Project. Further, 8 Y-SNPs from a small part of candidates were confirmed to be polymorphic in 2 or more Y-chromosomal haplogroups (sub-haplogroups) in the Chinese Han population. With a haplotype diversity value of 0.9367, the investigated subset of Y-SNPs with recurrent mutation shows a high discrimination power. Therefore, Y-SNPs with recurrent mutation should function as useful markers to provide information in the forensic applications.

Identification and Characterization of Nine Novel X-Chromosomal Short Tandem Repeats on Xp21.1, Xq21.31, and Xq23 Regions.

The application of X-chromosomal short tandem repeats (X-STRs) has been recognized as a powerful tool in complex kinship testing. To support further development of X-STR analysis in forensic use, we identified nine novel X-STRs, which could be clustered into three linkage groups on Xp21.1, Xq21.31, and Xq23. A multiplex PCR system was built based on the electrophoresis. A total of 198 unrelated Shanghai Han samples along with 168 samples from 43 families was collected to investigate the genetic polymorphism and forensic parameters of the nine loci. Allele numbers ranged from 5 to 12, and amplicon sizes ranged from 146 to 477 bp. The multiplex showed high values for the combined power of discrimination (0.99997977 in males and 0.99999999 in females) and combined mean exclusion chances (0.99997918 and 0.99997821 in trios, 0.99984939 in duos, and 0.99984200 in deficiency cases). The linkage between all pairs of loci was estimated via Kosambi mapping function and linkage disequilibrium test, and further investigated through the family study. The data from 43 families strongly demonstrated an independent transmission between LGs and a tight linkage among loci within the same LG. All these results support that the newly described X-STRs and the multiplex system are highly promising for further forensic use.

Variants in linkage status at D5S818 detected by multiple STR kits comparison and Sanger sequencing.

BACKGROUND: D5S818 discrepancies have been reported in forensic parental testing due to null alleles. However, more cases may be ignored since proportional null alleles were missed without detection of heredity discrepancy between parents and offspring. RESULTS: In this study, null allele 12 at D5S818 was detected by the PowerPlex® 21 System with a higher occurrence rate on the basis of review on 2824 samples from the 1282 routine cases in Chinese Han population. Sequencing results revealed novel variant of guanine (G) into adenine (A) in the 7th [AGAT] repeats in the core repeat region accompanied by rs1187948322 in the samples with null allele 12. CONCLUSIONS: Forensic STR typing may benefit from this discovery: (1) primer design of CE profiling system could be improved for sensitive population and (2) polymorphic information could be enriched for the accuracy and precision of NGS genotyping system. Peak area of D5S818 was also analyzed through different commercial STR kits. It is suggested that more attention should be paid on observed homozygosity with reduced peak area, especially for the samples from Chinese Han population.

Assessing Autosomal InDel Loci With Multiple Insertions or Deletions of Random DNA Sequences in Human Genome.

Multiple mutational events of insertion/deletion occurring at or around InDel sites could form multi-allelic InDels and multi-InDels (abbreviated as MM-InDels), while InDels with random DNA sequences could imply a unique mutation event at these loci. In this study, preliminary investigation of MM-InDels with random sequences was conducted using high-throughput phased data from the 1000 Genomes Project. A total of 3,599 multi-allelic InDels and 6,375 multi-InDels were filtered with multiple alleles. A vast majority of the obtained MM-InDels (85.59%) presented 3 alleles, which implies that only one secondary insertion or deletion mutation event occurred at these loci. The more frequent presence of two adjacent InDel loci was observed within 20 bp. MM-InDels with random sequences presented an uneven distribution across the genome and showed a correlation with InDels, SNPs, recombination rate, and GC content. The average allelic frequencies and prevalence of multi-allelic InDels and multi-InDels presented similar distribution patterns in different populations. Altogether, MM-InDels with random sequences can provide useful information for population resolution.

Calculation of the Paternity Index for STR with tri-allelic patterns in paternity testing.

The calculation of the paternity index (PI) value of common bi-allelic genotypes at STR loci has been standardized in paternity cases. However, for tri-allelic patterns, a rare category of genotyping aberration in forensic practice, the statistical analysis in paternity testing remains disputed. The Type 1 tri-allelic pattern generally results from somatic mutation in the early stage of individual development. The Type 2 tri-allelic pattern is commonly generated by segmental duplication in the genome. In this study, practical and theoretical aspects of the evaluation of evidence concerning the Type 1 and Type 2 tri-allelic patterns in healthy individuals are discussed based on the likelihood ratio (LR) in different categories of kinship cases. The calculation of the PI value concerning tri-allelic genotypes is formulated according to the generation and genetic transmission of tri-allelic patterns. Meanwhile, a package tool named TriPI is developed to assist the calculation of the PI value in paternity testing concerning tri-allelic subjects, which could benefit the evaluation of the weight of evidence in the interpretation of tri-allelic pattern in forensic practice.