Integrative lncRNA, circRNA, and mRNA analysis reveals expression profiles of six forensic body fluids/tissue.

RNAs have attracted much attention in forensic body fluid/tissue identification (BFID) due to their tissue-specific expression characteristics. Among RNAs, long RNAs (e.g., mRNA) have a higher probability of containing more polymorphic sites that can be used to assign the specific donor of the body fluid/tissue. However, few studies have characterized their overall profiles in forensic science. In this study, we sequenced the transcriptomes of 30 samples from venous blood, menstrual blood, semen, saliva, vaginal secretion, and skin tissue, obtaining a comprehensive picture of mRNA, lncRNA, and circRNA profiles. A total of 90,305 mRNAs, 102,906 lncRNAs (including 19,549 novel lncRNAs), and 40,204 circRNAs were detected. RNA type distribution, length distribution, and expression distribution were presented according to their annotation and expression level, and many novel body fluid/tissue-specific RNA markers were identified. Furthermore, the cognate relations among the three RNAs were analyzed according to gene annotations. Finally, SNPs and InDels from RNA transcripts were genotyped, and 21,611 multi-SNP and 4,471 multi-InDel transcriptomic microhaplotypes (tMHs) were identified. These results provide a comprehensive understanding of transcriptome profiles, which could provide new avenues for tracing the origin of the body fluid/tissue and identifying an individual.

Targeted S5 RNA sequencing assay for the identification and direct association of common body fluids with DNA donors in mixtures.

The evidentiary value of DNA profiles varies depending upon the context in which the DNA was found. Linking a DNA profile to a particular cellular phenotype in mixtures may aid in assessing its evidentiary relevance and value. We report the development of two dual-function high-resolution messenger RNA (mRNA) sequencing assays that can each identify the presence of 6 body fluids/tissues (blood, semen, saliva, vaginal secretions, menstrual blood, skin) and, via coding region SNPs (cSNPs) present in the body fluid-specific mRNA transcripts, directly associate particular body fluids with their specific DNA donors in mixtures. The original blood, semen, and saliva (BSS) assay contains 23 cSNPs for blood, semen, and saliva, while the expanded 6F (all 6 fluids/tissues) assay encompasses the BSS assay and also contains 23 additional cSNPs for vaginal secretions, menstrual blood, and skin. Software tools were developed to infer the identity of the body fluids present as well as providing the corresponding cSNP genotypes. Concomitant genomic DNA assays (BSS-d and 6F-d), required to genotype the same cSNPs from persons of interest/inferred contributors to the body fluid mixture, were also developed. Body fluid specificity was demonstrated by the ability to identify the body fluid origin of single-source and two-fluid admixtures. The discriminatory power (European Caucasians) for all body fluids is 0.957-0.997, with linkage disequilibrium considered. Reciprocal body fluid admixtures (mixture pairs with the same two donors but reversed body fluid types) were used to demonstrate the ability to identify the body fluid source of origin as well as associate the donor of each of the two fluids.

Analysis of mutations spectrum in the ATP7B gene in patients with Wilson disease using massively parallel sequencing.

The study aimed to search for mutations in the ATP7B gene using massively parallel sequencing in patients with Wilson disease in the Tomsk region. For 42 patients with suspected Wilson's disease (aged from 1 to 33 years) was performed molecular genetic analysis. Enrichment of the interest genome regions was carried out by the long-range PCR. DNA libraries with ligated adapters were constructed with Nextera DNA Flex (Illumina, USA) kit. Sequencing was performed on the Illumina MiSeq platform (Illumina, USA). As a result of this work, we identified 9 pathogenic genetic variants. All variants were previously described in the literature and were found in patients with Wilson's disease. Five missense mutations, one splice site mutation, and 3 frameshift mutations were identified. In patients with Wilson's disease in the Tomsk region, the most common variant was c.3207C>A, this variant is the most common both in the Russian Federation and in other European populations. Also, a pathogenic variant c.3036dupC was found, which is probably endemic to the Russian Federation.

Re-evaluation of single nucleotide variants and identification of structural variants in a cohort of 45 sudden unexplained death cases.

Sudden unexplained death (SUD) takes up a considerable part in overall sudden death cases, especially in adolescents and young adults. During the past decade, many channelopathy- and cardiomyopathy-associated single nucleotide variants (SNVs) have been identified in SUD studies by means of postmortem molecular autopsy, yet the number of cases that remain inconclusive is still high. Recent studies had suggested that structural variants (SVs) might play an important role in SUD, but there is no consensus on the impact of SVs on inherited cardiac diseases. In this study, we searched for potentially pathogenic SVs in 244 genes associated with cardiac diseases. Whole-exome sequencing and appropriate data analysis were performed in 45 SUD cases. Re-analysis of the exome data according to the current ACMG guidelines identified 14 pathogenic or likely pathogenic variants in 10 (22.2%) out of the 45 SUD cases, whereof 2 (4.4%) individuals had variants with likely functional effects in the channelopathy-associated genes SCN5A and TRDN and 1 (2.2%) individual in the cardiomyopathy-associated gene DTNA. In addition, 18 structural variants (SVs) were identified in 15 out of the 45 individuals. Two SVs with likely functional impairment were found in the coding regions of PDSS2 and TRPM4 in 2 SUD cases (4.4%). Both were identified as heterozygous deletions, which were confirmed by multiplex ligation-dependent probe amplification. In conclusion, our findings support that SVs could contribute to the pathology of the sudden death event in some of the cases and therefore should be investigated on a routine basis in suspected SUD cases.

A novel forensic panel of 186-plex SNPs and 123-plex STR loci based on massively parallel sequencing.

The MiSeq® FGX Forensic system and the HID-Ion AmpliSeq Panel were previously developed for massively parallel sequencing (MPS) for forensic casework. Among the three major sequencing platforms, BGISEQ-500TM, which is based on multiple PCRs, is still lacking in forensics. Here, a novel forensic panel was constructed to detect 186 single-nucleotide polymorphisms (SNPs) and 123 short tandem repeats (STRs) with MPS technology on the BGISEQ-500™ platform. First, the library preparation, sequencing process, and data analysis were performed, focusing on the average depth of coverage and heterozygote balance. We calculated the allelic frequencies and forensic parameters of STR and SNP loci in 73 unrelated Chinese Han individuals. In addition, performance was evaluated with accuracy, uniformity, sensitivity, PCR inhibitor, repeatability and reproducibility, mixtures, degraded samples, case-type samples, and pedigree analyses. The results showed that 100% accurate and concordant genotypes can be obtained, and the loci with an abundance in the interquartile range accounted for 92.90% of the total, suggesting reliable uniformity in this panel. We obtained a locus detection rate that was higher than 98.78% from 78 pg of input DNA, and the optimal amount was 1.25-10 ng. The maximum concentrations of hematin and humic acid were 200 and 100 µM, respectively (the ratios of detected loci were 96.52% and 92.41%), in this panel. As a mixture, compared with those of SNPs, minor-contributor alleles of STRs could be detected at higher levels. For the degraded sample, the ratio of detected loci was 98.41%, and most profiles from case-type samples were not significantly different in abundance in our studies. As a whole, this panel showed high-performance, reliable, robust, repeatable, and reproducible results, which are sufficient for paternity testing, individual identification, and use for potentially degraded samples in forensic science.

A microhap panel for kinship analysis through massively parallel sequencing technology.

It is widely recognized that microhaps are powerful markers for different forensic purposes, mainly due to their advantages of both short tandem repeats and single nucleotide polymorphisms, including multiple alleles, low mutation rate, and absence of stutter peaks. In the present study, a panel of 60 microhap loci was developed and utilized in forensic kinship analysis as a preliminary study. Genotyping of microhap was performed by massively parallel sequencing and haplotypes were directly achieved from sequence reads of 73 samples from Chinese Han population. We observed that 49 out of 60 loci have effective number of alleles greater than 3.0 and 10 out of 60 have values above 4.0, with an average value of 3.5598. The heterozygosity values were in a range from 0.5840 to 0.8546 with an average of 0.7268 and the cumulative power of exclusion value of the 60 loci is equal to 1-4.78 × 10-18 . Moreover, we demonstrated the applicability of this method by different relationship inference problems, including identification of single parent-offspring, full-sibling, and second-degree relative. The results indicated that the assembled microhap panel provided more power for relationship inference, than commonly used short tandem repeats or single nucleotide polymorphism system.

Assigning forensic body fluids to donors in mixed body fluids by targeted RNA/DNA deep sequencing of coding region SNPs.

Biological traces found at crime scenes are analysed not only to genetically identify the donor(s) but also to determine the composition of the stain. For some cases, it is essential to associate a body fluid with a donor. Especially in mixed body fluid stains, but also in body fluid stains that appear to be single-source, this may be of importance. Linking a DNA profile (sub-source level) with evidence from a presumptive test or mRNA analysis (source level) is not straightforward. Our results support that associating donors and body fluids by means of comparing mixture ratios in RNA and DNA is not recommended. We introduce a set of 35 coding region SNPs (cSNPs) in body fluid-specific mRNA transcripts that represent a direct link between the body fluids and their donors. The discrimination power of the cSNPs was estimated based on allele frequencies calculated from a population sample (n = 188), and we investigated the practical application of the cSNPs in different scenarios. The results demonstrate that more cSNPs are needed to improve the discrimination power. However, the findings are promising as we were able to associate donors with body fluids in mixtures of different body fluids as well as in stains where both donors have contributed the same body fluid, e.g. a blood-blood mixture. In addition, the cSNP assay can be used for body fluid identification. The results of this proof-of-concept study support the use of cSNPs to assign body fluids to the respective donors.

Optimizing DNA recovery and forensic typing of degraded blood and dental remains using a specialized extraction method, comprehensive qPCR sample characterization, and massively parallel sequencing.

Human dental remains encountered in criminal casework evidence, missing person cases, or mass disaster tragedies provide a valuable sample source for DNA typing when suitable soft tissue is unavailable. Using traditional methods, teeth samples can be challenging to process, resulting in low-quantity and/or quality nuclear DNA and insufficient profiles for comparisons. This study examines the performance of a three-part nuclear DNA analysis workflow for teeth samples based on (1) improved dental tissue recovery using the Dental Forensic Kit (DFKMR) (Universidad de los Andes) and DNA extraction with QuickExtract™ FFPE DNA Extraction Kit (Lucigen®), (2) quantification with InnoQuant® HY (InnoGenomics Technologies) for sensitive assessment of total human and male DNA quantity/quality, and (3) massively parallel sequencing for simultaneous genotyping of 231 short tandem repeat (STR) and single-nucleotide polymorphism (SNP) markers with the ForenSeq® DNA Signature Prep Kit (Verogen, Inc.). Initial evaluation of artificially degraded blood samples (n = 10) achieved highly sensitive and informative quantification results with InnoQuant® HY, enabling successful first pass genotyping with the MiSeq FGx® System. Twenty-three STR alleles (out of 85) and 70 identity informative SNP loci (out of 94) were recovered from two pg total long target DNA input (0.86 ng total short target input) and an InnoQuant degradation index (DI) of 460 (severely degraded). The three-part workflow was subsequently applied to teeth samples (dental pulp, root cement tissues; n = 13) with postmortem intervals (PMI) of the teeth ranging from 8 days to approximately 6 months. Informative SNP and STR DNA profiles were obtained, to include 78 STR alleles and 85 identity informative SNP loci typed (of 94 total SNP targets) in a 1 month, four-day PMI root cement sample with one pg total long target DNA input and a DI of 76. These data indicate successful performance of the proposed workflow from degraded DNA from teeth samples.

Mixture deconvolution by massively parallel sequencing of microhaplotypes.

Short tandem repeat polymorphisms (STRs) are the standard markers for forensic human identification. STRs are highly polymorphic loci analyzed using a direct PCR-to-CE (capillary electrophoresis) approach. However, STRs have limitations particularly when dealing with complex mixtures. These include slippage of the polymerase during amplification causing stutter fragments that can be indistinguishable from minor contributor alleles, preferential amplification of shorter alleles, and limited number of loci that can be effectively co-amplified with CE. Massively parallel sequencing (MPS), by enabling a higher level of multiplexing and actual sequencing of the DNA, provides forensic practitioners an increased power of discrimination offered by the sequence of STR alleles and access to new sequence-based markers. Microhaplotypes (i.e., microhaps or MHs) are emerging multi-allelic loci of two or more SNPs within < 300 bp that are highly polymorphic, have alleles all of the same length, and do not generate stutter fragments. The growing number of loci described in the literature along with initial mixture investigations supports the potential for microhaps to aid in mixture interpretation and the purpose of this study was to demonstrate that practically. A panel of 36 microhaplotypes, selected from a set of over 130 loci, was tested with the Ion S5™ MPS platform (Thermo Fisher Scientific) on single-source samples, synthetic two-to-six person mixtures at different concentrations/contributor ratios, and on crime scene-like samples. The panel was tested both in multiplex with STRs and SNPs and individually. The analysis of single-source samples showed that the allele coverage ratio across all loci was 0.88 ± 0.08 which is in line with the peak height ratio of STR alleles in CE. In mixture studies, results showed that the input DNA can be much higher than with conventional CE, without the risk of oversaturating the detection system, enabling an increased sensitivity for the minor contributor in imbalanced mixtures with abundant amounts of DNA. Furthermore, the absence of stutter fragments simplifies the interpretation. On casework-like samples, MPS of MHs enabled the detection of a higher number of alleles from minor donors than MPS and CE of STRs. These results demonstrated that MPS of microhaplotypes can complement STRs and enhance human identification practices when dealing with complex imbalanced mixtures.

TIA: algorithms for development of identity-linked SNP islands for analysis by massively parallel DNA sequencing.

BACKGROUND: Single nucleotide polymorphisms (SNPs) located within the human genome have been shown to have utility as markers of identity in the differentiation of DNA from individual contributors. Massively parallel DNA sequencing (MPS) technologies and human genome SNP databases allow for the design of suites of identity-linked target regions, amenable to sequencing in a multiplexed and massively parallel manner. Therefore, tools are needed for leveraging the genotypic information found within SNP databases for the discovery of genomic targets that can be evaluated on MPS platforms. RESULTS: The SNP island target identification algorithm (TIA) was developed as a user-tunable system to leverage SNP information within databases. Using data within the 1000 Genomes Project SNP database, human genome regions were identified that contain globally ubiquitous identity-linked SNPs and that were responsive to targeted resequencing on MPS platforms. Algorithmic filters were used to exclude target regions that did not conform to user-tunable SNP island target characteristics. To validate the accuracy of TIA for discovering these identity-linked SNP islands within the human genome, SNP island target regions were amplified from 70 contributor genomic DNA samples using the polymerase chain reaction. Multiplexed amplicons were sequenced using the Illumina MiSeq platform, and the resulting sequences were analyzed for SNP variations. 166 putative identity-linked SNPs were targeted in the identified genomic regions. Of the 309 SNPs that provided discerning power across individual SNP profiles, 74 previously undefined SNPs were identified during evaluation of targets from individual genomes. Overall, DNA samples of 70 individuals were uniquely identified using a subset of the suite of identity-linked SNP islands. CONCLUSIONS: TIA offers a tunable genome search tool for the discovery of targeted genomic regions that are scalable in the population frequency and numbers of SNPs contained within the SNP island regions. It also allows the definition of sequence length and sequence variability of the target region as well as the less variable flanking regions for tailoring to MPS platforms. As shown in this study, TIA can be used to discover identity-linked SNP islands within the human genome, useful for differentiating individuals by targeted resequencing on MPS technologies.

Ancestry inference of 96 population samples using microhaplotypes.

Microhaplotypes have become a new type of forensic marker with a great ability to identify and deconvolute mixtures because massively parallel sequencing (MPS) allows the alleles (haplotypes) of the multi-SNP loci to be determined directly for an individual. As originally defined, a microhaplotype locus is a short segment of DNA with two or more SNPs defining three or more haplotypes. The length is short enough, less than about 300 bp, that the read length of current MPS technology can produce a phase-known sequence of each chromosome of an individual. As part of the discovery phase of our studies, data on 130 microhaplotype loci with estimates of haplotype frequency data on 83 populations have been published. To provide a better picture of global allele frequency variation, we have now tested 13 more populations for 65 of the microhaplotype loci from among those with higher levels of inter-population gene frequency variation, including 8 loci not previously published. These loci provide clear distinctions among 6 biogeographic regions and provide some information distinguishing up to 10 clusters of populations.

An SNP panel for the analysis of paternally inherited alleles in maternal plasma using ion Torrent PGM.

Researchers have sought to develop an effective protocol for paternity analysis using cell-free DNA (cfDNA) in maternal plasma. The use of massively parallel sequencing (MPS) technology for SNP testing is attractive because of its high-throughput capacity and resolution to single-base precision. In this study, we designed a customized SNP panel for cfDNA sequencing that includes 720 short amplicons (< 140 bp) targeting SNPs on the autosome and Y chromosome. The systemic performance was evaluated using the Ion Torrent PGM, indicating balanced coverage among most of the included loci, except for 78 poorly performing SNPs that were observed to have an inconsistent allele balance, lower coverage reads or high background signals. Then, the custom panel was used to perform cfDNA genotyping in maternal plasma from 20 pregnancies in the first and second trimesters (9 to 21 weeks). By establishing an allele fraction cutoff of 2.0%, 53 to 128 autosomal SNP loci were considered informative for paternal origin. Validation results in foetal samples showed that 49.43% to 100% of the real paternal alleles were accurately identified, with incorrect alleles encountered in 3 cases. The concentration of foetal cfDNA ranged from 4.28% to 10.70%. Our results show that this amplicon-based sequencing strategy could be utilized in analysing paternally inherited alleles in maternal plasma. However, further studies and optimization are required for a more detailed and accurate interpretation of the cfDNA sequencing results based on MPS technology.

Comparison between magnetic bead and qPCR library normalisation methods for forensic MPS genotyping.

Massively parallel sequencing (MPS) is fast approaching operational use in forensic science, with the capability to analyse hundreds of DNA identity and DNA intelligence markers in multiple samples simultaneously. The ForenSeq™ DNA Signature Kit on MiSeq FGx™ (Illumina) workflow can provide profiles for autosomal short tandem repeats (STRs), X chromosome and Y chromosome STRs, identity single nucleotide polymorphisms (SNPs), biogeographical ancestry SNPs and phenotype (eye and hair colour) SNPs from a sample. The library preparation procedure involves a series of steps including target amplification, library purification and library normalisation. This study highlights the comparison between the manufacturer recommended magnetic bead normalisation and quantitative polymerase chain reaction (qPCR) methods. Furthermore, two qPCR chemistries, KAPA® (KAPA Biosystems) and NEBNext® (New England Bio Inc.), have also been compared. The qPCR outperformed the bead normalisation method, while the NEBNext® kit obtained higher genotype concordance than KAPA®. The study also established an MPS workflow that can be utilised in any operational forensic laboratory.

More comprehensive forensic genetic marker analyses for accurate human remains identification using massively parallel DNA sequencing.

BACKGROUND: Although the primary objective of forensic DNA analyses of unidentified human remains is positive identification, cases involving historical or archaeological skeletal remains often lack reference samples for comparison. Massively parallel sequencing (MPS) offers an opportunity to provide biometric data in such cases, and these cases provide valuable data on the feasibility of applying MPS for characterization of modern forensic casework samples. In this study, MPS was used to characterize 140-year-old human skeletal remains discovered at a historical site in Deadwood, South Dakota, United States. The remains were in an unmarked grave and there were no records or other metadata available regarding the identity of the individual. Due to the high throughput of MPS, a variety of biometric markers could be typed using a single sample. RESULTS: Using MPS and suitable forensic genetic markers, more relevant information could be obtained from a limited quantity and quality sample. Results were obtained for 25/26 Y-STRs, 34/34 Y SNPs, 166/166 ancestry-informative SNPs, 24/24 phenotype-informative SNPs, 102/102 human identity SNPs, 27/29 autosomal STRs (plus amelogenin), and 4/8 X-STRs (as well as ten regions of mtDNA). The Y-chromosome (Y-STR, Y-SNP) and mtDNA profiles of the unidentified skeletal remains are consistent with the R1b and H1 haplogroups, respectively. Both of these haplogroups are the most common haplogroups in Western Europe. Ancestry-informative SNP analysis also supported European ancestry. The genetic results are consistent with anthropological findings that the remains belong to a male of European ancestry (Caucasian). Phenotype-informative SNP data provided strong support that the individual had light red hair and brown eyes. CONCLUSIONS: This study is among the first to genetically characterize historical human remains with forensic genetic marker kits specifically designed for MPS. The outcome demonstrates that substantially more genetic information can be obtained from the same initial quantities of DNA as that of current CE-based analyses.

Massively Parallel Sequencing Detected a Mutation in the MFN2 Gene Missed by Sanger Sequencing Due to a Primer Mismatch on an SNP Site.

We describe a patient with early onset severe axonal Charcot-Marie-Tooth disease (CMT2) with dominant inheritance, in whom Sanger sequencing failed to detect a mutation in the mitofusin 2 (MFN2) gene because of a single nucleotide polymorphism (rs2236057) under the PCR primer sequence. The severe early onset phenotype and the family history with severely affected mother (died after delivery) was very suggestive of CMT2A and this suspicion was finally confirmed by a MFN2 mutation. The mutation p.His361Tyr was later detected in the patient by massively parallel sequencing with a gene panel for hereditary neuropathies. According to this information, new primers for amplification and sequencing were designed which bind away from the polymorphic sites of the patient's DNA. Sanger sequencing with these new primers then confirmed the heterozygous mutation in the MFN2 gene in this patient. This case report shows that massively parallel sequencing may in some rare cases be more sensitive than Sanger sequencing and highlights the importance of accurate primer design which requires special attention.

Implication of genomic characterization in synchronous endometrial and ovarian cancers of endometrioid histology.

OBJECTIVES: Synchronous endometrial and ovarian carcinomas (SEOCs) present gynecologic oncologists with a challenging diagnostic puzzle: discriminating between double primary cancers and single primary cancer with metastasis. We aimed to determine the clonal relationship between simultaneously diagnosed endometrial and ovarian carcinomas. METHODS: Fourteen pairs of SEOCs of endometrioid type and two pairs of SEOCs with disparate histologic types (control for dual primary tumors) were subjected to massively parallel sequencing (MPS) and molecular inversion probe microarrays. RESULTS: Thirteen of the 14 pairs of SEOCs harbored somatic mutations shared by both uterine and ovarian lesions, indicative of clonality. High degree of chromosomal instability in the tumors from 10 patients who received adjuvant chemotherapy, of whom 9 had synchronous carcinomas with significantly overlapping copy number alterations (CNAs), suggestive of single primary tumors with metastasis. The clonal relationship determined by genomic analyses did not agree with clinicopathological criteria in 11 of 14 cases. Minimal CNAs were identified in both ovarian and endometrial carcinomas in 4 patients, who did not receive adjuvant chemotherapy and experienced no recurrent diseases. In contrast, two of the 10 patients with chromosomally unstable cancers developed recurrent tumors. CONCLUSION: Our findings support a recent paradigm-shifting concept that most SEOCs originate from a single tumor. It also casts doubt on the clinicopathological criteria used to distinguish between dual primary tumors and single primary tumor with metastasis. Testing of CNAs on SEOCs may help determining the need of adjuvant therapy.

Effects of the Ion PGM™ Hi-Q™ sequencing chemistry on sequence data quality.

Massively parallel sequencing (MPS) offers substantial improvements over current forensic DNA typing methodologies such as increased resolution, scalability, and throughput. The Ion PGM™ is a promising MPS platform for analysis of forensic biological evidence. The system employs a sequencing-by-synthesis chemistry on a semiconductor chip that measures a pH change due to the release of hydrogen ions as nucleotides are incorporated into the growing DNA strands. However, implementation of MPS into forensic laboratories requires a robust chemistry. Ion Torrent's Hi-Q™ Sequencing Chemistry was evaluated to determine if it could improve on the quality of the generated sequence data in association with selected genetic marker targets. The whole mitochondrial genome and the HID-Ion STR 10-plex panel were sequenced on the Ion PGM™ system with the Ion PGM™ Sequencing 400 Kit and the Ion PGM™ Hi-Q™ Sequencing Kit. Concordance, coverage, strand balance, noise, and deletion ratios were assessed in evaluating the performance of the Ion PGM™ Hi-Q™ Sequencing Kit. The results indicate that reliable, accurate data are generated and that sequencing through homopolymeric regions can be improved with the use of Ion Torrent's Hi-Q™ Sequencing Chemistry. Overall, the quality of the generated sequencing data supports the potential for use of the Ion PGM™ in forensic genetic laboratories.

Whole mitochondrial genome genetic diversity in an Estonian population sample.

Mitochondrial DNA is a useful marker for population studies, human identification, and forensic analysis. Commonly used hypervariable regions I and II (HVI/HVII) were reported to contain as little as 25% of mitochondrial DNA variants and therefore the majority of power of discrimination of mitochondrial DNA resides in the coding region. Massively parallel sequencing technology enables entire mitochondrial genome sequencing. In this study, buccal swabs were collected from 114 unrelated Estonians and whole mitochondrial genome sequences were generated using the Illumina MiSeq system. The results are concordant with previous mtDNA control region reports of high haplogroup HV and U frequencies (47.4 and 23.7% in this study, respectively) in the Estonian population. One sample with the Northern Asian haplogroup D was detected. The genetic diversity of the Estonian population sample was estimated to be 99.67 and 95.85%, for mtGenome and HVI/HVII data, respectively. The random match probability for mtGenome data was 1.20 versus 4.99% for HVI/HVII. The nucleotide mean pairwise difference was 27 ± 11 for mtGenome and 7 ± 3 for HVI/HVII data. These data describe the genetic diversity of the Estonian population sample and emphasize the power of discrimination of the entire mitochondrial genome over the hypervariable regions.

Differentiating between monozygotic twins through next-generation mitochondrial genome sequencing.

Monozygotic (MZ) twins, considered to be genetically identical, cannot be distinguished from one another by standard forensic DNA testing. A recent study employed whole genome sequencing to identify extremely rare mutations and reported that mutation analysis could be used to differentiate between MZ twins. Compared with nuclear DNA, mitochondrial DNA (mtDNA) has higher mutation rates; therefore, minor differences theoretically exist in MZ twins' mitochondrial genome (mtGenome). However, conventional Sanger-type sequencing (STS) is neither amenable to, nor feasible for, the detection of low-level sequence variants. The recent introduction of massively parallel sequencing (MPS) has the capability to sequence many targeted regions of multiple samples simultaneously with desirable depth of coverage. Thus, the aim of this study was to assess whether full mtGenome sequencing analysis can be used to differentiate between MZ twins. Ten sets of MZ twins provided blood samples that underwent extraction, quantification, mtDNA enrichment, library preparation, and ultra-deep sequencing. Point heteroplasmies were observed in eight sets of MZ twins, and a single nucleotide variant (nt15301) was detected in five sets of MZ twins. Thus, this study demonstrates that ultra-deep mtGenome sequencing could be used to differentiate between MZ twins.

Optimization of virus detection in cells using massively parallel sequencing.

Massively parallel sequencing (MPS)-based virus detection has potential regulatory applications. We studied the ability of one of these approaches, based on degenerate oligonucleotide primer (DOP)-polymerase chain reaction (PCR), to detect viral sequences in cell lines known to express viral genes or particles. DOP-PCR was highly sensitive for the detection of small quantities of isolated viral sequences. Detected viral sequences included nodavirus, bracovirus, and endogenous retroviruses in High Five cells, porcine circovirus type 1 and porcine endogenous retrovirus in PK15 cells, human T-cell leukemia virus 1 in MJ cells, human papillomavirus 18 in HeLa cells, human herpesvirus 8 in BCBL-1 cells, and Epstein-Barr Virus in Raji cells. Illumina sequencing (for which primers were most efficiently added using PCR) provided greater sensitivity for virus detection than Roche 454 sequencing. Analyzing nucleic acids extracted both directly from samples and from capsid-enriched preparations provided useful information. Although there are limitations of these methods, these results indicate significant promise for the combination of nonspecific PCR and MPS in identifying contaminants in clinical and biological samples, including cell lines and reagents used to produce vaccines and therapeutic products.