Sequencing quality comparison based on MtDNA
The rCRS genome has accepted as a refercence genome for mitochondria.
We compaired the rCRS genome with KoRef mtDNA.
As the Table1, the two genomes has showed very similar in sequence composition.
Sequences of genes were aligned very clear (not appear stop codons). (figure 1)
Table 1. Comparison rCRS with KoRef mtDNA in sequence composition.
|
genome |
All |
A |
T |
G |
C |
AT% |
ATskew |
GCskew |
ATratio |
GCratio |
|
CRS |
16569 |
5124 |
4094 |
2169 |
5181 |
56 |
0.112 |
-0.41 |
1.252 |
0.419 |
|
Korean_MT |
16571 |
5113 |
4086 |
2180 |
5192 |
55.6 |
0.112 |
-0.405 |
1.251 |
0.42 |
Variation aspect of mitochondrial genome
13 mitochondrial protein coding genes, 12 rRNA genes, and 22 tRNA genes were identified by aligning with the rCRS mitochondrial genome sequences using Clustal W.

Figure 1. Annotation and visualization of KoRef mitochondria genome. Variations and genes of the genome were identified through sequence alignment against rCRS.
Table 2. Genes involved in KoRef mitochondria genome.
|
Gene Name |
start |
stop |
strand |
|
HVS2 |
1 |
57 |
|
|
tRNA-Phe |
579 |
649 |
plus |
|
12s_rRNA |
650 |
1603 |
plus |
|
tRNA-Val |
1604 |
1672 |
plus |
|
16s_rRNA |
1673 |
3230 |
plus |
|
tRNA-Leu |
3231 |
3305 |
plus |
|
ND1 |
3308 |
4264 |
plus |
|
tRNA-Ile |
4264 |
4332 |
plus |
|
tRNA-Gln |
4330 |
4401 |
minus |
|
tRNA-Met |
4403 |
4470 |
plus |
|
ND2 |
4471 |
5512 |
plus |
|
tRNA-Trp |
5513 |
5580 |
plus |
|
tRNA-Ala |
5588 |
5656 |
minus |
|
tRNA-Asn |
5658 |
5730 |
minus |
|
tRNA-Cys |
5762 |
5827 |
minus |
|
tRNA-Tyr |
5827 |
5892 |
minus |
|
COX1 |
5905 |
7446 |
plus |
|
tRNA-Ser |
7446 |
7517 |
minus |
|
tRNA-Asp |
7519 |
7586 |
plus |
|
COX2 |
7587 |
8270 |
plus |
|
tRNA-Lys |
8296 |
8365 |
plus |
|
ATP8 |
8367 |
8573 |
plus |
|
ATP6 |
8528 |
9208 |
plus |
|
COX3 |
9208 |
9991 |
plus |
|
tRNA-Gly |
9992 |
10059 |
plus |
|
ND3 |
10060 |
10405 |
plus |
|
tRNA-Arg |
10406 |
10470 |
plus |
|
ND4L |
10471 |
10767 |
plus |
|
ND4 |
10761 |
12138 |
plus |
|
tRNA-His |
12139 |
12207 |
plus |
|
tRNA-Ser2 |
12208 |
12266 |
plus |
|
tRNA-Leu2 |
12267 |
12337 |
plus |
|
ND5 |
12338 |
14149 |
plus |
|
ND6 |
14150 |
14674 |
minus |
|
tRNA-Glu |
14675 |
14743 |
minus |
|
CYTB |
14748 |
15888 |
plus |
|
tRNA-Thr |
15889 |
15954 |
plus |
|
tRNA-Pro |
15956 |
16024 |
minus |
|
HVS1 |
16025 |
16571 |
|
Table 3. Nucleotide substitutions in KoRef mtDNA against rCRS genome.
| # | Position | rCRS | KSJ | Class | Gene | Type | rCRS_aa | KSJ_aa | rCRS_nt | KSJ_nt |
| 1 | 73 | A | G | single | - | control_region | ||||
| 2 | 150 | C | T | single | - | control_region | ||||
| 3 | 195 | T | C | single | - | control_region | ||||
| 4 | 263 | A | G | single | - | control_region | ||||
| 5 | 310 | - | T | insertion | - | control_region | ||||
| 6 | 310 | T | C | single | - | control_region | ||||
| 7 | 311 | - | C | insertion | - | control_region | ||||
| 8 | 408 | T | A | single | - | control_region | ||||
| 9 | 750 | A | G | single | 12s_rRNA | rRNA | ||||
| 10 | 1438 | A | G | single | 12s_rRNA | rRNA | ||||
| 11 | 2352 | T | C | single | 16s_rRNA | rRNA | ||||
| 12 | 2483 | T | C | single | 16s_rRNA | rRNA | ||||
| 13 | 2706 | A | G | single | 16s_rRNA | rRNA | ||||
| 14 | 3107 | X | T | single | 16s_rRNA | rRNA | ||||
| 15 | 3109 | T | - | deletion | 16s_rRNA | rRNA | ||||
| 16 | 4769 | A | G | single | ND2 | synonymous | Met | Met | ATA | ATG |
| 17 | 5580 | T | C | single | tRNA-Trp | tRNA | ||||
| 18 | 7028 | C | T | single | COX1 | synonymous | Ala | Ala | GCC | GCT |
| 19 | 8701 | A | G | single | ATP6 | non-synonymous | Thr | Ala | ACC | GCC |
| 20 | 8860 | A | G | single | ATP6 | non-synonymous | Thr | Ala | ACA | GCA |
| 21 | 9377 | A | G | single | COX3 | synonymous | Trp | Trp | TGA | TGG |
| 22 | 9540 | T | C | single | COX3 | synonymous | Leu | Leu | TTA | CTA |
| 23 | 10398 | A | G | single | ND3 | non-synonymous | Thr | Ala | ACC | GCC |
| 24 | 10819 | A | G | single | ND4 | synonymous | Lys | Lys | AAA | AAG |
| 25 | 10873 | T | C | single | ND4 | synonymous | Pro | Pro | CCT | CCC |
| 26 | 11017 | T | C | single | ND4 | synonymous | Ser | Ser | AGT | AGC |
| 27 | 11719 | G | A | single | ND4 | synonymous | Gly | Gly | GGG | GGA |
| 28 | 11722 | T | C | single | ND4 | synonymous | Leu | Leu | CTT | CTC |
| 29 | 12705 | C | T | single | ND5 | synonymous | Ile | Ile | ATC | ATT |
| 30 | 12850 | A | G | single | ND5 | non-synonymous | Ile | Val | ATC | GTC |
| 31 | 14212 | T | C | single | ND6 | synonymous | Val | Val | GTA | GTG |
| 32 | 14580 | A | G | single | ND6 | synonymous | Leu | Leu | TTG | CTG |
| 33 | 14766 | C | T | single | CYTB | non-synonymous | Thr | Ile | ACT | ATT |
| 34 | 14905 | G | A | single | CYTB | synonymous | Met | Met | ATG | ATA |
| 35 | 15301 | G | A | single | CYTB | synonymous | Leu | Leu | TTG | TTA |
| 36 | 15326 | A | G | single | CYTB | non-synonymous | Thr | Ala | ACA | GCA |
| 37 | 15932 | T | C | single | tRNA-Thr | tRNA | ||||
| 38 | 16172 | T | C | single | HVS1 | control_region | ||||
| 39 | 16183 | A | C | single | HVS1 | control_region | ||||
| 40 | 16189 | T | C | single | HVS1 | control_region | ||||
| 41 | 16193 | - | C | insertion | HVS1 | control_region | ||||
| 42 | 16223 | C | T | single | HVS1 | control_region | ||||
| 43 | 16320 | C | T | single | HVS1 | control_region | ||||
| 44 | 16519 | T | C | single | HVS1 | control_region |
MtDNA ethno-genographic lineage
Fig XX. mtDNA ethno-genogeographic lineage
Korean mitochondrial DNA had on XX SNPs. According to 22 haplogroups defined by sequence differences from the rCRS, haplogroup D4, prevalent in Asia, mapping into our mtDNA.
Mapping disease information
We have been able to detect the expected phenotype based on the literature given the according to published literature (ref 1)
Position: A1438G
Disease: Diabetes Mellitus
Gene: 12S rRNA
Reference: Tawata, M., Ohtaka, M., Iwase, E., Ikegishi, Y., Aida, K. and Onaya, T. (1998). "New mitochondrial DNA homoplasmic mutations associated with Japanese patients with type 2 diabetes." Diabetes 47(2):276-277.
Pubmed ID: 9519725
Heteroplasmy
possibility
- paternal inheritance (ref 2)
- by aging, MtDNA variation in one indivisual is increasing.
Observed heteroplasmy
10 =< MHP < 50% : 6
5 <= MHP < 10 : 16
1 <= MHP < 5 : 1188
MHP: Minor-heteroplasmy
Table 4. Heteroplasmic mtDNA changes in KoRef mtDNA.
| Position | Nucleotide substitutions | Gene location: amino acid change | %heteroplasmy |
References
1. Tawata, M., Ohtaka, M., Iwase, E., Ikegishi, Y., Aida, K. and Onaya, T. (1998). "New mitochondrial DNA homoplasmic mutations associated with Japanese patients with type 2 diabetes." Diabetes 47(2):276-277. 9519725
2. Gustafson A. W., Heckerling P. S., Vissing J., Schwartz M. (2002). Paternal Inheritance of Mitochondrial DNA. N Engl J Med 347: 2081-2082
댓글 0