Locus-specific Sequencing Research Articles

Genetic Diversity of the Species of the Genus Deschampsia P.Beauv. (Poaceae) Based on the Analysis of the ITS Region: Polymorphism Proves Distant Hybridization.

The species of the genus Deschampsia are difficult for identification, and the genus is difficult for taxonomic treatment. The regions of 35S rRNA genes were studied for the species of the genus Deschampsia of different geographical origin with a method of sequencing by Sanger (ITS1-5.8S rRNA gene-ITS2, 14 species) and with a method of a locus-specific next-generation sequencing (NGS) on the Illumina platform (ITS1-5.8S rRNA, 7 species). All species of Deschampsia formed one clade; the species, referred by some authors on the basis of morphological characters to the species D. cespitosa s.l., entered one subclade. Subantarctic species formed a separate subclade and their ribotypes formed their own subnetwork. Avenella flexuosa, earlier referred to Deschampsia, entered the other clade, though this species contains some ribotypes common with some Deschampsia species. Deschampsia pamirica and related mountain species have their own specific ribotype groups. On the network of the ribotypes, one can see that D. cespitosa from Great Britain forms a network with some species, but D. cespitosa from the USA forms its own network. Ribotype analysis of each sample revealed traces of introgression with Deyeuxia/Calamagrostis in D. cespitosa and with A. flexuosa and probable introgression of Northern and subantarctic species.

Personalized recurrence risk assessment following the birth of a child with a pathogenic de novo mutation

Following the diagnosis of a paediatric disorder caused by an apparently de novo mutation, a recurrence risk of 1–2% is frequently quoted due to the possibility of parental germline mosaicism; but for any specific couple, this figure is usually incorrect. We present a systematic approach to providing individualized recurrence risk. By combining locus-specific sequencing of multiple tissues to detect occult mosaicism with long-read sequencing to determine the parent-of-origin of the mutation, we show that we can stratify the majority of couples into one of seven discrete categories associated with substantially different risks to future offspring. Among 58 families with a single affected offspring (representing 59 de novo mutations in 49 genes), the recurrence risk for 35 (59%) was decreased below 0.1%, but increased owing to parental mixed mosaicism for 5 (9%)—that could be quantified in semen for paternal cases (recurrence risks of 5.6–12.1%). Implementation of this strategy offers the prospect of driving a major transformation in the practice of genetic counselling.

HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 allele frequencies of North Tanzanian Maasai

Locus-specific amplicon sequencing was used to HLA type 336 participants of Maasai ethnicity at the HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 loci. Participants were recruited from three study villages in North Tanzania, for the purpose of investigating risk factors for trachomatous scarring in children. Other than HLA-A, all loci significantly deviated from Hardy-Weinberg equilibrium, possibly due to high relatedness between individuals: 238 individuals shared a house with at least one another participant. The most frequent allele for each locus were A*68:02 (14.3 %), B*53:01 (8.4 %), C*06:02 (19.2 %), DRB1*13:02 (17.7 %), DQB1*02:01 (16.9 %) and DPB1*01:01 (15.7 %), while the most common inferred haplotype was A*68:02 ∼ B*18:01 ∼ C*07:04 ∼ DRB1*08:04 ∼ DQB1*04:02 ∼ DPB1*04:01 (1.3 %).

Sequence based HLA-DRB1, -DQB1 and -DPB1 allele and haplotype frequencies in The Gambia

Class II HLA loci DRB1, DQB1 and DPB1 were typed for a total of 939 Gambian participants by locus-specific amplicon sequencing. Participants were from multiple regions of The Gambia and drawn from two studies: a family study aiming to identify associations between host genotype and trachomatous scarring (N = 796) and a cohort study aiming to identify correlates of immunity to trachoma (N = 143). All loci deviated from Hardy-Weinberg equilibrium, likely due to the family-based nature of the study: 608 participants had at least one other family member included in the study population. The most common alleles for HLA-DRB1, DQB1 and DPB1 respectively were DRB1*13:04 (18.8 %), DQB1*03:19 (27.9 %) and DPB1*01:01 (25.4 %). Participants belonged to a variety of ethnicities, including the Mandinka, Fula, Wolof and Jola ethnic groups.

CYP2D6 gene resequencing in the Malagasy, a population at the crossroads between Asia and Africa: a pilot study.

Background:Plasmodium vivax malaria is endemic in Madagascar, where populations have genetic inheritance from Southeast Asia and East Africa. Primaquine, a drug of choice for vivax malaria, is metabolized principally via CYP2D6. CYP2D6 variation was characterized by locus-specific gene sequencing and was compared with TaqMan™ genotype data.Materials & methods:Long-range PCR amplicons were generated from 96 Malagasy samples and subjected to next-generation sequencing.Results:The authors observed high concordance between TaqMan™-based CYP2D6 genotype calls and the base calls from sequencing. In addition, there are new variants and haplotypes present in the Malagasy.Conclusion:Sequencing unique admixed populations provides more detailed and accurate insights regarding CYP2D6 variability, which may help optimize primaquine treatment across human genetic diversity.

Rats exhibit age-related mosaic loss of chromosome Y

Mosaic loss of the Y chromosome (LOY) is the most frequent chromosomal aberration in aging men and is strongly correlated with mortality and disease. To date, studies of LOY have only been performed in humans, and so it is unclear whether LOY is a natural consequence of our relatively long lifespan or due to exposure to human-specific external stressors. Here, we explored whether LOY could be detected in rats. We applied a locus-specific PCR and target sequencing approach that we used as a proxy to estimate LOY in 339 samples covering eleven tissues from young and old individuals. We detected LOY in four tissues of older rats. To confirm the results from the PCR screening, we re-sequenced 60 full genomes from old rats, which revealed that the Y chromosome is the sole chromosome with low copy numbers. Finally, our results suggest that LOY is associated with other structural aberrations on the Y chromosome and possibly linked to the mosaic loss of the X chromosome. This is the first report, to our knowledge, demonstrating that the patterns of LOY observed in aging men are also present in a rodent, and conclude that LOY may be a natural process in placental mammals.

Genomic constitution of diploid species of the genus Avena L. with A-genome according to the data of the next-generation sequencing (NGS)

For diploid (2x) species with the A-genome, as well as for hexaploid (6x) from the genus Avena, a locus-specific next-generation sequencing (NGS) of the sequence of the region of the internal transcribed spacer ITS1 and the beginning of the 5.8S rRNA gene was carried out on the Illumina platform. The high diversity and heterogeneity of the genomes of diploid species are shown. It was revealed that the genomes of modern diploid oat species are relatively far removed from the hexaploid species. It was found that A. canariensis occupies an isolated position among other diploid species, and also takes only an insignificant role in the formation of hexaploid genomes.

Mutation Associated with Orange Fruit Color Increases Concentrations of β-Carotene in a Sweet Pepper Variety (Capsicumannuum L.).

Pepper is the second most important vegetable crop in Bulgarian agriculture and has become the subject of extensive breeding programs that frequently employ induced mutagenesis. The success of breeding programs can be enhanced by the efficient and integral application of different biochemical and molecular methods to characterize specific mutant alleles. On the other hand, identifying new cost-effective methods is important under a limited-resources environment. In this paper we compare the levels of five health-related carotenoid compounds of fruits (α-carotene, β-carotene, lutein, β-cryptoxanthin, zeaxanthin) between a mutant variety Oranzheva kapia (possessing high ß-carotene concentration) and a corresponding initial pepper variety Pazardzhishka kapia 794. Both varieties are intended for fresh consumption. Pepper is a major natural source of β-carotene. It was observed that fruit at both commercial and botanical maturity from mutant variety had greater α-carotene and β-carotene concentrations to the initial variety (7.49 and 1.94 times higher, respectively) meaning that the mutant was superior in fruit quality to the initial genotype. Two hydroxylase enzymes, converting α- and β-carotene to lutein and zeaxanthin, respectively, are known to exist in pepper and are encoded by two genes on chromosomes 3 and 6-CrtZchr03 and CrtZchr06. The molecular characterization of the mutant variety through locus-specific Polymerase chain reaction amplification, gene cloning and sequencing as well as expression was performed. Our results suggest that the increased ß-carotene accumulation in the mutant variety Oranzheva kapia results from a biosynthetic pathway breakdown due to deletion of CrtZchr03 gene.

LINE-1 Evasion of Epigenetic Repression in Humans.

Epigenetic silencing defends against LINE-1 (L1) retrotransposition in mammalian cells. However, the mechanisms that repress young L1 families and how L1 escapes to cause somatic genome mosaicism in the brain remain unclear. Here we report that a conserved Yin Yang 1 (YY1) transcription factor binding site mediates L1 promoter DNA methylation in pluripotent and differentiated cells. By analyzing 24 hippocampal neurons with three distinct single-cell genomic approaches, we characterized and validated a somatic L1 insertion bearing a 3' transduction. The source (donor) L1 for this insertion was slightly 5' truncated, lacked the YY1 binding site, and was highly mobile when tested invitro. Locus-specific bisulfite sequencing revealed that the donor L1 and other young L1s with mutated YY1 binding sites were hypomethylated in embryonic stem cells, during neurodifferentiation, and in liver and brain tissue. These results explain how L1 can evade repression and retrotranspose in the human body.

Targets and genomic constraints of ectopic Dnmt3b expression.

DNA methylation plays an essential role in mammalian genomes and expression of the responsible enzymes is tightly controlled. Deregulation of the de novo DNA methyltransferase DNMT3B is frequently observed across cancer types, yet little is known about its ectopic genomic targets. Here, we used an inducible transgenic mouse model to delineate rules for abnormal DNMT3B targeting, as well as the constraints of its activity across different cell types. Our results explain the preferential susceptibility of certain CpG islands to aberrant methylation and point to transcriptional state and the associated chromatin landscape as the strongest predictors. Although DNA methylation and H3K27me3 are usually non-overlapping at CpG islands, H3K27me3 can transiently co-occur with DNMT3B-induced DNA methylation. Our genome-wide data combined with ultra-deep locus-specific bisulfite sequencing suggest a distributive activity of ectopically expressed Dnmt3b that leads to discordant CpG island hypermethylation and provides new insights for interpreting the cancer methylome.

High-density genetic map construction and QTLs identification for plant height in white jute (Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing

BackgroundGenetic mapping and quantitative trait locus (QTL) detection are powerful methodologies in plant improvement and breeding. White jute (Corchorus capsularis L.) is an important industrial raw material fiber crop because of its elite characteristics. However, construction of a high-density genetic map and identification of QTLs has been limited in white jute due to a lack of sufficient molecular markers. The specific locus amplified fragment sequencing (SLAF-seq) strategy combines locus-specific amplification and high-throughput sequencing to carry out de novo single nuclear polymorphism (SNP) discovery and large-scale genotyping. In this study, SLAF-seq was employed to obtain sufficient markers to construct a high-density genetic map for white jute. Moreover, with the development of abundant markers, genetic dissection of fiber yield traits such as plant height was also possible. Here, we present QTLs associated with plant height that were identified using our newly constructed genetic linkage groups.ResultsAn F8 population consisting of 100 lines was developed. In total, 69,446 high-quality SLAFs were detected of which 5,074 SLAFs were polymorphic; 913 polymorphic markers were used for the construction of a genetic map. The average coverage for each SLAF marker was 43-fold in the parents, and 9.8-fold in each F8 individual. A linkage map was constructed that contained 913 SLAFs on 11 linkage groups (LGs) covering 1621.4 cM with an average density of 1.61 cM per locus. Among the 11 LGs, LG1 was the largest with 210 markers, a length of 406.34 cM, and an average distance of 1.93 cM between adjacent markers. LG11 was the smallest with only 25 markers, a length of 29.66 cM, and an average distance of 1.19 cM between adjacent markers. ‘SNP_only’ markers accounted for 85.54% and were the predominant markers on the map. QTL mapping based on the F8 phenotypes detected 11 plant height QTLs including one major effect QTL across two cultivation locations, with each QTL accounting for 4.14–15.63% of the phenotypic variance.ConclusionsTo our knowledge, the linkage map constructed here is the densest one available to date for white jute. This analysis also identified the first QTL in white jute. The results will provide an important platform for gene/QTL mapping, sequence assembly, genome comparisons, and marker-assisted selection breeding for white jute.

Y chromosome haplotype diversity in Mongolic-speaking populations and gene conversion at the duplicated STR DYS385a,b in haplogroup C3-M407.

Y chromosome microsatellite (Y-STR) diversity has been studied in different Mongolic-speaking populations from South Siberia, Mongolia, North-East China and East Europe. The results obtained indicate that the Mongolic-speaking populations clustered into two groups, with one group including populations from eastern part of South Siberia and Central Asia (the Buryats, Barghuts and Khamnigans) and the other group including populations from western part of Central Asia and East Europe (the Mongols and Kalmyks). High frequency of haplogroup C3-M407 (>50%) is present in the Buryats, Barghuts and Khamnigans, whereas in the Mongols and Kalmyks its frequency is much lower. In addition, two allelic combinations in DYS385a,b loci of C3-M407 haplotypes have been observed: the combination 11,18 (as well as 11,17 and 11,19) is frequent in different Mongolic-speaking populations, but the 11,11 branch is present mainly in the Kalmyks and Mongols. Results of locus-specific sequencing suggest that the action of gene conversion is a more likely explanation for origin of homoallelic 11,11 combination. Moreover, analysis of median networks of Y-STR haplotypes demonstrates that at least two gene conversion events can be revealed-one of them has probably occurred among the Mongols, and the other event occurred in the Barghuts. These two events give an average gene conversion rate range of 0.24-7.1 × 10(-3) per generation.

Single-cell, locus-specific bisulfite sequencing (SLBS) for direct detection of epimutations in DNA methylation patterns.

Stochastic epigenetic changes drive biological processes, such as development, aging and disease. Yet, epigenetic information is typically collected from millions of cells, thereby precluding a more precise understanding of cell-to-cell variability and the pathogenic history of epimutations. Here we present a novel procedure for directly detecting epimutations in DNA methylation patterns using single-cell, locus-specific bisulfite sequencing (SLBS). We show that within gene promoter regions of mouse hepatocytes the epimutation rate is two orders of magnitude higher than the mutation rate.

Estimating absolute methylation levels at single-CpG resolution from methylation enrichment and restriction enzyme sequencing methods

Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to map complete DNA methylomes. These include whole-genome bisulfite sequencing (WGBS, MethylC-seq, or BS-seq), reduced-representation bisulfite sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. These methods yield largely comparable results but differ significantly in extent of genomic CpG coverage, resolution, quantitative accuracy, and cost, at least while using current algorithms to interrogate the data. None of these existing methods provides single-CpG resolution, comprehensive genome-wide coverage, and cost feasibility for a typical laboratory. We introduce methylCRF, a novel conditional random fields–based algorithm that integrates methylated DNA immunoprecipitation (MeDIP-seq) and methylation-sensitive restriction enzyme (MRE-seq) sequencing data to predict DNA methylation levels at single-CpG resolution. Our method is a combined computational and experimental strategy to produce DNA methylomes of all 28 million CpGs in the human genome for a fraction (<10%) of the cost of whole-genome bisulfite sequencing methods. methylCRF was benchmarked for accuracy against Infinium arrays, RRBS, WGBS sequencing, and locus-specific bisulfite sequencing performed on the same human embryonic stem cell line. methylCRF transformation of MeDIP-seq/MRE-seq was equivalent to a biological replicate of WGBS in quantification, coverage, and resolution. We used conventional bisulfite conversion, PCR, cloning, and sequencing to validate loci where our predictions do not agree with whole-genome bisulfite data, and in 11 out of 12 cases, methylCRF predictions of methylation level agree better with validated results than does whole-genome bisulfite sequencing. Therefore, methylCRF transformation of MeDIP-seq/MRE-seq data provides an accurate, inexpensive, and widely accessible strategy to create full DNA methylomes.

DMEAS: DNA methylation entropy analysis software

Summary: DMEAS is the first user-friendly tool dedicated to analyze the distribution of DNA methylation patterns for the quantification of epigenetic heterogeneity. It supports the analysis of both locus-specific and genome-wide bisulfite sequencing data. DMEAS progressively scans the mapping results of bisulfite sequencing reads to extract DNA methylation patterns for contiguous CpG dinucleotides. It determines the DNA methylation level and calculates methylation entropy for genomic segments to enable the quantitative assessment of DNA methylation variations observed in cell populations.Availability and implementation: DMEAS program, user guide and all the testing data are freely available from http://sourceforge.net/projects/dmeas/files/Contact: davidxie@vt.eduSupplementary Information: Supplementary data are available at Bioinformatics online.

Tet-assisted bisulfite sequencing of 5-hydroxymethylcytosine

A complete understanding of the potential function of 5-hydroxymethylcytosine (5-hmC), a DNA cytosine modification in mammalian cells, requires an accurate single-base resolution sequencing method. Here we describe a modified bisulfite-sequencing method, Tet-assisted bisulfite sequencing (TAB-seq), which can identify 5-hmC at single-base resolution, as well as determine its abundance at each modification site. This protocol involves β-glucosyltransferase (β-GT)-mediated protection of 5-hmC (glucosylation) and recombinant mouse Tet1(mTet1)-mediated oxidation of 5-methylcytosine (5-mC) to 5-carboxylcytosine (5-caC). After the subsequent bisulfite treatment and PCR amplification, both cytosine and 5-caC (derived from 5-mC) are converted to thymine (T), whereas 5-hmC reads as C. The treated genomic DNA is suitable for both whole-genome and locus-specific sequencing. The entire procedure (which does not include data analysis) can be completed in 14 d for whole-genome sequencing or 7 d for locus-specific sequencing.

Abstract 4647: Clinical Implications of Rabphillin-3A-Like (RPH3AL/NOC2) Gene Alterations in Breast Cancer

Abstract Background: The Rabphillin-3A-like (RPH3AL) gene, located at the 17p13.3 locus, has been reported to be involved in progression and metastasis of colorectal cancer. Since there have been no studies of this gene in breast cancer (BC), mutations and loss of heterozygosity (LOH) of RPH3AL in BCs were evaluated. Methods: BCs (n=104) and matching benign tissues were evaluated for LOH by use of two locus-specific polymorphic markers (D17S1866 and D17S643) and ABI-3100 sequencing technology. Of 104 cases, 16 had corresponding snap-frozen BCs and benign epithelial tissues. The cDNA samples of these frozen cases were analyzed for RPH3AL mutations by RT-PCR and direct sequencing. Since most of the missense mutations were located in the exon-6 region, the genomic DNA samples of 68 (of 104) BCs and their corresponding benign tissues were assessed only for exon 6 by the use of exon-6 specific primers. Additionally, 96 (of 104) BCs and the corresponding control tissues were analyzed for germ line mutations at 43 base pairs downstream from exon-6. Mutations and LOH status were correlated with clinicopathologic features of BCs and patient overall survival by chi-square and Kaplan-Meier analyses, respectively. Results: LOH analysis demonstrated that 69 of 99 (70%) informative cases exhibited LOH in at least one marker. In total, 19% (13 of 68) of BCs exhibited missense point mutations in the coding region (exon-6) of RPH3AL. Six BCs exhibited mutations at codon 175, two BCs at codon 200, one BC had double mutations at codon 178 and 189, and four BCs had single mutations at codons 165, 174, 181 and 196. There was a germ-line mutation (G&amp;gt;T) at 43 base pairs downstream from exon-6. Genotyping of this alteration showed that 31% (30 of 96) of BCs were homozygous for the G/G allele, 56% (54 of 96) were heterozygous for G/T, and 13% (12 of 96) were homozygous for T/T. In our correlative analysis, most BCs with mutations were associated with lymph-node metastases (9 of 13, 69%) and with an allelic loss (8 of 13, 62%). The G/G alleles were more commonly found in African-Americans than Caucasians (χ2, p=0.003), in ER/PR positive BCs than in ER/PR negative BCs (χ2, p=0.03), and in BCs with advanced stage (χ2, p=0.02). The incidence of LOH increased with nodal metastasis (χ2, p=0.04) and tumor size (χ2, p=0.05) of the BCs. Furthermore, LOH was significantly associated with poor patient survival (log rank, p=0.002). Conclusions: The results suggest that genetic alterations in the RPH3AL gene are associated with aggressive behavior of breast cancer. (Grant Support: UAB-Breast SPORE-Pilot and Susan G. Komen Breast Cancer Foundation, POP138306). Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 4647.

Imported bancroftian filariasis: Diethylcarbamazine response and benzimidazole susceptibility of Wuchereria bancrofti in dynamic cross-border migrant population targeted by the National Program to Eliminate Lymphatic Filariasis in South Thailand

The implementation on the Thailand–Myanmar border of annual mass drug administration (MDA) of a single 6 mg/kg dose of diethylcarbamazine (DEC) plus 400 mg albendazole, part of the National Program to Eliminate Lymphatic Filariasis (PELF), has been challenging. In particular, chain migration of cross-border Myanmar workers at risk for nocturnally periodic Wuchereria bancrofti infection can lead to imported bancroftian filariasis (IBF) in Thailand. IBF is targeted for multiple-dose MDA with 300 mg DEC, in addition to what is recommended by the World Health Organization (WHO). The dynamic Myanmar migrants in Phang-nga, southern Thailand were sampled to test whether the responsible W. bancrofti has a genetic predisposition of benzimidazole exposure, and IBF exhibits DEC susceptibility. The long-term migrants had more access to DEC. IBF in W. bancrofti antigenemic (microfilaremic vs. amicrofilaremic) short-term migrants exhibited susceptibility to a 300-mg single-dose DEC treatment. During the course of a 3-month follow-up, antigenemia was significantly reduced, but microfilaremia was fluctuated. Surprisingly, a newly recognized Mansonella infection co-existing among W. bancrofti-affected Myanmar migrants elicited microfilaremia clearance within a month after treatment. As a result of the presence of genetically stable W. bancrofti β-tubulin (Wb tubb) gene responsible for benzimidazole susceptibility, IBF did not possess a genetic predisposition for benzimidazole exposure. Point mutations at positions Phe167Tyr and Phe200Tyr were not detected by Wb tubb locus-specific nested PCR and sequencing. This study has the potential to help guide not only the Thai/Myanmar PELF surveillance and monitoring of mass treatment impacts on W. bancrofti, but also the other endemic countries allied with the Global Program to Eliminate Lymphatic Filariasis (GPELF).

Tracing isolates of bacterial species by multilocus variable number of tandem repeat analysis (MLVA)

All bacterial genomes contain multiple loci of repetitive DNA. Repeat unit sizes and repeat sequences may vary when multiple loci are considered for different isolates of an individual microbial species. Moreover, it has been documented on many occasions that the number of repeat units per locus is a strain-defining parameter. Consequently, there is isolate-specificity in the number of repeats per locus when different strains of a given bacterial species are compared. The experimental assessment of this variability for a number of different loci has been called 'multilocus variable number of tandem repeat analysis' (MLVA). The approach can be supported or extended by locus-specific DNA sequencing for establishing mutations in the individual repeat units, which usually enhances the resolution of the approach considerably. Essentially, MLVA with or without supportive sequencing has been developed for all of the medically relevant bacterial species and can be used effectively for tracing outbreaks or other forms of bacterial dissemination. MLVA is a modern, timely and versatile bacterial typing methodology.

Simultaneous extraction of nuclear and mitochondrial DNA from human blood

A method of simultaneous isolation of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) from human blood has been proposed by improvising Lahiri's method of isolation of nuclear DNA. The approach presented here provides selectively enriched fractions and eliminates the need for two different methods or separate reagent sets for the extraction of nDNA and mtDNA. It employs an initial nuclear/ cytoplasm partitioning, followed by the similar procedural steps for the two fractions separately. It gives good quality and quantity of the nDNA as well as the mtDNA, suitable for processes like PCR amplification and sequencing and may prove to be useful for people studying population genetics and evolution using molecular markers maximizing the available resources, especially in cases where a large database needs to be generated from limited amount of blood sample. From 3 ml of blood, the yields of mtDNA salvaged from the supernatant were sufficient to set approximately 4x10(5) reactions (starting with 250 fg DNA per reactions) of mtDNA loci which otherwise would have been discarded as per original Lahiri's procedure. The quality of mtDNA from the mitochondrial fraction was suitable for all major downstream processes as confirmed by locus specific PCR amplifications and sequencing. Through this procedure, the wastage of nDNA can be avoided when mtDNA loci is studied.