RESUMO
The identification of tissue-specific differentially methylated regions has significantly contributed to the field of forensic genetics, particularly in body fluid identification crucial for linking evidence to crimes. Among the various approaches to analyzing DNA methylation, the SNaPshot assay has been popularly studied in numerous researches. However, there is a growing interest in exploring alternative methods such as the use of massively parallel sequencing (MPS), which can process a large number of samples simultaneously. This study compares SNaPshot and MPS multiplex assays using nine cytosine-phosphate-guanine markers for body fluid identification. As a result of analyzing 112 samples, including blood, saliva, vaginal fluid, menstrual blood, and semen, both methods demonstrated high sensitivity and specificity, indicating their reliability in forensic investigations. A total of 92.0% samples were correctly identified by both methods. Although both methods accurately identified all blood, saliva, and semen samples, some vaginal fluid samples showed unexpected methylation signals at nontarget loci in addition to the target loci. In the case of menstrual blood samples, due to their complexity, independent typing criteria were applied, and successful menstrual blood typing was possible, whereas a few samples showed profiles similar to vaginal fluid. The MPS method worked better in vaginal fluid samples, and the SNaPshot method performed better in menstrual blood samples. This study offers valuable insights into body fluid identification based on the characteristics of the SNaPshot and MPS methods, which may help in more efficient forensic applications.
RESUMO
In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance in determining cell identity and their pivotal role in numerous disorders, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 60 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways.
RESUMO
Forensic Biology is contingent upon matching DNA profiles between a crime sample and a reference sample. There are several capillary electrophoresis kits available to generate a short tandem repeat (STR) profile from DNA samples, while newer methods using massively parallel sequencing are slowly being implemented in forensic laboratories worldwide. During evaluation of a newer capillary electrophoresis kit, Applied Biosystems™ VeriFiler™ Plus, a discordance was observed in the Penta D locus. The previous kit, Promega PowerPlex 21® System produced a 13.4,14 genotype, whilst VeriFiler™ Plus produced a 14,14 genotype. An expanded investigation into Penta D microvariant alleles revealed that multiple discordances were observed for DNA profiles containing larger x.4 variants. There was full concordance between PowerPlex® 21 and QIAGEN Investigator® 26plex, however discordances were observed between VeriFiler™ Plus and the other three kits tested, including the massively parallel sequencing kit, Verogen ForenSeq® MainstAY. Notably, four of these discordances resulted in null alleles with the VeriFiler™ Plus kit. A review of the Penta D DNA sequences in MainstAY revealed fully concordant microvariant alleles involved deletions within the repeat region, whilst variability in the discordances observed were dependent on the location of the variation outside the repeat region and the analysis method used. Variations observed within the 5' flanking region produced the same allele designation across all capillary electrophoresis kits. However, deletions within the 3' region either produced a null allele for VeriFiler™ Plus where the deletion is thought to overlap the primer binding site, or microvariant alleles for the PowerPlex® 21 and Investigator 26plex kits, which produced longer Penta D amplicons. The discovery of these variations in the Penta D flanking sequences is informative as it increases the awareness of Penta D discordances between different kit chemistries in nominated reference DNA profile comparisons and DNA database searching and matching alike, and provides support for this phenomenon when providing evidence as to the admissibility of such results in trial proceedings.
RESUMO
This study evaluates the performance of analysing surface DNA samples using massively parallel sequencing (MPS) compared to traditional capillary electrophoresis (CE). A total of 30 samples were collected from various surfaces in an office environment and were analysed with CE and MPS. These were compared against 60 reference samples (office inhabitants). To identify contributors, likelihood ratios (LRs) were calculated for MPS and CE data using the probabilistic genotyping software MPSproto and EuroForMix respectively. Although a higher number of sequences/peaks were observed per DNA profile in MPS compared to CE, LR values were found to be lower for MPS data formats. This might be the result of the increased complexity of MPS data, along with a possible elevation of unknown alleles and/or artefacts. The study highlights avenues for improving MPS data quality and analysis to facilitate more robust interpretation of challenging casework-like samples.
RESUMO
Investigative leads are not generated by traditional forensic DNA testing, if the source of the forensic evidence or a 1st degree relative of unidentified human remains is not in the DNA database. In such cases, forensic genetic genealogy (FGG) can provide valuable leads. However, FGG generated genetic data contain private and sensitive information. Therefore, it is essential to deploy approaches that minimize unnecessary disclosure of these data to mitigate potential risks to individual privacy. We recommend protective practices that need not impact effective reporting of relationship identifications. Examples include performing one-to-one comparisons of DNA profiles of third-party samples and evidence samples offline with an "air gap" to the internet and shielding the specific shared single nucleotide polymorphisms (SNP) states and locations by binning adjacent SNPs in forensic reports. Such approaches reduce risk of unwanted access to or reverse engineering of third-party individuals' genetic data and can give these donors greater confidence to support use of their DNA profiles in FGG investigation.
[Box: see text].
RESUMO
Microhaplotypes (MHs), comprising two or more single-nucleotide polymorphisms in a short fragment, are promising forensic markers owing to their remarkable polymorphic nature. Several studies have demonstrated the utility of MHs through massively parallel sequencing (MPS). Nevertheless, the background noise level associated with MHs in MPS, which imposes a practical detection limit for the system, remains uninvestigated. Currently, unique molecular identifier (UMI) systems are known to effectively mitigate background noise by tracking original DNA molecules and facilitating PCR and MPS error corrections. Hence, this study aimed to design a UMI-based amplicon sequencing system, designated MH-UMIseq, which can amplify 46 MHs simultaneously and generate MPS libraries in four steps: barcoding PCR, nuclease reaction, boosting PCR, and indexing PCR. The performance of the MH-UMIseq system was evaluated using the Illumina NextSeq 550 and MiniSeq systems with 31 sets for 5â¯ng, 1â¯ng, and 200â¯pg of input DNA. The fgbio toolkit was used in conjunction with STRait Razor 3.0 and Visual Microhap to analyze the UMI data on MHs. The corresponding average not suppressed noise proportion of MH-UMIseq were 0.1â¯%, 0.3â¯%, and 0.7â¯% for 5â¯ng, 1â¯ng, and 200â¯pg of DNA, respectively, which substantially suppressed the background noise for more than 1â¯ng of DNA. Interestingly, the proportion of not suppressed noise in MH-UMIseq notably decreased as the amount of input DNA increased. The number of UMI families was proportional to the copy number of the template DNA and closely correlated with the system resolution. Therefore, the resolution of MH-UMIseq system is expected to be higher than that of conventional MPS for the deconvolution of mixtures containing more than 1â¯ng of DNA.
Assuntos
Haplótipos , Sequenciamento de Nucleotídeos em Larga Escala , Reação em Cadeia da Polimerase , Polimorfismo de Nucleotídeo Único , Humanos , Análise de Sequência de DNA , Impressões Digitais de DNA/métodosRESUMO
In mass disaster events, forensic DNA laboratories may be called upon to quickly pivot their operations toward identifying bodies and reuniting remains with family members. Ideally, laboratories have considered this possibility in advance and have a plan in place. Compared with traditional short tandem repeat (STR) typing, single nucleotide polymorphisms (SNPs) may be better suited to these disaster victim identification (DVI) scenarios due to their small genomic target size, resulting in an improved success rate in degraded DNA samples. As the landscape of technology has shifted toward DNA sequencing, many forensic laboratories now have benchtop instruments available for massively parallel sequencing (MPS), facilitating this operational pivot from routine forensic STR casework to DVI SNP typing. Herein, we present the commercially available SNP sequencing assays amenable to DVI, we use data simulations to explore the potential for kinship prediction from SNP panels of varying sizes, and we give an example DVI scenario as context for presenting the matrix of considerations: kinship predictive potential, cost, and throughput of current SNP assay options. This information is intended to assist laboratories in choosing a SNP system for disaster preparedness.
RESUMO
O'Donnell-Luria-Rodan (ODLURO) syndrome is an autosomal dominant genetic disorder caused by mutations in the KMT2E (lysine methyltransferase 2E) gene. The Third Xiangya Hospital of Central South University admitted a 12-year and 9-month-old male patient who presented with growth retardation, intellectual disability, and distinctive facial features. Peripheral blood was collected from the patient, and DNA was extracted for genetic testing. Chromosome karyotyping showed 46XY. Whole-exome sequencing and low-coverage massively parallel copy number variation sequencing (CNV-seq) revealed a 506 kb heterozygous deletion in the 7q22.3 region, which includes 6 genes, including KMT2E. The patient was diagnosed with ODLURO syndrome. Both the patient's parents and younger brother had normal clinical phenotypes and genetic test results, indicating that this deletion was a de novo mutation. The clinical and genetic characteristics of this case can help increase clinicians' awareness of ODLURO syndrome.
Assuntos
Deficiência Intelectual , Humanos , Masculino , Deficiência Intelectual/genética , Criança , Histona-Lisina N-Metiltransferase/genética , Mutação , Transtornos do Crescimento/genética , Anormalidades Múltiplas/genética , Cromossomos Humanos Par 7/genética , Cariotipagem , Fenótipo , Variações do Número de Cópias de DNA , Sequenciamento do Exoma , Heterozigoto , Contratura , Microcefalia , FáciesRESUMO
The present study proposes an atomistic molecular dynamics model system of a magnetite (Fe3O4) {111} surface. The effect of temperature on the adsorption process of ammonia (NH3) at low concentrations in the aqueous phase has been considered. The molecular dynamics simulations were carried out using the Clay force field (Clay FF) with a modification for the iron atoms in the NPT ensemble at a pressure of 90 bar. The considered system was heated in a temperature range from 293 to 473 K, and additional relaxations were performed at temperatures of interest. Within the scope of this study, the basic parameters of the magnetite surface were calculated and the distances between the ammonia molecules and the surface were determined. A general idea of the degree and rate of adsorption at specific temperatures was obtained. The calculation results were compared to the experimental data where possible and to other available simulations of adsorption processes on metal oxides.
RESUMO
The large-scale experimental measures of variant functional assays submitted to MaveDB have the potential to provide key information for resolving variants of uncertain significance, but the reporting of results relative to assayed sequence hinders their downstream utility. The Atlas of Variant Effects Alliance mapped multiplexed assays of variant effect data to human reference sequences, creating a robust set of machine-readable homology mappings. This method processed approximately 2.5 million protein and genomic variants in MaveDB, successfully mapping 98.61% of examined variants and disseminating data to resources such as the UCSC Genome Browser and Ensembl Variant Effect Predictor.
RESUMO
Type I diabetes is an autoimmune disease mediated by T-cell destruction of ß cells in pancreatic islets. Currently, there is no known cure, and treatment consists of daily insulin injections. Genome-wide association studies and twin studies have indicated a strong genetic heritability for type I diabetes and implicated several genes. As most strongly associated variants are noncoding, there is still a lack of identification of functional and, therefore, likely causal variants. Given that many of these genetic variants reside in enhancer elements, we have tested 121 CD4+ T-cell enhancer variants associated with T1D. We found four to be functional through massively parallel reporter assays. Three of the enhancer variants weaken activity, while the fourth strengthens activity. We link these to their cognate genes using 3D genome architecture or eQTL data and validate them using CRISPR editing. Validated target genes include CLEC16A and SOCS1. While these genes have been previously implicated in type 1 diabetes and other autoimmune diseases, we show that enhancers controlling their expression harbor functional variants. These variants, therefore, may act as causal type 1 diabetic variants.
Assuntos
Linfócitos T CD4-Positivos , Diabetes Mellitus Tipo 1 , Elementos Facilitadores Genéticos , Predisposição Genética para Doença , Diabetes Mellitus Tipo 1/genética , Diabetes Mellitus Tipo 1/imunologia , Humanos , Linfócitos T CD4-Positivos/imunologia , Linfócitos T CD4-Positivos/metabolismo , Elementos Facilitadores Genéticos/genética , Proteína 1 Supressora da Sinalização de Citocina/genética , Estudo de Associação Genômica Ampla , Lectinas Tipo C/genética , Variação Genética , Polimorfismo de Nucleotídeo Único , Locos de Características QuantitativasRESUMO
DNA mixtures are a common sample type in forensic genetics, and we typically assume that contributors to the mixture are unrelated when calculating the likelihood ratio (LR). However, scenarios involving mixtures with related contributors, such as in family murder or incest cases, can also be encountered. Compared to the mixtures with unrelated contributors, the kinship within the mixture would bring additional challenges for the inference of the number of contributors (NOC) and the construction of probabilistic genotyping models. To evaluate the influence of potential kinship on the individual identification of the person of interest (POI), we conducted simulations of two-person (2â¯P) and three-person (3â¯P) DNA mixtures containing unrelated or related contributors (parent-child, full-sibling, and uncle-nephew) at different mixing ratios (for 2â¯P: 1:1, 4:1, 9:1, and 19:1; for 3â¯P: 1:1:1, 2:1:1, 5:4:1, and 10:5:1), and performed massively parallel sequencing (MPS) using MGIEasy Signature Identification Library Prep Kit on MGI platform. In addition, in silico simulations of mixtures with unrelated and related contributors were also performed. In this study, we evaluated 1): the MPS performance; 2) the influence of multiple genetic markers on determining the presence of related contributors and inferring the NOC within the mixture; 3) the probability distribution of MAC (maximum allele count) and TAC (total allele count) based on in silico mixture profiles; 4) trends in LR values with and without considering kinship in mixtures with related and unrelated contributors; 5) trends in LR values with length- and sequence-based STR genotypes. Results indicated that multiple numbers and types of genetic markers positively influenced kinship and NOC inference in a mixture. The LR values of POI were strongly dependent on the mixing ratio. Non- and correct-kinship hypotheses essentially did not affect the individual identification of the major POI; the correct kinship hypothesis yielded more conservative LR values; the incorrect kinship hypothesis did not necessarily lead to the failure of POI individual identification. However, it is noteworthy that these considerations could lead to uncertain outcomes in the identification of minor contributors. Compared to length-based STR genotyping, using sequence-based STR genotype increases the individual identification power of the POI, concurrently improving the accuracy of mixing ratio inference using EuroForMix. In conclusion, the MGIEasy Signature Identification Library Prep kit demonstrated robust individual identification power, which is a viable MPS panel for forensic DNA mixture interpretations, whether involving unrelated or related contributors.
Assuntos
Impressões Digitais de DNA , DNA , Sequenciamento de Nucleotídeos em Larga Escala , Humanos , DNA/genética , Funções Verossimilhança , Análise de Sequência de DNA , Repetições de Microssatélites , Genótipo , Genética Forense/métodosRESUMO
Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility.
Assuntos
Predisposição Genética para Doença , Estudo de Associação Genômica Ampla , Neoplasias Pulmonares , Polimorfismo de Nucleotídeo Único , Humanos , Neoplasias Pulmonares/genética , Neoplasias Pulmonares/patologia , Desequilíbrio de Ligação , Herança Multifatorial/genética , Linhagem Celular Tumoral , Alelos , Células A549RESUMO
In recent years, probabilistic genotyping software has been adapted for the analysis of massively parallel sequencing (MPS) forensic data. Likelihood ratios (LR) are based on allele frequencies selected from populations of interest. This study provides an outline of sequence-based (SB) allele frequencies for autosomal short tandem repeats (aSTRs) and identity single nucleotide polymorphisms (iSNPs) in 371 individuals from Southern Norway. 27 aSTRs and 94 iSNPs were previously analysed with the ForenSeq™ DNA Signature Prep Kit (Verogen). The number of alleles with frequencies less than 0.05 for sequenced-based alleles was 4.6 times higher than for length-based alleles. Consistent with previous studies, it was observed that sequence-based data (both with and without flanks) exhibited higher allele diversity compared to length-based (LB) data; random match probabilities were lower for SB alleles confirming their advantage to discriminate between individuals. Two alleles in markers D22S1045 and Penta D were observed with SNPs in the 3´ flanking region, which have not been reported before. Also, a novel SNP with a minor allele frequency (MAF) of 0.001, was found in marker TH01. The impact of the sample size on minor allele frequency (MAF) values was studied in 88 iSNPs from Southern Norway (n = 371). The findings were then compared to a larger Norwegian population dataset (n = 15,769). The results showed that the smaller Southern Norway dataset provided similar results, and it was a representative sample. Population structure was analyzed for regions within Southern Norway; FST estimates for aSTR and iSNPs did not indicate any genetic structure. Finally, we investigated the genetic differences between Southern Norway and two other populations: Northern Norway and Denmark. Allele frequencies between these populations were compared, and we found no significant frequency differences (p-values > 0.0001). We also calculated the pairwise FST values per marker and comparisons between Southern and Northern Norway showed small differences. In contrast, the comparisons between Southern Norway and Denmark showed higher FST values for some markers, possibly driven by distinct alleles that were present in only one of the populations. In summary, we propose that allele frequencies from each population considered in this study could be used interchangeably to calculate genotype probabilities.
Assuntos
Impressões Digitais de DNA , Frequência do Gene , Genética Populacional , Sequenciamento de Nucleotídeos em Larga Escala , Repetições de Microssatélites , Polimorfismo de Nucleotídeo Único , Humanos , Noruega , Análise de Sequência de DNA , Funções Verossimilhança , GenótipoRESUMO
Forensic Investigative Genetic Genealogy, a recent sub discipline of forensic genomics, leverages the high throughput and sensitivity of detection of next generation sequencing and established genetic and genealogical approaches to support the identification of human remains from missing persons investigations and investigative lead generation in violent crimes. To facilitate forensic DNA evidence analysis, the ForenSeq® Kintelligence multiplex, consisting of 10,230 SNPs, was developed. Design of the ForenSeq Kintelligence Kit, the MiSeq FGx® Sequencing System and the ForenSeq Universal Analysis Software is described. Developmental validation in accordance with SWGDAM guidelines and forensic quality assurance standards, using single source samples, is reported for the end-to-end workflow from library preparation to data interpretation. Performance metrics support the conclusion that more genetic information can be obtained from challenging samples compared to other commercially available forensic targeted DNA assays developed for capillary electrophoresis (CE) or other current next generation sequencing (NGS) kits due to the higher number of markers, the overall shorter amplicon sizes (97.8% <150â¯bp), and kit design. Data indicate that the multiplex is robust and fit for purpose for a wide range of quantity and quality samples. The ForenSeq Kintelligence Kit and the Universal Analysis Software allow transfer of the genetic component of forensic investigative genetic genealogy to the operational forensic laboratory.
Assuntos
Impressões Digitais de DNA , Sequenciamento de Nucleotídeos em Larga Escala , Polimorfismo de Nucleotídeo Único , Análise de Sequência de DNA , Software , HumanosRESUMO
A large-scale search for the genetic variants with a bias in the representation of alleles in transcriptome data (AE SNPs) and the binding sites in microRNA 3'-UTRs was performed and their functional significance was assessed using massively parallel reporter assay (MPRA). Of the 629,559 associated "SNP-gene" pairs (eQTLs) discovered in the human liver tissue according to the GTEx Analysis V8 data, 4394 polymorphic positions in the 3'-UTRs of the genes, which represent the eQTLs for these genes were selected. The TargetScanHuman 7.0 algorithm and PolymiRTS database were searched for the potential microRNA-binding sites. Of the predicted microRNA sites affected by eQTL-SNPs, we selected 51 sites with the best evidence of functionality according to Ago2-CLIP-seq, CLEAR-CLIP, and eCLIP-seq for RNA-binding proteins. For MPRA, a library of the plasmids carrying the main and alternative alleles for each AE SNP (in total, 102 constructs) was created. Allele-specific expression for 6 SNPs was detected by transfection of the HepG2 cell line with the constructed plasmid library and sequencing of target DNA and RNA sequences using the Illumina (MiSeq) platform.
Assuntos
Regiões 3' não Traduzidas , Alelos , MicroRNAs , Polimorfismo de Nucleotídeo Único , Humanos , Polimorfismo de Nucleotídeo Único/genética , MicroRNAs/genética , MicroRNAs/metabolismo , Células Hep G2 , Sítios de Ligação/genética , Regiões 3' não Traduzidas/genética , Sequenciamento de Nucleotídeos em Larga Escala/métodos , Genes Reporter/genética , Fígado/metabolismo , Proteínas Argonautas/genética , Proteínas Argonautas/metabolismo , Transcriptoma/genéticaRESUMO
Numerous hereditary ophthalmic diseases display significant genetic diversity. Consequently, the utilization of gene panel sequencing allows a greater number of patients to receive a genetic diagnosis for their clinical manifestations. We investigated how to improve the yield of genetic diagnosis through additional gene panel sequencing in hereditary ophthalmic diseases. A gene panel sequencing consisting of a customized hereditary retinopathy panel or hereditary retinitis pigmentosa (RP) panel was prescribed and referred to a CAP-accredited clinical laboratory. If no significant mutations associated with hereditary retinopathy and RP were detected in either panel, additional gene panel sequencing was requested for research use, utilizing the remaining panel. After additional gene panel sequencing, a total of 16 heterozygous or homozygous variants were identified in 15 different genes associated with hereditary ophthalmic diseases. Of 15 patients carrying any candidate variants, the clinical symptoms could be tentatively accounted for by genetic mutations in seven patients. However, in the remaining eight patients, given the in silico mutation predictive analysis, variant allele frequency in gnomAD, inheritance pattern, and genotype-phenotype correlation, fully elucidating the clinical manifestations with the identified rare variant was challenging. Our study highlights the utility of gene panel sequencing in achieving accurate diagnoses for hereditary ophthalmic diseases and enhancing the diagnostic yield through additional gene panel sequencing. Thus, gene panel sequencing can serve as a primary tool for the genetic diagnosis of hereditary ophthalmic diseases, even in cases where a single genetic cause is suspected. With a deeper comprehension of the genetic mechanisms underlying these diseases, it becomes feasible.
RESUMO
PURPOSE OF REVIEW: The purpose of this review is to outline the principles of clinical genetic testing and to provide practical guidance to clinicians in navigating genetic testing for patients with suspected monogenic forms of osteoporosis. RECENT FINDINGS: Heritability assessments and genome-wide association studies have clearly shown the significant contributions of genetic variations to the pathogenesis of osteoporosis. Currently, over 50 monogenic disorders that present primarily with low bone mass and increased risk of fractures have been described. The widespread availability of clinical genetic testing offers a valuable opportunity to correctly diagnose individuals with monogenic forms of osteoporosis, thus instituting appropriate surveillance and treatment. Clinical genetic testing may identify the appropriate diagnosis in a subset of patients with low bone mass, multiple or unusual fractures, and severe or early-onset osteoporosis, and thus clinicians should be aware of how to incorporate such testing into their clinical practices.
Assuntos
Densidade Óssea , Testes Genéticos , Osteoporose , Humanos , Osteoporose/genética , Densidade Óssea/genética , Estudo de Associação Genômica Ampla , Fraturas por Osteoporose/genética , Fraturas Ósseas/genética , Predisposição Genética para DoençaRESUMO
Massively parallel sequencing (MPS) is increasingly applied in forensic short tandem repeat (STR) analysis. The presence of stutter artefacts and other PCR or sequencing errors in the MPS-STR data partly limits the detection of low DNA amounts, e.g., in complex mixtures. Unique molecular identifiers (UMIs) have been applied in several scientific fields to reduce noise in sequencing. UMIs consist of a stretch of random nucleotides, a unique barcode for each starting DNA molecule, that is incorporated in the DNA template using either ligation or PCR. The barcode is used to generate consensus reads, thus removing errors. The SiMSen-Seq (Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing) method relies on PCR-based introduction of UMIs and includes a sophisticated hairpin design to reduce unspecific primer binding as well as PCR protocol adjustments to further optimize the reaction. In this study, SiMSen-Seq is applied to develop a proof-of-concept seven STR multiplex for MPS library preparation and an associated bioinformatics pipeline. Additionally, machine learning (ML) models were evaluated to further improve UMI allele calling. Overall, the seven STR multiplex resulted in complete detection and concordant alleles for 47 single-source samples at 1â¯ng input DNA as well as for low-template samples at 62.5â¯pg input DNA. For twelve challenging mixtures with minor contributions of 10â¯pg to 150â¯pg and ratios of 1-15% relative to the major donor, 99.2% of the expected alleles were detected by applying the UMIs in combination with an ML filter. The main impact of UMIs was a substantially lowered number of artefacts as well as reduced stutter ratios, which were generally below 5% of the parental allele. In conclusion, UMI-based STR sequencing opens new means for improved analysis of challenging crime scene samples including complex mixtures.