Mitochondria analysis

0(0명)
문서 역사

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.

KSJ mitochondria genome.

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