PCR conditions.
Abstract
Mitochondrial DNA is an important tool for human identification and is used to differentiate between human and animal blood at the crime scene, because in extreme conditions nuclear DNA is severely destroyed while Mitochondrial DNA contains multiple copies (200–2000) per cell and resists harsh and more stable conditions. Seventy-two blood samples were collected from humans (Homo sapiens), sheep (Ovis aries), goats (Capra hircus), and cows (Bos taurus) (18 blood samples for each). All blood samples were withdrawn by a technician and 5 ml were aspirated using an aseptic technique and transferred to EDTA-Na2 tubes. They were mixed well and stored in a refrigerator. The collection took 2 weeks (May 15, 2019–May 30, 2019). All samples were collected from Al-Diwanyia city. The results of PCR testing revealed that the primer pairs were specific and non-specific products did not appear for all samples. The amplification of Homo sapiens mitochondrial DNA with primer pairs of other (Ovis aries, Capra hircus, and Bos taurus) and amplification of each with primer pairs of another genus gave negative results, and this is primary evidence for primer pair specificity. The amplicon of 16S rRNA gene of Homo sapiens was 1200 bp, Ovis aries was 1060 bp, Capra hircus was 820 bp, and Bos taurus was 1300 bp. The sequencing revealed that no cross-reactivity of designed primer pairs and the PCR assay based on the designed primer pairs will be simple, fast, sensitive, specific, and cost-effective. There is sensitivity, specificity, and accuracy in the designed species-specific primer pairs and applicability of the designed primer pairs in forensics to investigate blood spots or evidence belonging for human, sheep, goat, and cow.
Keywords
- Homo sapiens
- Ovis aries
- Capra hircus
- Bos taurus
- forensic
1. Introduction
Mitochondria possess a small spherical genome, mtDNA, which encodes for the 13 important subunits of the electron transport chain and ATP synthase together with 22 tRNAs and 2 rRNAs necessary for mitochondrial protein synthesis [1, 2]. Mitochondrial DNA presents several characteristics which have the potential to be valuable for forensic studies, especially attendant to the absence of recombination, to a large copy number, and to matrilineal inheritance. Mitochondrial DNA typing founded on sequences of the control region otherwise filled genomic sequence is used to examine a variation of forensic mtDNA profiling methods used for human proof of identity and present their use in the chief cases of human identification from non-human [3, 4, 5]. Mitochondrial markers that are used for species identification are as follows: cytb gene, cytochrome c oxidase subunit I gene, 12S and 16S rRNA segment, and control region in wildlife [6, 7, 8]. A short fragment of the 12S rDNA was employed for DNA amplification leading to species identification. The mitochondrial DNA 16S rRNA gene is an advanced genetic marker for animal genetic diversity. Polymorphism sites, nucleotide variation, and haplotype variety were determined using whole sequences of the mitochondrial DNA 16S rDNA gene [9, 10]. Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization [11]. The aim of this study is to design valuable species-specific-PCR tool to discriminate blood of humans from non-human using a species-specific primer design.
2. Methodology
2.1 Study design
The study design was experimental to design species-specific primer pairs for typing the blood samples and their assignment to human (
2.2 Blood sample collection
Seventy-two blood samples were collected from humans, sheep, goats, and cows (18 blood samples for each). All blood samples were withdrawn by technicians and 5 ml were aspirated using an aseptic technique and transferred to EDTA-Na2 tubes and mixed well and stored in a refrigerator. The collection took 2 weeks (May 15, 2019–May 30, 2019). All samples were collected from Al-Diwanya city.
2.3 Primer design
The gene selected for this study is the mitochondrial 16S rRNA gene. The NCBI data base was used to recover the sequences chosen for a primer design. The sequence ID of
2.4 Mitochondrial DNA extraction
G-spin™ Total DNA Extraction Kit (50 Preps) (REF: 17045) was used to extract mitochondrial DNA from blood of different species according to the manufacturer’s protocol instructions.
2.5 Agarose gel electrophoresis
Agarose gel was prepared by dissolving agarose powder in 1X TBE buffer. The amount of agarose can be dissolved depending upon the purpose in which agarose sheet used. About 0.7% agarose gel was used for visualization of the DNA after extraction while 1.5–2% agarose sheet was used for visualization of PCR product (amplicon). RedSafe (alternative for ethidium bromide) stock solution concentration was 10 mg/ml. Only 5 μl of RedSafe stock solution were added to 100 ml of melted agarose gel to get the final concentration of 0.5 μg/ml [13, 14].
2.6 Primer pairs preparation and PCR conditions
The primers were synthesized at Macrogen/Korea, were provided in a lyophilized form, which were re-dissolved with 300 nuclease-free water according to the institution of the manufacture company to reach to the final concentration (100 pmoles/μl). The working solution will be 10 pmoles/μl to be used directly in PCR [15, 16]. The PCR conditions were calculated using online Protocol Optimize writer software. The conditions were illustrated in Table 1.
Primer | Conditions | References | ||
---|---|---|---|---|
Homo 16S-F Homo 16S-R | 1 | 95°C | 2 min | This study |
30 | 95°C | 30 sec | ||
59.3 °C | 30 sec | |||
72°C | 130 sec | |||
1 | 72°C | 5 min | ||
Ovis 16S-F Ovis 16S-R | 1 | 95°C | 2 min | This study |
30 | 95°C | 30 sec | ||
58.3 °C | 30 sec | |||
72°C | 130 sec | |||
1 | 72°C | 5 min | ||
Capra 16S-F Capra 16S-R | 1 | 95°C | 2 min | This study |
30 | 95°C | 30 sec | ||
58.3°C | 30 sec | |||
72°C | 130 sec | |||
1 | 72°C | 5 min | ||
Bos 16S-F Bos 16S-R | 1 | 95°C | 2 min | This study |
30 | 95°C | 30 sec | ||
59.3°C | 30 sec | |||
72°C | 130 sec | |||
1 | 72°C | 5 min |
3. Result and discussion
The four sets of designed primer pairs were submitted to specificity using Primer-Blast and the results revealed that, they are specific to amplify the 16S rRNA gene of humans (
Gene | Primer sequence 5′ to 3′ | Sequence ID of isolate | Identity |
---|---|---|---|
F:GCCTGGTGATAGCTGGTTGT R:ATCATTTACGGGGGAAGGCG | MN115376.1 MN053904.1 MN125706.1 MN163828.1 MN163832.1 MN125705.1 MN163282.1 MN125704.1 MN124446.1 MK069579.1 | 100% | |
16S rRNA | F:AGGCCTAAAAGCAGCCATCA R:GCCCTTTTCTAGGGCAGGTT | KP998473.1 KP998472.1 KP998470.1 KP702285.1 MH841968.1 MH841967.1 MH841966.1 MG837554.1 MG837553.1 KU681224.1 | 100% |
F:GCCTGGTGATAGCTGGTTGT R:TCACCCCAACCAAAACTGCT | LS992662.1 LS992661.1 LS992659.1 LS992658.1 LS992656.1 LS992655.1 LS992654.1 LS992653.1 LS992652.1 LS992651.1 | 100% | |
F:CTAAGCAGCCCGAAACCAGA R:GGGCAGGGTTTTGTGTTGTC | EU177866.1 EU177865.1 EU177864.1 EU177863.1 EU177862.1 EU177861.1 EU177860.1 EU177859.1 EU177858.1 EU177856.1 | 100% |
The results of PCR testing revealed that the primer pairs were specific and non-specific products did not appear for all samples. The amplification of Homo sapiens mtDNA with primer pairs of other (Ovis aries, Capra hircus, and Bos taurus) and amplification of each with primer pairs of another genus gave negative results. This was primary evidence for primer pair specificity. The amplicon of 16S rRNA gene of Homo sapiens was 1200 bp (Figure 1A), Ovis aries was 1060 bp (Figure 1B), Capra hircus was 820 bp (Figure 1C), and Bos taurus was 1300 bp (Figure 1D). PCR amplification and sequence analysis of the mitochondrial 16S rRNA gene was utilized for differentiation/identification and subsequently evaluation of their application in solving forensic cases [21]. Mitochondrial 16S is suitable for the differentiation of 300 mammalian species. The 16S rDNA gene is a common mitochondrial gene for detection of blend mutton and pork at high sensitivity. The mitochondrial 16S rRNA genes have been used as molecular markers to identify mammals, birds, shrimps, and other species using species-specific primers that amplify the 12S rRNA or 16S rRNA gene regions from mtDNA [17, 22]. Gene loci on the mitochondrial genome have been used in species identification. These include the 12S and 16S rRNA loci. The D-loop (displacement loop) has been used less in species identification but more in intraspecies identification. Due to the greater sequence variation at this non-coding locus, it is now being used as a tool for identifying the presence of particular species within mixture of many species [23, 24].
The secondary and confirmatory assay for specificity of primer pairs used in the study was sequences of PCR products. Eight amplicons from each were sent for sequencing using the Sanger technique (Macrogen/Korea). The retrieved sequences firstly must be trimmed to remove unwanted sequences before submitting them for BLASTN. The trimming performed by Bioedit was utilized to obtain the finally processed sequences. Abbreviations of Homo sapiens sequences were used as (HIS-1 to HIS-8), Ovis aries sequences be (IOA-1–IOA-8), Capra hircus sequences (IBCH-1–IBCH-8), and Bos taurus sequences (IBT-1–IBT-8).
The identity percentage and alignment results of the amplified 16S rRNA gene of Homo sapiens, Ovis aries, Capra hircus, and Bos taurus with database were illustrated in Tables 3–6, respectively.
Isolate | Sequence ID | Expect | Identities | Gaps | Strand |
---|---|---|---|---|---|
IHS-1 | MH444415.1 | 0.0 | 98.98% | 0/885(0%) | Plus/Plus |
IHS-2 | MK069579.1 | 0.0 | 99.74% | 0/771(0%) | Plus/Plus |
IHS-3 | MK069579.1 | 0.0 | 99.57% | 2/697(0%) | Plus/Plus |
IHS-4 | MK069579.1 | 0.0 | 99.08% | 0/654(0%) | Plus/Plus |
IHS-5 | MK059695.1 | 0.0 | 99.86% | 0/701(0%) | Plus/Plus |
IHS-6 | MK069579.1 | 0.0 | 99.39% | 0/657(0%) | Plus/Plus |
IHS-7 | MK295855.1 | 0.0 | 99.50% | 0/599(0%) | Plus/Plus |
IHS-8 | MK069579.1 | 0.0 | 99.69% | 0/653(0%) | Plus/Plus |
Isolate | Sequence ID | Expect | Identities | Gaps | Strand |
---|---|---|---|---|---|
IOA-1 | MG489885.1 | 0.0 | 98.70% | 0/769(0%) | Plus/Plus |
IOA-2 | MG489885.1 | 4e-170 | 100.00% | 0/329(0%) | Plus/Plus |
IOA-3 | MG489885.1 | 0.0 | 98.98% | 0/586(0%) | Plus/Plus |
IOA-4 | MG489885.1 | 0.0 | 99.73% | 0/749(0%) | Plus/Plus |
IOA-5 | MG489885.1 | 0.0 | 99.63% | 0/542(0%) | Plus/Plus |
IOA-6 | MG489885.1 | 2e-157 | 99.36% | 0/312(0%) | Plus/Plus |
IOA-7 | MG489885.1 | 0.0 | 99.17% | 2/483(0%) | Plus/Plus |
IOA-8 | MG489885.1 | 0.0 | 99.80% | 0/489(0%) | Plus/Plus |
Isolate | Sequence ID | Expect | Identities | Gaps | Strand |
---|---|---|---|---|---|
IBCH-1 | KP271023.1 | 0.0 | 99.02% | 0/614(0%) | Plus/Plus |
IBCH-2 | KP271023.1 | 0.0 | 98.58% | 0/633(0%) | Plus/Plus |
IBCH-3 | KP271023.1 | 2e-174 | 99.74% | 0/378(0%) | Plus/Plus |
IBCH-4 | KP271023.1 | 0.0 | 100.00% | 0/729(0%) | Plus/Plus |
IBCH-5 | KP271023.1 | 0.0 | 99.83% | 0/595(0%) | Plus/Plus |
IBCH-6 | KP271023.1 | 0.0 | 100.00% | 0/480(0%) | Plus/Plus |
IBCH-7 | KP271023.1 | 0.0 | 99.84% | 0/618(0%) | Plus/Plus |
IBCH-8 | KP271023.1 | 0.0 | 99.17% | 0/481(0%) | Plus/Plus |
Isolate | Sequence ID | Expect | Identities | Gaps | Strand |
---|---|---|---|---|---|
IBT-1 | MF169214.1 | 0.0 | 99.50% | 3/601(0%) | Plus/Plus |
IBT-2 | KT184455.1 | 7e-139 | 100.00% | 0/273(0%) | Plus/Plus |
IBT-3 | KT184466.1 | 0.0 | 99.32% | 7/1177(0%) | Plus/Plus |
IBT-4 | KT184466.1 | 0.0 | 99.90% | 0/979(0%) | Plus/Plus |
IBT-5 | KT184466.1 | 0.0 | 99.90% | 0/965(0%) | Plus/Plus |
IBT-6 | KT184466.1 | 0.0 | 100.00% | 0/512(0%) | Plus/Plus |
IBT-7 | KT184466.1 | 0.0 | 99.81% | 1/1077(0%) | Plus/Plus |
IBT-8 | KT184466.1 | 0.0 | 100.00% | 0/1092(0%) | Plus/Plus |
The sequencing of the 16S rRNA has revolutionized the study and identification of human and non-human samples in forensic science. A simple method was developed using universal primers for species identification based on direct PCR sequencing using primer sets that were designed based on the conserved regions of the 16S rRNA loci detected by the comprehensive sequence comparison among 30 animals whole [25]. The mitochondrial DNA method could be a dominant tool for mammalian species identification, especially in forensic cases in which many unidentified biological samples need to be analyzed such as blood spots [25].
The 16S and 12S sequences allow identification of most species to the genus level. Faster-evolving DNA regions are required to identify closely related animal species [26]. The successfully used forensically informative nucleotide sequencing analysis of the 16S rRNA mitochondrial DNA were very valuable to identify before unknown biological specimens of human and animals [27]. The mitochondrial 12S rRNA and 16S rRNA genes, including those from fish and amphibians to mammals including human beings. Therefore, universal primers were designed to amplify sequences in the fast-evolving animal mtDNA [17]. The PCR amplifications of mitochondrial 16S rRNA followed by sequencing and analysis were demonstrated to be very efficient for identification of species origin [21]. The 12S rRNA and 16S rRNA gene sequences of animals reveal the fitting level of interspecific variation but the great level of intraspecific homogeneity [7]. The results showed no cross-reactivity of designed primer pairs and the PCR assay based on the designed primer pairs will be simple, fast, sensitive, specific, and cost-effective.
3.1 Strength of PCR-dependent 16S mtDNA gene
The extraction of mtDNA and amplification of mtDNA genes seem to be accessible, very easy, and cheap. Additionally when the sequences of amplified genes analyzed the results were very clear and no confusion with other genes. The accuracy is very high due to no cross-amplification between species-specific primers that were observed. The sensitivity is also high due to that the mixed blood at a very small amount (10%) can be detected. Their strength over another technique like real-time PCR was, post real time, cannot perform the sequencing when needed for verification.
3.2 Weakness of PCR-dependent 16S mtDNA gene
The only weakness in the technique is requirement for more time, preparation, possibility of contamination, and more machine when compared with real-time PCR.
4. Conclusions
There is sensitivity, specificity, and accuracy of the designed species-specific primer pairs and applicability of the designed primer pairs in forensics to investigate blood spots or evidence belonging for human, sheep, goat, and cow.
5. Registration of sequences in GenBank
All the 32 sequences of the 16S rRNA gene were submitted to GenBank for registration. After checking and revision, the following accession numbers were donated:
16S rRNA homo sapiens (Human): MN192057, MN192058, MN192059, MN192060, MN192061, MN192062, MN192063, MN192064(Appendix 4–69–4–76).
16S rRNA Ovis
16S rRNA Capra hircus (Goat): MN173285, MN173286, MN173287, MN173288, MN173289, MN173290, MN173291, MN173292(Appendix 4–85–4–92).
16S rRNA Bos taurus (Cow): MN197611, MN197612, MN197613, MN197614, MN197615, MN197616, MN197617, MN197618 (Appendix 4–93–4–100).
References
- 1.
D’Souza AR, Minczuk M. Mitochondrial transcription and translation: Overview. Essays in Biochemistry. 2018; 62 (3):309-320 - 2.
Herbers E, Kekäläinen NJ, Hangas A, Pohjoismäki JL, Goffart S. Tissue specific differences in mitochondrial DNA maintenance and expression. Mitochondrion. 2019; 44 :85-92 - 3.
Indo HP, Suenaga S, Tomita K, Ito H, Matsui H, Majima HJ. Analysis of oxidative stress marker, mitochondrial DNA copy numbers and mitochondrial DNA oxidation among 135 persons who live in Amami islands, a high centenarian population district in Kagoshima. Free Radical Biology and Medicine. 2018; 120 :S134 - 4.
Amorim A, Fernandes T, Taveira N. Mitochondrial DNA in human identification: A review. PeerJ Preprints. 2019; 7 :e27500v1 - 5.
Lee WC, Lin CS, Ko FC, Cheng W, Lee MH, Wei YH. Low mitochondrial DNA copy number of resected cecum appendix correlates with high severity of acute appendicitis. Journal of the Formosan Medical Association. 2019; 118 (1):406-413 - 6.
Mitra I, Roy S, Haque I. Application of molecular markers in wildlife DNA forensic investigations. Journal of Forensic Science and Medicine. 2018; 4 (3):156 - 7.
Mahmoodi M, Afshari KP, Seyedabadi HR, Aboozari M. Sequence analysis of 12S rRNA and 16S rRNA mitochondrial genes in Iranian Afshari sheep. Banat’s Journal of Biotechnology. 2018; 9 (18):5-11 - 8.
Andrejevic M, Markovic MK, Bursac B, Mihajlovic M, Tanasic V, Kecmanovic M, et al. Identification of a broad spectrum of mammalian and avian species using the short fragment of the mitochondrially encoded cytochrome b gene. Forensic Science, Medicine and Pathology. 2019; 15 (2):169-177 - 9.
Horreo JL, Fitze PS, Jiménez-Valverde A, Noriega JA, Pelaez ML. Amplification of 16S rDNA reveals important fish mislabeling in Madrid restaurants. Food Control. 2019; 96 :146-150 - 10.
Yan L, She Y, Elzo MA, Zhang C, Fang X, Chen H. Exploring genetic diversity and phylogenic relationships of Chinese cattle using gene mtDNA 16S rRNA. Archives Animal Breeding. 2019; 62 (1):325-333 - 11.
Lavrov DV, Pett W. Animal mitochondrial DNA as we do not know it: mt-genome organization and evolution in nonbilaterian lineages. Genome Biology and Evolution. 2016; 8 (9):2896-2913 - 12.
Untergrasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—New capabilities and interfaces. Nucleic Acids Research. 2012 August;40(15):e115 - 13.
Ahmed NS, Hadi YA, Dhefer IH. Polymorphism study of TCF7L2 gene and related to some biochemical parameters in DM2 females Iraqi patients. Research Journal of Science and Technology. 2019; 11 (1):01-08 - 14.
Jarrar YB, Ghishan M. The Nudix hydrolase 15 (NUDT15) gene variants among Jordanian Arab population. Asian Pacific Journal of Cancer Prevention: APJCP. 2019; 20 (3):801-808 - 15.
Cseke LJ, Kirakosyan A, Kaufman PB, Westfall MV. Handbook of Molecular and Cellular Methods in Biology and Medicine. CRC Press; 2016 - 16.
Heuvel JP. PCR Protocols in Molecular Toxicology. CRC Press; 2019 - 17.
ang L, Tan Z, Wang D, Xue L, Guan MX, Huang T, et al. Species identification through mitochondrial rRNA genetic analysis. Scientific Reports. 2014; 4 :4089 - 18.
Luo S, Valencia CA, Zhang J, Lee NC, Slone J, Gui B, et al. Biparental inheritance of mitochondrial DNA in humans. Proceedings of the National Academy of Sciences. 2018; 115 (51):13039-13044 - 19.
Barshad G, Marom S, Cohen T, Mishmar D. Mitochondrial DNA transcription and its regulation: An evolutionary perspective. Trends in Genetics. 2018; 34 (9):682-692 - 20.
Sharifi RS, Sofla SS, Seyedabadi HR. Genetic diversity and molecular phylogeny of iranian goats based on cytochrome oxidase I (COXI) gene sequences. Jurnal Veteriner. 2018; 18 (4):565-570 - 21.
Mane BG, Mendiratta SK, Tiwari AK, Narayan R. Sequence analysis of mitochondrial 16S rRNA gene to identify meat species. Journal of Applied Animal Research. 2013; 41 (1):77-81 - 22.
Xu J, Zhao W, Zhu M, Wen Y, Xie T, He X, et al. Molecular identification of adulteration in mutton based on mitochondrial 16S rRNA gene. Mitochondrial DNA Part A. 2016; 27 (1):628-632 - 23.
Mitani T, Akane A, Tokiyasu T, Yoshimura S, Okii Y, Yoshida M. Identification of animal species using the partial sequences in the mitochondrial 16S rRNA gene. Legal Medicine. 2009; 11 :S449-S450 - 24.
Linacre A, Tobe SS. An overview to the investigative approach to species testing in wildlife forensic science. Investigative Genetics. 2011; 2 (1):2 - 25.
Kitano T, Umetsu K, Tian W, Osawa M. Two universal primer sets for species identification among vertebrates. International Journal of Legal Medicine. 2007; 121 (5):423-427 - 26.
Cawthorn DM, Steinman HA, Witthuhn RC. Evaluation of the 16S and 12S rRNA genes as universal markers for the identification of commercial fish species in South Africa. Gene. 20121; 491 (1):40-48 - 27.
Guha S, Goyal SP, Kashyap VK. Genomic variation in the mitochondrially encoded cytochrome b (MT-CYB) and 16S rRNA (MT-RNR2) genes: Characterization of eight endangered Pecoran species. Animal Genetics. 2006; 37 (3):262-265