International Journal of Research and Innovation in Social Science

Submission Deadline-29th November 2024
November 2024 Issue : Publication Fee: 30$ USD Submit Now
Submission Deadline-05th December 2024
Special Issue on Economics, Management, Sociology, Communication, Psychology: Publication Fee: 30$ USD Submit Now
Submission Deadline-20th November 2024
Special Issue on Education, Public Health: Publication Fee: 30$ USD Submit Now

Comparison of Mitochondrial Genome Arrangement of Labeorohita with Some Selected Vertebrae Species.

  • Khandakar Arshadul Bari
  • Hawa Jahan
  • Reza Md. Shahjahan
  • Rowshan Ara Begum
  • 467-485
  • Apr 30, 2024
  • Science & Technology

Comparison of Mitochondrial Genome Arrangement of Labeorohita with Some Selected Vertebrae Species.

1Khandakar Arshadul Bari, 2Hawa Jahan, 2Reza Md. Shahjahan,*2Rowshan Ara Begum

1Pirganj Govt. College, Thakurgaon 

2University of Dhaka.

*Corresponding Author

DOI: https://dx.doi.org/10.47772/IJRISS.2024.804036

Received: 26 February 2024; Revised: 25 March 2024; Accepted: 30 March 2024; Published: 30 April 2024

ABSTRACT

The mitochondrial genome of Labeorohita serves as a microcosm for exploring genetic diversity and evolutionary relationships within vertebrates. Through comprehensive analysis of the mitochondrial genome, encompassing 37 genes including protein-coding genes (PCGs), transfer RNA (tRNA), and ribosomal RNA (rRNA), we elucidated its organization and compared it with other vertebrate taxa.

Our investigation revealed a compact mitochondrial genome of approximately 16 kb, comprising 13 PCGs, 22 tRNA genes, and 2 rRNA genes, with a distinctive dispersed loop region (D-loop). Notably, the arrangement of PCGs formed three distinct clusters, shared across Labeo species, suggesting evolutionary conservation within the genus. Comparative analyses with other osteichthyes species corroborated this conservation, underscoring L. rohita‘s representative status within the class.

Furthermore, comparisons with chondrichthyes, lungfish, amphibians, reptiles, birds, and mammals unveiled evolutionary trends in mitochondrial genome organization. Noteworthy findings include the unique replication origin in lungfish and amphibians, and transposition events in birds. Phylogenetic reconstructions based on mitochondrial gene sequences affirmed the close relationship between lungfish and amphibians, supporting lungfish as the closest relatives of class Amphibia.

In summary, our study provides comprehensive insights into the mitochondrial genome organization of L. rohita and its implications for vertebrate evolution. The findings contribute to a deeper understanding of genetic diversity and evolutionary relationships within vertebrates, with potential applications in taxonomy, conservation, and evolutionary biology.

Keywords: Labeorohita, mitochondrial genome, genetic diversity, evolutionary relationships, vertebrates, phylogenetic analysis, protein-coding genes, transfer RNA, ribosomal RNA, conservation.

INTRODUCTION

The genus Labeo contains 103 species, with four species found in Bangladesh (Rahman, 2005). Labeo species constitute an important group of fish with intense diversity and potential for commercial aquaculture in many Southeast Asian nations, including the Indian subcontinent.

There are only two types of DNA in animals: nuclear DNA and mitochondrial DNA (mtDNA).

Mitochondria are cellular organelles that have the function of oxidative phosphorylation and the formation of ATP. In humans, mtDNA is a double-stranded, circular, covalently closed molecule of 16.5 kb. MtDNA is inherited as a haploid from the mother, and heteroplasmy has been rarely found. From a population perspective, it could be considered as a system of small, sexually isolated demes or clonal lineages, with an evolutionary rate 5 to 10 times faster than the nuclear genome (José A. Castro, Antònia Picornell, &Misericòrdia Ramon, 1998).

Complete mitochondrial DNA analysis has proven useful in clarifying relationships among closely related species (Mayer, 1993). The 16S, 12S, Co1, and Cytb genes of the mitochondrial genome have been used as useful molecular markers for species-specific identification (Sperling & March, 1994; Girish et al., 2005).

Animal mitochondrial DNA showing position of all the genes.

Figure 1:  Animal mitochondrial DNA showing position of all the genes. Figure collected from Wikipedia (Mitochondrial DNA; http://commons.wikimedia.org).

The animal mitochondrial DNA (mtDNA) is a closed circular molecule (see Figure 1). Only the cnidarian classes Cubozoa, Scyphozoa, and Hydrozoa have been found to have linear mtDNA chromosomes (Bridge et al., 1992). In general, the vertebrate mt genome is highly compact, around 16 kb in length, and contains only 37 genes: 13 protein-coding genes, two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs), along with one long noncoding region known as the control region containing the signals for regulation and initiation of mtDNA replication and transcription (e.g., Wolstenholme, 1992). Mt gene arrangement tends to be conserved among vertebrates, with all 37 genes and the control region organized in relatively the same order from teleost fish to eutherian mammals (Boore, 1994a).

Complete mitochondrial DNA (mtDNA) sequences have been determined for 177 vertebrate species and 75 invertebrate species, with partial sequences for hundreds of others (Boore et al., 1999). All are circular DNA molecules containing 37 genes: 13 for proteins (COI-III, ND1-6, ND4L, Cytb, A6, A8); two for rRNAs (srRNA and lrRNA); and 22 for tRNAs (designated by the one-letter amino acid code, with the two S and two L tRNAs differentiated by the codons recognized [AGN/UCN and CUN/UUR, respectively]). The genes are arranged very compactly, with no introns and few intergenic nucleotides.

Studies on the genome rearrangement of Labeorohita and its comparison with other species have not been conducted yet. Genome arrangement of Labeorohita has not been compared with other chondrichthyes and osteichthyes species. Furthermore, lower vertebrate species have not been compared with genome arrangement with that of Labeorohita. Similarly, higher vertebrate species have not been compared with the genome arrangement of Labeorohita. Despite the extensive work on the complete mitochondrial genome or partial sequencing of mtDNA and identification of species based on 16S, CO1, and Cyt b genes, there has been no research on rearrangement. Therefore, I attempted to study this topic.

METHODOLOGY

Sample Collection and DNA Extraction:

  1. Specimens of Labeorohita and other vertebrate species were collected from diverse geographical locations.
  2. Genomic DNA was extracted from tissue samples using the CTAB method.

PCR Amplification and Sequencing:

  1. A short region of the mitochondrial genome, such as the 16S rRNA gene, was amplified using polymerase chain reaction (PCR) with species-specific primers.
  2. PCR products were purified and sequenced using Sanger sequencing technology.

Data Retrieval and Assembly:

  1. Mitochondrial DNA sequences of 17 species, namely L. rohita, L. bata, L. calbasu, Channa marulius, Sicamugilcascasia, Tenualosailisha, Heteropneustesfossilis, Macrognathus aculeatus, Mastacembelusarmatus, Neoceratodusforsteri, Chimaera monstrosa, Dasyatiszugei, Scoliodonlaticaudus, Bufo melanostictus, Gecko gekko, Acridotheres tristis, and Cavia porcellus, were collected from the GenBank Database (NCBI).
  2. The genes of the collected sequences were manually arranged according to their position in the collected sequence using Microsoft Office. To construct a phylogenetic tree from the nucleotide sequences of the 16S and COI genes of the 17 species mentioned above, different software programs were used, including Notepad, Serial Cloner, and SeaView.
  3. Sequences were assembled and annotated using bioinformatics tools to identify protein-coding genes, tRNA, rRNA, and the D-loop region.

Genome Arrangement Analysis:

  1. The arrangement of protein-coding genes, tRNA, rRNA, and the D-loop region in the mitochondrial genome of rohita was analyzed and compared with other vertebrate species.
  2. Manual arrangement was performed using software like Microsoft Office or specialized bioinformatics tools.

Phylogenetic Analysis:

  1. Nucleotide sequences of conserved mitochondrial genes (e.g., CoI, 16S rRNA, cytochrome b) from rohita and representative vertebrate species were aligned.
  2. Phylogenetic trees were constructed using maximum likelihood (ML), neighbor-joining (NJ), or maximum parsimony (MP) methods.
  3. Bootstrap analysis was conducted to assess the robustness of tree topologies.

Comparative Genomics:

  1. Comparative analyses were performed to identify conserved regions, gene clusters, and unique features in the mitochondrial genomes of different vertebrate taxa.
  2. Differences in gene arrangements, sequence lengths, and structural elements like replication origins were examined.

Data Interpretation and Conclusion:

  1. Results from genome arrangement analysis, phylogenetic reconstruction, and comparative genomics were interpreted to infer evolutionary relationships and genetic diversity.
  2. Conclusions regarding the representative nature of rohita mitochondrial genome, evolutionary trends, and relationships among vertebrates were drawn.

Validation and Review:

  1. Results were validated through rigorous statistical analyses and peer review processes.
  2. Feedback from experts in the field was incorporated to refine interpretations and ensure the robustness of conclusions.

RESULT & DISCUSSION

Comparison of mtDNA arrangement in different vertebrate

Figure 2: Schematic diagram of the mitochondrial genome of Labeorohita(AP 011201).Different color showing the position of 13 protein coding, 22 tRNA and 2 rRNA genes with D-loop region.

Mitochondrial Genome Arrangement:

From the NCBI, it was revealed that the complete mitochondrial genome of Labeorohita is 16 kb long. It consists of 37 genes with one dispersed loop region or D-loop region. Among them, 13 protein-coding genes, 22 tRNAs, and 2 rRNA genes are found.

1. Protein coding genes:

There are 13 protein-coding genes in the mitochondrial genome of Labeorohita, namely NADH1, NADH2, Cytochrome Oxidase I, Cytochrome Oxidase II, ATPase 8, ATPase 6, Cytochrome Oxidase III, NADH 3, NADH 4L, NADH 4, NADH 5, NADH 6, Cytochrome b, and the D-loop or control region. These 13 protein-coding genes are interrupted by 22 tRNAs except for three clusters, namely ATPase 6 – ATPase 8 – COIII, ND4L – ND4, ND5 – ND6. The arrangement of these three clusters of protein-coding genes is commonly observed in vertebrate mitochondrial genomes.

Comparing the mitochondrial genome of three Labeo species, namely Labeorohita, L. calbasu, and L. bata, it was observed that these three clusters of protein-coding genes are the same in all three species of Labeo. However, differences were observed only in gene length, which varied from 1 to 10 bp. From this observation, it was concluded that the genome arrangement of Labeorohita is a representative or typical genome arrangement for Labeo species in the case of protein-coding genes.

Typical Labeo species were compared with the genome arrangement of Tenualosailisha, Channa marulius, and Heteropneustesfossilis, and with Mastacembelusarmatus, Macrognathus aculeatus, and Sicamugilcascasia. The same result was found as before, that is, all three clusters are similar in these species. From those comparisons, it was eventually concluded that the mitochondrial genome arrangement of L. rohita can represent the class Osteichthyes for protein-coding genes. Thus, the arrangement of the mitochondrial genome of L. rohita can be described as the typical mitochondrial genome arrangement for class Osteichthyes.

Labeorohita Labeocalbasu Labeo bata Typical Labeosp
Phe, 69 Phe, 69 Phe, 69 Phe,
12S,954 12S, 954 12S,954 12S
Val,71 Val, 71 Val,71 Val
16S,1686 16S, 1689 16S,1686 16S
Leu,75 Leu, 75 Leu,75 Leu
ND1,974 ND1, 974 ND1,974 ND1
Ile,71 Ile, 71 Ile,71 Ile
Gln,70 Gln, 70 Gln,70 Gln
Met,68 Met, 68 Met, 68 Met
ND2,1044 ND2, 1044 ND2,1044 ND2
Trp,70 Trp, 70 Trp, 70 Trp,
Ala,68 Ala, 68 Ala, 68 Ala
Asn,72 Asn, 72 Asn, 72 Asn
 Cys,66 Cys, 66 Cys, 66 Cys
Tyr,70 Tyr, 70 Tyr, 70 Tyr
COI,1550 COI,  1550 COI,1550 COI
Ser,70 Ser, 70 Ser, 70 Ser
Asp,71 Asp, 71 Asp,71 Asp
COII,680 COII,  690 COII,680 COII
Lys,75 Lys, 73 Lys, 75 Lys
ATP8,164 ATP8, 164 ATP8, 164 ATP8
ATP6, 682 ATP6,  682 ATP6, 682 ATP6
COIII, 784 COIII, 785 COIII, 784 COIII
Gly,71 Gly, 71 Gly,71 Gly
ND3, 348 ND3, 348 ND3, 348 ND3
Arg, 69 Arg, 69 Arg, 69 Arg,
ND4L,296 ND4L, 296 ND4L, 296 ND4L
ND4, 1380 ND4, 1380 ND4, 1380 ND4
His, 68 His, 68 His, 68 His
Ser, 68 Ser, 68 Ser, 68 Ser
Leu, 72 Leu, 72 Leu, 72 Leu
ND5, 1823 ND5, 1823 ND5, 1823 ND5
ND6, 521 ND6, 521 ND6, 521 ND6
Glu, 68 Glu, 68 Glu, 68 Glu
Cyt b, 1140 Cyt b, 1140 Cyt b, 1140 Cyt b
Thr,71 Thr, 71 Thr,71 Thr
Pro,69 Pro, 69 Pro, 69 Pro
D-loop, 940 D-loop, 937 D-loop, 940 D-loop

Figure 3: Mitochondrial genome arrangement of three species of Labeo namely Labeorohita(AB 011201), L. calbasu (AP012143) and  L. bata (AP011198). Number showing the total length of particular genes in the boxes. Protein-coding genes are represented bold. tRNA genes are represented as three letter code.

  1. When the typical mitochondrial genome arrangement of fish is compared with chondrichthyes such as Chimaera monstrosa, Dasyatiszugei, and Scoliodonlaticaudus (see Figure 6), it was found that the arrangement of protein-coding genes is the same with the three clusters. So, it can be concluded that in the case of protein-coding genes, the genome arrangement of L. rohita is typical for chondrichthyes.
  2. When the typical mitochondrial genome arrangement of fish is compared with the lungfish Neoceratodusforsteri (see Figure 6), it was found that the arrangement of protein-coding genes is the same, and the three clusters are also the same. So, it can be concluded that in the case of protein-coding genes, the genome arrangement of L. rohita is typical for Neoceratodusforsteri.
  3. When the typical mitochondrial genome arrangement of fish is compared with amphibian Bufo melanostictus (see Figure 7), it was revealed that the protein-coding genes are the same, but the positions are different, and the number of base pairs also varies from that of the fish mitochondrial genome. Additionally, it was observed that in fish, between the clusters ND4L-ND4 and ND5-ND6, there are three tRNAs present (tRNAH, tRNAS, tRNAL2), but in the case of Amphibia (Bufo), only tRNAH and tRNAS are present, and tRNAL2 is absent. So, it can be concluded that amphibians are distant from the fish mitochondrial genome.
  4. When the typical mitochondrial genome arrangement of fish is compared with reptiles (G. gekko), birds (A. tristis), and mammals (C. porcellus), it was revealed that the protein-coding genes are the same in all cases. Additionally, the three clusters ATPase 6 – ATPase 8 – COIII, ND4L – ND4, ND5 – ND6 are similar to that of reptiles and mammals genome arrangement, but in the case of birds, there is an ND5-Cyt b cluster present instead of the ND5 – ND6 cluster. This is due to the transposition of the ND6 and tRNAGlu gene.
  5. The typical mitochondrial genome arrangement of fish was also compared with Drosophila melanogaster, and it was revealed that tRNAH, tRNAT, and tRNAP were present between the ND5 and ND6 genes (see Figure 9). The same result was found in the sessile barnacle Tetraclita japonica (Rowshan et al., 2004). Therefore, it can be concluded that the mitochondrial genome arrangement between fish and invertebrates is significantly different when protein-coding genes are considered.
Typical Labeosp Tenualosa Channa Heteropneustes
Phe, Phe, 69 Phe, 70 Phe, 69
12S 12S, 953 12S, 944 12S, 954
Val Val, 71 Val, 71 Val, 71
16S 16S, 1679 16S, 1685 16S, 1663
Leu Leu, 75 Leu, 75 Leu, 74
ND1 ND1, 974 ND1, 974 ND1, 974
Ile Ile, 71 Ile, 69 Ile, 71
Gln Gln, 70 Gln, 70 Gln,70
Met Met, 68 Met, 69 Met,69
ND2 ND2, 1049 ND2, 1045 ND2, 1044
Trp, Trp, 71 Trp, 70 Trp, 70
Ala Ala, 68 Ala, 68 Ala, 68
Asn Asn, 73 Asn, 72 Asn, 72
Cys Cys, 65 Cys, 66 Cys, 65
Tyr Tyr, 69 Tyr, 69 Tyr, 69
COI COI, 1550 COI, 1550 COI, 1550
Ser Ser, 70 Ser, 70 Ser,70
Asp Asp, 68 Asp, 71 Asp, 70
COII COII, 690 COII, 690 COII, 690
Lys Lys, 73 Lys, 74 Lys, 73
ATP8 ATP8, 167 ATP8, 167 ATP8, 167
ATP6 ATP6, 683 ATP6,  682 ATP6, 682
COIII COIII, 785 COIII, 784 COIII, 783
Gly Gly, 70 Gly, 68 Gly, 72
ND3 ND3, 350 ND3, 348 ND3, 348
Arg, Arg, 68 Arg,69 Arg, 68
ND4L ND4L, 296 ND4L,296 ND4L, 296
ND4 ND4, 1380 ND4, 1380 ND4, 1380
His His, 68 His, 69 His,69
Ser Ser,  66 Ser, 69 Ser, 66
Leu Leu, 70 Leu, 72 Leu, 72
ND5 ND5, 1847 ND5, 1838 ND5, 1826
ND6 ND6, 521 ND6, 521 ND6, 518
Glu Glu, 68 Glu, 67 Glu, 68
Cyt b Cyt b, 1140 Cyt b, 1140 Cyt b, 1137
Thr Thr, 71 Thr,71 Thr, 72
Pro Pro, 69 Pro,69 Pro, 69
D-loop D-loop,  1168 D-loop, 911 D-loop, 858

Figure 4: Compare the mitochondrial genome arrangement of Labeo with those of Tenualosailisha[Clupeiformes, (AP 011610)], Channa marulius [Channiformes(AB 968638)], Heteropneustesfossilis [Siluriformes (AP 012013)].

2. tRNA:

There are 22 tRNA genes observed in L. rohita, consistent with the observation of Wolfe (1993). The tRNA genes of L. rohita range in size from 66 to 76 bp, and they form 5 clusters, namely WANCY, IQM, SD, HSL2, and TP (see Figure 3).

  1. The arrangement of the 22 tRNA genes of three Labeo species, L. rohita, L. calbasu, and L. bata, was studied, and the aforementioned 5 clusters are the same in all three species. However, the lengths differ from each other by only 1-2 bp (see Figure 3). From this observation, we can conclude that the genome arrangement of Labeorohita is a representative or typical genome arrangement for Labeo species for tRNA genes.
  2. Subsequently, the typical Labeo species’ tRNA gene arrangement was compared with Tenualosailisha, Channa marulius, and Heteropneustesfossilis (see Figure 4), and with Mastacembelusarmatus, Macrognathus aculeatus, and Sicamugilcascasia (see Figure 5). The same result was found as before; all 5 clusters are similar in these species. From these comparisons, it was concluded that the mitochondrial genome arrangement of L. rohita can represent the class Osteichthyes for tRNA gene arrangement. Thus, the arrangement of the mitochondrial genome of L. rohita can be described as the typical mitochondrial genome arrangement for class Osteichthyes or the typical mitochondrial genome arrangement for fish.
Typical Labeosp Mastacembelus Macrognathus Sicamugil
Phe, Phe, 69 Phe, 69 Phe, 69
12S 12S, 948 12S, 944 12S, 950
Val Val, 71 Val, 72 Val, 70
16S 16S, 1670 16S, 1669 16S, 1696
Leu Leu, 73 Leu, 74 Leu, 73
ND1 ND1, 971 ND1, 971 ND1, 974
Ile Ile, 70 Ile, 70 Ile, 69
Gln Gln, 70 Gln, 70 Gln, 70
Met Met, 68 Met, 68 Met, 69
ND2 ND2, 1045 ND2, 1045 ND2, 1044
Trp, Trp, 71 Trp, 71 Trp, 71
Ala Ala, 68 Ala, 68 Ala, 68
Asn Asn, 72 Asn, 72 Asn, 72
Cys Cys, 66 Cys, 66 Cys, 65
Tyr Tyr, 66 Tyr, 66 Tyr, 67
COI COI, 1550 COI, 1550 COI, 1592
Ser Ser, 70 Ser, 70 Ser, 70
Asp Asp, 71 Asp, 71 Asp, 72
COII COII, 690 COII, 690 COII, 690
Lys Lys, 74 Lys, 74 Lys, 74
ATP8 ATP8, 167 ATP8, 167 ATP8, 167
ATP6 ATP6, 682 ATP6, 683 ATP6, 683
COIII COIII, 784 COIII, 784 COIII, 783
Gly Gly,70 Gly, 69 Gly, 72
ND3 ND3, 345 ND3, 345 ND3, 348
Arg, Arg, 68 Arg, 69 Arg, 68
ND4L ND4L, 297 ND4L, 296 ND4L, 296
ND4 ND4, 1380 ND4, 1380 ND4, 1380
His His, 64 His, 69 His, 68
Ser Ser, 66 Ser, 66 Ser, 67
Leu Leu, 72 Leu, 72 Leu, 72
ND5 ND5, 1832 ND5, 1838 ND5, 1841
ND6 ND6, 515 ND6, 521 ND6, 521
Glu Glu, 68 Glu, 68 Glu, 68
Cyt b Cyt b, 1137 Cyt b, 1137 Cyt b, 1134
Thr Thr, 71 Thr, 71 Thr, 76
Pro Pro, 69 Pro, 69 Pro, 69
D-loop D-loop, 871 D-loop, 917 D-loop, 1389

Figure 5: Compare the mitochondrial genome arrangement of Labeo with Mastacembelusarmatus [synbranchiformes (NC_023977)], Macrognathus aculeatus [synbranchiformes (NC_027435)], Sicamugilcascasia [mugiliformes (NC_017898)].

  1. When the typical mitochondrial genome arrangement of osteichthyes is compared with that of chondrichthyes such as Chimaera monstrosa, Dasyatiszugei, and Scoliodonlaticaudus (see Figure 27), it was found that the arrangement of tRNA was almost the same. However, a replication origin containing 37 base pairs was observed within the cluster WANCY, more specifically, between tRNAN and tRNAC. So, it can be concluded that in the case of tRNA genes, chondrichthyes differ from osteichthyes.
  2. When the typical mitochondrial genome arrangement of fish is compared with lungfish Neoceratodusforsteri (see Figure 6), it was found that the arrangement of tRNA is the same, but the difference is that in the cluster WANCY, there lies a replication origin of 35 base pairs. So, it can be concluded that in the case of tRNA, lungfish is slightly different from osteichthyes and resembles that of chondrichthyes.
  3. When the typical mitochondrial genome arrangement of fish is compared with amphibian Bufo melanostictus (see Figure 7), it was revealed that the tRNA clusters are different. In amphibians, the tRNA clusters are L1TPF, IQM, WANCY, SD, and HS. A replication origin of 10 base pairs is present between tRNAN and tRNAC. Lungfish also has a replication origin in the same clusters between tRNAN and tRNAC (see Figure 7). The position of the origin of replication is a unique identifying characteristic of lungfish, which is similar to that of amphibian mtDNA and dissimilar to that of other osteichthyesmtDNA.
  4. When the typical mitochondrial genome arrangement of fish is compared with reptiles (G. gekko), birds (A. tristis), and mammals (C. porcellus), it was revealed that the 5 clusters (WANCY, IQM, SD, HSL2, TP) are similar in all cases. However, the position of the TP cluster is different in the case of birds.
  5. By comparing the typical mitochondrial genome arrangement of osteichthyes with that of Drosophila melanogaster, it was revealed that the tRNA clusters of Drosophila are different. In Drosophila, the tRNA clusters ARNS1EP, IQP, WCY, KD, and TP were found (see Figure 9). The same result was found in the sessile barnacle Tetraclita japonica (Begum, R.A. et al., 2004). So, it can be concluded that the mitochondrial genome arrangement of osteichthyes, more commonly known as fish, and invertebrates are far different from each other whenever vertebrate and invertebrate mtDNA is compared.
Typical Osteichthyes Neoceratodus Chimaera Dasyatis Scoliodon
Phe, Phe,  68 Phe, 70 Phe,  71 Phe,  70
12S 12S, 951 12S, 948 12S, 963 12S, 951
Val Val, 70 Val, 70 Val,71 Val, 71
16S 16S, 1679 16S, 1664 16S, 1704 16S, 1669
Leu Leu, 74 Leu, 74 Leu, 74 Leu, 74
ND1 ND1, 971 ND1, 971 ND1, 977 ND1, 974
Ile Ile, 76 Ile, 72 Ile, 68 Ile, 69
Gln Gln, 70 Gln, 69 Gln, 72 Gln, 71
Met Met, 68 Met, 68 Met, 70 Met, 68
ND2 ND2, 1043 ND2, 1041 ND2, 1046 ND2, 1044
Trp, Trp,  71 Trp, 68 Trp, 69 Trp, 70
Ala Ala, 68 Ala, 68 Ala, 68 Ala, 68
Asn Asn, 72 Asn, 72 Asn, 72 Asn, 72
Cys Rep org, 35 Rep org, 52 Rep org, 35 Rep org, 37
Tyr Cys, 70 Cys, 73 Cys, 66 Cys, 66
COI Tyr, 69 Tyr, 67 Tyr, 69 Tyr, 68
Ser COI, 1556 COI, 1559 COI, 1556 COI, 1556
Asp Ser, 70 Ser, 71 Ser, 71 Ser, 70
COII Asp, 68 Asp, 71 Asp, 70 Asp, 68
Lys COII, 712 COII, 690 COII, 69 0 COII, 690
ATP8 Lys, 71 Lys, 72 Lys, 73 Lys, 73
ATP6 ATP8, 167 ATP8, 167 ATP8, 167 ATP8, 167
COIII ATP6, 683 ATP6, 683 ATP6, 683 ATP6, 683
Gly COIII, 785 COIII, 785 COIII, 785 COIII, 785
ND3 Gly, 70 Gly, 68 Gly, 70 Gly, 69
Arg, ND3, 348 ND3, 348 ND3, 348 ND3, 348
ND4L Arg, 69 Arg, 69 Arg, 68 Arg, 69
ND4 ND4L, 296 ND4L, 296 ND4L, 296 ND4L, 296
His ND4, 1380 ND4, 1374 ND4, 1380 ND4, 1380
Ser His, 68 His, 68 His, 68 His, 68
Leu Ser, 66 Ser, 68 Ser, 66 Ser, 66
ND5 Leu, 68 Leu, 71 Leu, 71 Leu, 71
ND6 ND5, 1829 ND5, 1838 ND5, 1844 ND5, 1829
Glu ND6, 515 ND6, 521 ND6, 515 ND6, 521
Cyt b Glu, 68 Glu, 69 Glu, 68 Glu, 69
Thr Cyt b, 1143 Cyt b, 1143 Cyt b, 1142 Cyt b, 1144
Pro Thr, 72 Thr, 70 Thr, 72 Thr, 69
D-loop Pro, 69 Misc_feature, Pro, 69 Pro, 68
D-loop, 925 1533(CRI) D-loop, 2522 D-loop,  1062
Pro, 68
Misc_feature,
1434(CRI)

Figure 6: Compare the mitochondrial genome arrangement of Labeo with other chondrichthyesnamely  Neoceratodusforsteri (AJ 584642), Chimaera monstrosa (AJ 310140), DasyatisZugei (NC_019643), and ScoliodonLaticaudus(KP 336547).

3. rRNA:

Two ribosomal RNA genes were found in the LabeorohitamtDNA, one encoding the large subunit, 16S rRNA (1686 bp), and the other encoding the small subunit, 12S rRNA (954 bp). The tRNAV is present between the 12S rRNA and 16S rRNA genes (see Figure 24), representing a unique arrangement among vertebrate and invertebrate mitochondrial genomes (Begum, R.A. et al., 2004).

  1. Comparing L. rohita with the other two Labeo species, it was observed that the number of base pairs is almost similar, differing only by 1-3 bp (see Figure 3).
  2. Comparing with other osteichthyes, it was found that the number of base pairs differs by up to 22 bp (see Figure 4 & 5).
  3. Similar results were found when comparing with chondrichthyes and lungfish.
  4. When compared with amphibians, it was revealed that the number of base pairs differs by 82 bp (see Figure 7).
  5. When the typical mitochondrial genome arrangement of fish was compared with reptiles (G. gekko), birds (A. tristis), and mammals (C. porcellus), it was found that tRNAVal is present between the small rRNA (12S) and large rRNA (16S) like in other vertebrates and invertebrates. However, the number of base pairs differs by 3 to 25 bp.
  6. Comparing with Drosophila (invertebrate), it was found that there is a large difference between vertebrate and invertebrate rRNA base pair numbers. In Drosophila, 12S rRNA has 785 bp, whereas in vertebrates, it is around 954 bp, and 16S rRNA has 1324 bp in Drosophila, whereas in vertebrates, it is around 1686 bp.
Typical chondrichthyes Typical osteichthyes (Labeo) Neoceratodus(Lung Fish) Bufo melanostictus(Amphibia)
Phe,  70 Phe, Phe,  68 Phe, 67
12S, 951 12S 12S, 951 12S, 936
Val, 71 Val Val, 70 Val, 68
16S, 1669 16S 16S, 1679 16S, 1597
Leu, 74 Leu Leu, 74 Leu, 72
ND1, 974 ND1 ND1, 971 ND1, 960
Ile, 69 Ile Ile, 76 Ile, 70
Gln, 71 Gln Gln, 70 Gln, 68
Met, 68 Met Met, 68 Met, 72
ND2, 1044 ND2 ND2, 1043 ND2, 1032
Trp, 70 Trp, Trp,  71 Trp, 69
Ala, 68 Ala Ala, 68 Ala, 68
Asn, 72 Asn Asn, 72 Asn, 72
Rep org, 37 Cys Rep org, 35 Rep org, 10
Cys, 66 Tyr Cys, 70 Cys, 65
Tyr, 68 COI Tyr, 69 Tyr, 69
COI, 1556 Ser COI, 1556 COI, 1541
Ser, 70 Asp Ser, 70 Ser, 67
Asp, 68 COII Asp, 68 Asp, 68
COII, 690 Lys COII, 712 COII, 687
Lys, 73 ATP8 Lys, 71 Lys, 70
ATP8, 167 ATP6 ATP8, 167 ATP8, 164
ATP6, 683 COIII ATP6, 683 ATP6, 683
COIII, 785 Gly COIII, 785 COIII, 783
Gly, 69 ND3 Gly, 70 Gly, 67
ND3, 348 Arg, ND3, 348 ND3, 339
Arg, 69 ND4L Arg, 69 Arg, 70
ND4L, 296 ND4 ND4L, 296 ND4L, 299
ND4, 1380 His ND4, 1380 ND4, 1364
His, 68 Ser His, 68 His, 68
Ser, 66 Leu Ser, 66 Ser, 66
Leu, 71 ND5 Leu, 68 ND5, 1788
ND5, 1829 ND6 ND5, 1829 ND6, 494
ND6, 521 Glu ND6, 515 Glu, 67
Glu, 69 Cyt b Glu, 68 Cyt b, 1145
Cyt b, 1144 Thr Cyt b, 1143 D-loop, 1967
Thr, 69 Pro Thr, 72 Leu, 71
Pro, 68 D-loop Pro, 69 Thr, 71
D-loop,  1062 D-loop, 925 Pro, 68

Figure 7: Mitochondrial genome arrangement showing the evolutionary relation among fishes and amphibian. Rearrangement shown in arrow marks.

4. D-loop region:

The D-loop region or control region is observed between two tRNA genes, namely tRNAP and tRNAF, in the mitochondrial genome of L. rohita (Accession Number: AB011201).

  1. Comparing the mitochondrial genome of three Labeo species – Labeorohita, L. calbasu, and L. bata (see Figure 3) – it was observed that the position of the D-loop region is the same. Its length differs by only (1-3) bp. Therefore, it can be said that in the case of the D-loop region, the mitochondrial genome arrangement of L. rohita can represent the Labeo species.
  2. When the mitochondrial genome arrangement of L. rohita is compared with that of Tenualosailisha, Channa marulius, and Heteropneustesfossilis (see Figure 4), and with Mastacembelusarmatus, Macrognathus aculeatus, and Sicamugilcascasia (see Figure 5), it was found that although the position is the same, the length of the D-loop region differs from 82 bp to 228 bp, which suggests placing these osteichthyes in different orders. However, from the arrangement comparison, it can be concluded that the arrangement is typical.
  3. When the typical mitochondrial genome arrangement of fish is compared with chondrichthyes fish such as Chimaera monstrosa, Dasyatiszugei, and Scoliodonlaticaudus (see Figure 6), it was found that the arrangement and position of the D-loop region are the same. However, the length differs from (122-1582) bp, which indicates that these fish cannot be placed in the same class. It was observed that they belong to the class Chondrichthyes.
Typical Fish Amphibia Bufo melanostictus Reptile Gekko gecko Bird Acridotheres tristis Mammal Cavia porcellus
Phe, Leu,  71 Phe,  73 Phe, 69 Phe, 65
12S Thr, 71 12S, 957 12S, 979 12S, 942
Val Pro, 68 Val, 69 Val, 70 Val, 69
16S Phe, 67 16S, 1585 16S, 1598 16S, 1565
Leu 12S, 936 Leu, 75 Leu, 75 Leu, 74
ND1 Val, 68 ND1, 962 ND1, 978 ND1, 960
Ile 16S, 1597 Ile, 72 Ile, 72 Ile, 65
Gln Leu, 72 Gln, 72 Gln, 71 Gln, 71
Met ND1, 960 Met, 69 Met, 69 Met, 69
ND2 Ile, 70 ND2, 1035 ND2, 1040 ND2, 1044
Trp, Gln, 68 Trp, 71 Trp, 70 Trp, 70
Ala Met, 72 Ala, 68 Ala, 69 Ala, 69
Asn ND2, 1032 Asn, 72 Asn, 73 Asn, 73
Cys Trp, 69 Cys, 62 Cys, 67 Rep_origin, 37
Tyr Ala, 68 Tyr, 65 Tyr, 71 Cys, 67
COI Asn, 72 COI, 1544 COI, 1551 Tyr, 67
Ser Rep org, 10 Ser, 74 Ser, 73 COI, 1542
Asp Cys, 65 Asp, 67 Asp, 69 Ser, 69
COII Tyr, 69 COII, 692 COII, 684 Asp, 69
Lys COI, 1541 Lys, 68 Lys, 68 COII, 684
ATP8 Ser, 67 ATP8, 169 ATP8, 168 Lys, 67
ATP6 Asp, 68 ATP6, 686 ATP6, 684 ATP8, 204
COIII COII, 687 COIII, 783 COIII, 784 ATP6, 681
Gly Lys, 70 Gly, 69 Gly, 69 COIII, 784
ND3 ATP8, 164 ND3, 339 ND3, 351 Gly, 68
Arg, ATP6, 683 Arg, 68 Arg, 70 ND3, 347
ND4L COIII, 783 ND4L, 296 ND4L, 297 Arg, 68
ND4 Gly, 67 ND4, 1380 ND4, 1378 ND4L, 297
His ND3, 339 His, 71 His, 70 ND4, 1378
Ser Arg, 70 Ser, 64 Ser, 64 His, 69
Leu ND4L, 299 Leu, 70 Leu, 71 Ser, 57
ND5 ND4, 1364 ND5, 1829 ND5, 1818 Leu, 70
ND6 His, 68 ND6, 509 Cyt b, 1143 ND5, 1812
Glu Ser, 66 Glu, 65 Thr, 70 ND6, 531
Cyt b ND5, 1788 Cyt b, 1134 Pro, 70 Glu, 69
Thr ND6, 494 Thr, 70 ND6, 519 Cyt b, 1140
Pro Glu, 67 Pro, 67 Glu, 72 Thr, 69
D-loop Cyt b, 1145 D-loop, 1172 Misc_feature, 1250 Pro, 68
D-loop, 1967 D-loop, 1357

Figure 8: Comparison of mitochondrial genome arrangement of typical fish with other vertebrate such as Amphibian (Bufo melanostictus, A/N: NC_005794), Reptile (Gekko gecko, A/N: HM 370130), Bird (Acridotheres tristis, A/N: NC_015195), Mammal (Cavia porcellus, A/N: AJ222767).number showing the total length of particular gene in the boxes, Protein coding genes are represent bold. tRNA gene are represent as three-letter code.

  1. When the typical mitochondrial genome arrangement of fish is compared with lung fish Neoceratodusforsteri (see Figure 6), it was found that the arrangement and position of the D-loop region are the same. The length differs from that of typical fish by only about 15 bp. This finding suggests that this lungfish can be placed with typical fish in the class Osteichthyes.
  2. When the typical mitochondrial genome arrangement of fish is compared with amphibian Bufo melanostictus (see Figure 7), it was revealed that the D-loop region is 1967 bp long and positioned between the cytochrome b gene and tRNAleu, which differs from the typical mitochondrial genome of fish. Specifically, it can be said that this position of the D-loop region in amphibian Bufo melanostictus is different from that in fish. Therefore, it can be concluded that amphibians belong to a different class.
  3. When the typical mitochondrial genome arrangement of fish is compared with reptiles (G. gecko), birds (A. tristis), and mammals (C. porcellus), it was revealed that the position of the D-loop region of reptiles and mammals is similar to that of typical fish genome. However, the position of the D-loop region in birds is between tRNAglu and tRNAphe. This is due to the transposition of the cyt b gene.
  4. When the typical mitochondrial genome arrangement of fish is compared with Drosophila melanogaster (Invertebrate) (see Figure 9), it was found that the D-loop region of Drosophila is located between the small ribosomal RNA (12S rRNA) and tRNAile, with a length ranging from 14917 bp to 19517 bp. As the position and length vary significantly from vertebrates, it can be concluded that Drosophila should be placed in another subkingdom. The same result was found in the case of barnacles Tetraclita japonica (Rowshan et al., 2004).

Thus, we can conclude that the vertebrate order is first Osteichthyes, then Chondrichthyes, followed by lungfish, and immediately after lungfish is the position of amphibians. Therefore, it can be said that lungfishes are the closest relatives of the class Amphibia.

  The comparison of mitochondrial genome arrangement between typical Osteichthyes (Labeo) and Drosophila.

Figure 9: The comparison of mitochondrial genome arrangement between typical Osteichthyes (Labeo) and Drosophila. Rearrangement are shown in arrow marks (black= a.a; yellow=CO gene; blue=ND gene).

Phylogenetic Tree:

The relationships among Labeo species, Labeo species with other Osteichthyes, Chondrichthyes, lungfish, and amphibians were investigated using a phylogenetic tree based on the nucleotide sequences of CoI, 16S, and cytochrome b genes. The molecular phylogenetic tree with bootstrap values was constructed with particular attention to different groups of animals.

1. Phylogenetic Tree based on Nucleotide Sequence of CoI Gene:

  1. For the nucleotide dataset of the concatenated CoI, the Maximum Likelihood (ML) method produced distinct branching orders. Three species of Labeo formed a monophyletic group with the highest bootstrap value of 100% in the ML tree. This relationship is consistent with the genome arrangement and the phylogeny based on the CoI gene (see Figure 10).
  2. Osteichthyes such as Tenualosailisha, Channa marulius, Heteropneustesfossilis, Mastacembelusarmatus, Macrognathus aculeatus, and Sicamugilcascasia form a monophyletic group according to ML trees, consistent with the phylogeny based on taxonomy.
  3. Neoceratodusforsteri (lungfish) is placed just after the Osteichthyes group in the Maximum Likelihood tree.
  4. The phylogenetic tree reveals that Chondrichthyes fish such as Chimaera monstrosa, Dasyatiszugei, and Scoliodonlaticaudus form a monophyletic group, positioned after Amphibia (Bufo melanostictus).
  5. ML shows that Bufo melanostictus (amphibian) is positioned after the lungfish. This relationship is consistent with the genome arrangement and the phylogeny based on the CoI gene (see Figure 10).
  6. The Maximum Likelihood tree shows that birds are situated just before mammals, with a bootstrap value of 0.93, indicating the relevance of this position. This relationship also supports the taxonomic hierarchy of vertebrates.
  7. In the case of Drosophila, the ML tree shows Drosophila as an outgroup. This finding is consistent with the comparison of mitochondrial genome arrangements.

Phylogenic tree (Maximum Likelihood) based on nucleotide sequence of CO1 gene, number showing the bootstrap values

Figure 10: Phylogenic tree (Maximum Likelihood) based on nucleotide sequence of CO1 gene, number showing the bootstrap values. Nucleotide sequence of Drosophila was used as outgroup.

2. Phylogenetic Tree based on Nucleotide Sequence of 16S rRNA Gene:

  1. Analysis of the phylogenetic tree based on the 16S rRNA gene of 15 vertebrate species revealed that the Maximum Parsimony and Neighbor Joining methods produce different branching orders. However, a common finding in both trees is that Labeo species form a monophyletic group with a bootstrap value of 100.
  2. According to both Maximum Parsimony and Neighbor Joining trees, Tenualosailisha, Channa marulius, Heteropneustesfossilis, Mastacembelusarmatus, Macrognathus aculeatus, Sicamugilcascasia, including Labeo species, form a group, although there are some sub-groups within this larger group. This grouping can be referred to as the Osteichthyes fish group.
  3. The phylogenetic tree reveals that Chondrichthyes fish such as Chimaera monstrosa, Dasyatiszugei, and Scoliodonlaticaudus form a monophyletic group, positioned with the lungfish Neoceratodusforsteri.
  4. Neoceratodusforsteri (lungfish) is placed just after the Chondrichthyes group in the Maximum Likelihood tree.
  5. When a phylogenetic tree (maximum likelihood) was constructed from the 16S nucleotide sequence of different selected species, it was observed that higher vertebrates form a monophyletic group.
  6. In the case of Drosophila, the Neighbor Joining tree shows Drosophila as an outgroup, while the Maximum Parsimony tree shows both Drosophila and Bufo as monophyletic. Bufo’s position is after the fish group (as Osteichthyes, Neoceratodus, and Chondrichthyes are all fish and sequentially situated in the phylogenetic tree, they can be referred to as the fish group).

Phylogenic tree (Maximum Likelihood) based on nucleotide sequence of 16S gene, number showing the bootstrap values.

Figure 11: Phylogenic tree (Maximum Likelihood) based on nucleotide sequence of 16S gene, number showing the bootstrap values. Nucleotide sequence of Drosophila was used as outgroup. CONCLUSION

  1. The genome arrangement of Labeorohita is conserved in most osteichthyes species, indicating a common evolutionary history among these species.
  2. The phylogenetic tree constructed from the nucleotide sequence of 16S rRNA supports the conclusion that lungfish is the closest relative of amphibians.
  3. The arrangement of the 37 genes in L. rohita’s genome could serve as representative data for vertebrates, excluding birds due to a distinct arrangement.
  4. The genome arrangement of invertebrates differs significantly from that of vertebrates, highlighting a clear distinction in evolutionary paths between these two groups.

REFERENCE

  1. Ahmed, M. S., Datta, S. K., &Zhilik, A. A. (2020). Molecular diversity of freshwater fishes of Bangladesh assessed by DNA barcoding. Bangladesh Journal of Zoology, 48(1), 1-19.
  2. Rahman, A. K. A. (2005). Freshwater fishes of Bangladesh (2nd ed.). Zoological Society of Bangladesh, University of Dhaka.
  3. Wolstenholme, D. R. (1992). Animal mitochondrial DNA: structure and evolution. International Review of Cytology, 141, 173-216.
  4. Castro, J. A., Picornell, A., & Ramon, M. (1998). Mitochondrial DNA: a tool for populational genetics studies. International Microbiology, 1(4), 327-332.
  5. Girish, P. S., Anjaneyulu, A. S. R., Viswas, K. N., Shivakumar, B. M., Anand, M., Patel, M., & Sharma, B. (2005). Meat species identification by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) of mitochondrial 12S rRNA gene. Meat Science, 70(1), 107-112.
  6. Bridge, J. S. (1992). A revised model for water flow, sediment transport, bed topography and grain size sorting in natural river bends. Water Resources Research, 28(4), 999-1013.
  7. Boore, J. L., & Brown, W. M. (1994). Complete DNA sequence of the mitochondrial genome of the black chiton, Katharina tunicata. Genetics, 138(2), 423-443.
  8. Boore, J. L. (1999). Animal mitochondrial genomes. Nucleic Acids Research, 27(8), 1767-1780.
  9. Wolfe, K. H., & Sharp, P. M. (1993). Mammalian gene evolution: nucleotide sequence divergence between mouse and rat. Journal of Molecular Evolution, 37, 441-456.
  10. Begum, R. A., Yamaguchi, T., & Watabe, S. (2004). Molecular phylogeny of thoracican barnacles based on the mitochondrial 12S and 16S rRNA genes. Sessile Organisms, 21(2), 47-54.

ACKNOWLEDGEMENT

First and foremost, I extend my deepest gratitude to the Almighty ALLAH for His boundless blessings, guidance, and protection throughout my academic journey.

I am profoundly indebted to my esteemed Professor, Dr. Rowshan Ara Begum, from the Department of Zoology, University of Dhaka, for her unwavering guidance, invaluable insights, and continuous support during the course of this study. Her mentorship has been instrumental in shaping this research.

I extend my heartfelt appreciation to Professor Dr. Reza Md. Shahjahan for his amiable demeanor and generous counsel, which have enriched this work.

Special thanks are due to Professor Mohammad Shamimul Alam Ph.D., for his encouragement, scholarly guidance, and insightful advice, which have been invaluable in the development of this thesis.

I am grateful to my co-supervisor, Associate Professor Hawa Jahan, from the Department of Zoology, University of Dhaka, for her unwavering cooperation, assistance, and friendly disposition throughout this endeavor.

I acknowledge with gratitude the support and guidance provided by Professor Dr. Md. Anowarul Islam, Former Chairman of the Department of Zoology, University of Dhaka.

I extend my deepest appreciation to my parents, whose boundless love, encouragement, and blessings have been my source of strength and inspiration. Their unwavering support has been instrumental in my academic pursuits.

I am also indebted to my family members for their unwavering encouragement and support throughout my academic journey.

Lastly, I extend my sincere thanks to all my friends who have made my time at Dhaka University memorable and enjoyable. I am also grateful to all the members of the laboratory for their assistance and camaraderie.

February,2024                                                                                             Authors

Article Statistics

Track views and downloads to measure the impact and reach of your article.

2

PDF Downloads

4 views

Metrics

PlumX

Altmetrics

Paper Submission Deadline

GET OUR MONTHLY NEWSLETTER

Subscribe to Our Newsletter

Sign up for our newsletter, to get updates regarding the Call for Paper, Papers & Research.

    Subscribe to Our Newsletter

    Sign up for our newsletter, to get updates regarding the Call for Paper, Papers & Research.