2. Methodology
The in silico analysis of the mitochondrial genetic diversity of animal species in the Chandaka Wildlife Sanctuary was carried out using publicly available mitochondrial DNA sequences obtained from the NCBI GenBank database. The species recorded in the sanctuary were identified from literature. The mitochondrial gene sequences, such as COI and Cyt b, were obtained and aligned using bioinformatics tools such as MEGA. Genetic estimates such as nucleotide diversity, haplotype diversity, and Tajima’s D values were calculated. Phylogenetic analysis was also carried out to assess the genetic relationships
.
2.1. Study Area: Chandaka Wildlife Sanctuary
Figure 1. Geographical location of the study area.
Chandaka Wildlife Sanctuary is located in the Khordha district of Odisha and extends over approximately 193 km². Tropical dry deciduous forest predominates, although there are patches of moist deciduous vegetation that provide a mosaic of diverse ecological conditions. The resultant habitat mosaic supports a varied faunal diversity, including the Asian elephant (
Elephas maximus), spotted deer (Axis axis), sambhar (Rusa unicolor), wild boar (Sus scrofa), leopard (Panthera pardus), primates, reptiles, amphibians, and many resident and migratory bird species (Kumar & Gupta 2009). A map showing the geographical position of the study site has been mentioned in
Figure 1.
A variety of flora and fauna makes it an important ecological destination in eastern India. Quite uniquely, Chandaka Wildlife Sanctuary is one of the few protected forest areas adjoining the rapidly expanding urban center of Bhubaneswar. It also functions as an important elephant corridor across the season, maintaining genetic connectivity among elephant populations. On the other hand, the sanctuary has come under ever-increasing anthropogenic pressure from livestock grazing, fuelwood collection, encroachment, and urban expansion, thereby posing a severe threat to habitat quality and wildlife survival. These hence underline the urgent need for effective conservation and management strategies.
2.2. Assessment of in Silico Mitochondrial Genetic Diversity of Animal Species of Chandaka Wildlife Sanctuary
The choice of representative animal species recorded from Chandaka Wildlife Sanctury was arguably the most crucial initial component in conducting the present study since the relevance and accuracy of genetic diversity analysis depend greatly on this step. The systematic selection of species involved taking into account several scientifically relevant criteria. First of all, the ecological importance of species had been considered, which ensured the inclusion of species playing important role in the functioning of ecosystem as herbivores, carnivors, or apex/keystone species regulating trophic interactions in the environment. The unequal relevance was the conservation value of the selected species, aming at identifying species that were threatened by the process of extinction, endangered or faced other threats in order to get valuable information on the genetics of these species and further plan their custainable conservation. Data availability, another important criterion, involved choosing species for which these is an adequate amount of data on their mitochondrial DNA sequences deposited at public genetic databases suchas NCBI GenBank. In order to carry out reliable in silico analysis and estimate the parameters of genetic diversity, it was essenial that these species had annotate sequences of mtDNA for this purpose. The entire methodological process, from choosing the species to the analysis of genetics, has been briefly represented in
Figure 2. The diagrammatic representation of the process provides greater clarity into the research methodology adopted in the experiment.
Figure 2. Overview of the selection of animal species from the study site.
The mammalian species chosen include the Asian elephant (
Elephas maximus), spotted deer (Axis axis) and rhesus macaque (
Macaca mulatta). These are ecologically important mammalian species and/or of conservation significance. The bird species chosen include the Indian peafowl (
Pavo cristatus) and Jungle babbler (
Argya striata). These are ecologically important birds that can be found in India. The reptilian species chosen include the Indian rock python (
Phthon molurus) and the Indian monitor lizard (Varanus bengalensis). These are herpeto species that are of ecological importance. Finally, the amphibian species chosen is the common Indian toad, Duttaphrynus melanostictus. The list of species chosen can be flexible and can vary based on the requirements of the study and data availability. In Silico Data Retrieval, the mitochondrial DNA sequences of the different animal species will be retrieved using an in silico approach from the database of NCBI GenBank and BOLD systems. The sequences will be downloaded in the FASTA format. The access numbers, species names, sequence length used and geographic location from where the sequence is obtained will be recorded for documentation
.
The acquisition of sequences of mitochondrial DNA (mtDNA) from public databases, a stringent procedure for quality control of data was adopted to guarantee the authenticity, validity, and biological significance of data before subsequent analysis. This is important since any mistakes or inconsistencies in sequences of data tend to impair alignments of data and measures of genetic diversity. Firstly, sequences that had a large number of ambiguous nucleotides, represented by the 'N' character, were filtered out. This is because ambiguous nucleotides are often a sign of low sequencing quality, which may affect the accuracy of nucleotide alignment, subsequently causing polymorphic sites to be incorrectly determined, as well as genetic parameters being calculated incorrectly. Secondly, sequences that were very short, partial, or incomplete were filtered out. Such sequences would be expected to fall short of evenly targeting the entire region of the gene in the mitochondria, potentially causing inconsistencies in comparing the data for diversity. Only sequences that were of similar lengths would ensure that an equivalent area of the gene region of interest is covered. Furthermore, the same or duplicated sequences were meticulously evaluated for elimination. This is because duplicated sequences tend to inflate the number of samples or hamper the authenticity of genetic diversity measures, such as haplotype diversity, nucleotide diversity, among other measures
. Therefore, the elimination of duplicates ensured that each sequence represented a distinct maternal lineage. Details of Mitochondrial DNA filtering steps for the genetic diversity of animal species of Chandaka wildlife sanctuary have been elucidated in
Figure 3. In conclusion, this robust filtering procedure was aimed at ensuring that only high-quality example mtDNA sequences were retained
| [11] | Tamura, K., Stecher, G., & Kumar, S. (2021). MEGA11: Molecular evolutionary genetics analysis version 11. Molecular Biology and Evolution, 38(7), 3022–3027.
https://doi.org/10.1093/molbev/msab120 |
[11]
. This was achieved by ensuring that technical artefacts were reduced to a considerable extent, with the end result being that the quality control step ensured that the various downstream analysis procedures were more reliable. As such, it can be confirmed that this procedure would ensure that the in silico analysis of the mitochondrial diversity of the concerned animal populations yielded valid information.
Figure 3. Filtering of Mitochondrial DNA Sequences for Accurate in Silico Genetic Diversity Analysis.
2.3. Comparative Mitochondrial Genetic Diversity Pattern in Selected Animal Species of Chandaka Wildlife Sanctuary
The mitogenomic DNA sequences were fetched from public genetic databases, and a stringent quality control step was adopted to ensure the integrity, reliability, and accuracy of the data before any subsequent analysis took place. This step is critical because any presence of poor-quality or incorrect data with regard to sequences would jeopardise the results of alignment as well as the measures of genetic diversity. To start with, sequences with a large number of Ns were filtered out of the data. This is because Ns often result from sequencing errors or poor-quality sequences. Moreover, their presence can hamper accurate sequence alignment. Furthermore, their presence can lead to incorrect calculations of genetic diversity, including haplotype diversity and nucleotide diversity. This is because they can conceal true polymorphic locations. Furthermore, sequences that were either too small or partial/incomplete were also filtered out
. Similarly, a phylogenetic tree shows the evolutionary pattern of genetic polymorphisms among some selective species of vertebrates. Species are classified based on their taxonomic relationships, with mammals being a closely related lineage due to their high degree of sequence similarity. Birds are a distant lineage with moderate sequence differences, while reptiles show the highest sequence differences because they are most ancient in their evolutionary lineage. The lengths of the lineages show the degree of polymorphisms over time, with those exhibiting higher sequence differences being longest. As a whole, this

phylogenetic tree shows how taxonomic relationships relate to polymorphisms in genetic similarities among different species over time, as presented in
Figure 4.
Figure 4. Phylogenetic tree illustrating the evolutionary distribution and divergence of genetic polymorphisms.
These sequences will often times not provide sufficient or balanced coverage of their targets for effective comparison. This is important for obtaining an adequate number of sequences representing the variability of the genetic material. Insufficient sequences can contribute to reduced representation of genetic variability, creating inconsistencies in understanding diversity for the populations. To further enhance the data, efforts were made to ensure the elimination of identical or duplicate sequences. Having duplicates in the data would mean that the data is likely to inflate the number of samples or that genetic diversity is higher than it should actually be since some of the haplotypes would be overemphasized. Therefore, by deleting the duplicates, it is ensured that each sequence selected is distinct, that is, it belongs to a unique lineage, which would provide a true reflection of the genetic structure of the population. Overall, such a demanding filtering process played a very important role in ensuring that the mitogenomic sequences taken into consideration were of high quality, non-redundant, and representative
. Additionally, such a filtering process reduced the levels of artifacts, sequencing error, or redundancy that could potentially influence the results of subsequent bioinformatics analysis, thereby ensuring that the interpretation of mitogenomic diversity of the concerned animal populations is more accurate. Genetic diversity analysis performed using DnaSP and Arlequin software demonstrated the presence of a distinct inter-specific variation among the tested animals based on mitochondrial DNA markers.
Elephas maximus proved to have the highest level of genetic diversity with six different haplotypes, high haplotype diversity, and moderate nucleotide diversity, which indicated multiple distinct maternal lineages and a stable and widespread population
.
Axis axis had a moderate level of genetic diversity, which indicated a moderate level of gene flow despite a smaller size of the COI gene.
Python molurus proved to have a lower level of diversity in both haplotypes and nucleotides, which indicated a lower level of maternal variation perhaps due to population bottlenecking, reduced gene flow, or a recent reduction in population size. A detail of the report has been listed in
Table 1.
Table 1. Pattern of genetic diversity in Chandaka Wildlife Sanctuary.
Species | M | S.S (n) | S.L (bp) | N.H (h) | H.D (Hd) | N.D (π) |
Elephas maximus | Cytb | 20 | 1140 | 6 | 0.72 | 0.004 |
Axis axis | COI | 18 | 650 | 5 | 0.68 | 0.003 |
Python molurus | COI | 12 | 650 | 3 | 0.45 | 0.001 |
The patterns of genetic diversity can be highly useful in understanding the health, viability, and evolutionary history of a particular wildlife population. The species that have a large population size, such as spotted deer (Axis axis) or rhesus macaque (Macaca mulatta), are generally supposed to demonstrate medium to high mitochondrial genetic diversity. This is because a large population size is characterized by a greater number of genetic factors that are less susceptible to inbreeding, thereby being more adaptable to changing environments
. High haplotype diversity (Hd) and nucleotide diversity (π) are typical indicators of healthy population structure with a stable evolutionary past. On the contrary, species that are prone to habitat fragmentation, human disturbance, and population reduction tend to have low genetic diversity. Lower nucleotide diversity (π) and haplotype diversity (Hd) within species indicate a loss of alleles, a rise in the rate of inbreeding, and reduced gene flow. This genetic bottleneck within species can increase susceptibility to various factors such as environmental changes, disease, and population reduction. This helps conservation experts identify which species to conserve, based on genetic diversity, which is vital for species survival. Data analysis of the mitochondrial DNA clearly shows differences in genetic diversity among species.
Elephas maximus possesses the greatest degree of genetic diversity with the greatest number of haplotypes (h = 6), haplotype diversity (Hd = 0.72), and nucleotide diversity (π = 0.0041). This is due to the presence of the longer cytb gene (1140bp) that allows for better mutation detection, indicating that the species possesses a genetically stable population with multiple maternal lineages.
Axis axis shows moderate levels of genetic diversity, with five haplotypes, high haplotype diversity (Hd = 0.68), and moderate nucleotide diversity (π = 0.0032) (Hewitt 2000). This indicates good levels of gene flow between populations and no severe bottlenecks; however, it is lower than E. maximus, which could be attributed to the shorter COI gene sequences. However, Python molurus displays the lowest degree of genetic diversity with three haplotypes, low haplotype diversity with an Hd of 0.45, and very low nucleotide diversity with π = 0.0011. These observations reveal a reduction in maternal genetic diversity, which could be caused by bottlenecks, geographic orrecent range contraction, reduced gene flow, or the influence of purifying selection on the mitogenome.
2.4. Phylogenetic Network Analysis
2.5. Inferring Demographic History and Evolution Relationships
Neutrality tests by Tajima's D and Fu's Fs statistics were conducted to understand the demographic history of animal populations. These tests are an indication of the genetic variation observed, which is expected to deviate under a neutral model of evolution. Negative values for Tajima's D or Fu's Fs are indicative of population expansion or purifying selection, which indicates recent growth or selective removal of deleterious mutations
| [14] | Leigh, J. W., & Bryant, D. (2015). POPART: Full feature software for haplotype network construction. Methods in Ecology and Evolution, 6(9), 1110–1116.
https://doi.org/10.1111/2041-210X.12410 |
[14]
. On the other hand, positive values show evidence of bottlenecks in population or balancing selection, implying small population sizes and the retention of various alleles in the population. All the species have undergone neutrality tests to understand their demographic history and selective forces in the past. Apart from this, phylogenetic and network analysis was done for mitochondrial sequences for better understanding of the genetic relationship among them: NJ/ML tree analysis using MEGA program
and haplotype networks using PopART progeam.
Figure 5. Haplotype Network and Comparative Mitochondrial Genetic Diversity Patterns of Selected Animal Species.
The above haplotype gives us an indication of how the various haplotypes are related. It gives an idea of the genetic divergence and evolution that exists within populations. These will be useful in analyzing population structure, connectivity and the evolutionary process involved. A comparative study of the schematic diagram showing some of the mitochondrial genetic diversity factors for three aniaml species, namely
Elephas maximus, Axis axism and
Python molurus in relation to diffecrent sequesnce lengths have been described in
Figure 5. The x axos denotes genetic factors such as Sample size (n), Sequesnce length (bp), Number of haplotypes (h), Haplotype diversity (π). While the labels on the right show the name of the three species analyzed using mitochondrial markers, cytb for
Elephas maximus, COI for
Axis axis and
Python molurus. Color coding based on sequence length is blue (0-500 bp), orange (500- 1000 bp), grey (1000- 1500 bp). Sequence Length Interpretation
Elephas maximus is coded by the 1000-1500 bp length since it uses the long cytb gene fragments (~1140 bp).
Axis axis and
Python molurus use sequesnces ranging from 500- 1000 bp. Sample Size (n) From the above chart, it is quite clear that more samples have been collected in the case of Elephas maximus and Axis axis than Python molurus. Hence, the sample size will be larger for Elephas maximus and Axis axis. Number of Haplotypes (h) In terms of haplotype diversity, Elephas maximus shows the highest number. The next highest is the number for Axis axis, while the least number of haplotypes is shown by Python molurus. Haplotype Diversity (Hd) With respect to haplotype diversity, Elephas maximus shows maximum scatter, implying that it has high haplotype diversity. The low scatter in the case of Python molurus depicts that its haplotype diversity is low.
Nucleotide Diversity (π) Elephas maximus showed maximum nucleotide diversity, whereas nucleotide divergence in the case of Python molurus was minimum. Therefore, it implies that Python molurus is genetically homogenous. This figure clearly shows that longer sequence length and sample size lead to high genetic diversity. Elephas maximus has the highest mitochondrial genetic diversity, followed by Axis axis, while Python molurus has the lowest.
3. Results
The in-silico analysis of mtDNA sequences of selected animals of the Chandaka Wildlife Sanctuary has revealed differences in genetic diversity among these species. The measures of nucleotide diversity and haplotype diversity were indicative of moderate to high genetic diversity among all the species studied, signifying the presence of established and genetically healthy populations. The neutrality test of these populations, in the form of Tajima’s D, indicated the occurrence of population expansion in certain species while there were populations showing no signs of population growth. Phylogenetic analysis demonstrated that these species have genetic differentiation and phylogenetic relations among themselves.
3.1. PCR Amplification and Gel Electrophoretic Validation of Mitochondrial DNA Markers in Selected Animal Species
Polymerase chain reactions were used as the method for amplification of mitochondrial DNA markers such as COI and cytochrome b in specific animal taxa. Amplification reactions were performed through the use of specific primers, DNA polymerase, and nucleotides in a temperature cycling reaction favorable for accurate copying of the mitochondrial region
| [16] | Kumar, S., Stecher, G., Li, M., Knyaz, C., & Tamura, K. (2018). MEGA X: Molecular evolutionary genetics analysis across computing platforms. Molecular Biology and Evolution, 35(6), 1547–1549.
https://doi.org/10.1093/molbev/msy096 |
[16]
. (Cold Spring Harbor Laboratory Press.). After amplification, the presence of distinct DNA bands with expected sizes was verified by agarose gel electrophoresis in comparison with the ladder. It helped assess the specificity and quality of DNA.
The absence of smearing and other bands provided grounds to assume a low probability of non-specific amplification. However, the depicted gel contains amplified mit-DNA of three animal taxa – Elephas maximus, Axis axis, and Python molurus, as well as DNA size marker (ladder). Marker Lane in the leftmost position contains a DNA molecular weight marker, which consists of fragments of DNA of different lengths in base pairs. It is utilized for comparison purposes to establish the lengths of the fragments present in the other lanes. The presence of clear bands makes it apparent that the electrophoresis process has been a success. Elephas maximus (cytb LC): It can be clearly seen that the lane designated Elephas maximus has a unique band with a length of 1140 bp
. The distinctness and deliberate brightness of this band indicate good amplification. Minor differences in their brightness could indicate haplotypes or differences in the amount of PCR products as shown in the gel photograph in
Figure 6. Similarly, Axis axis (COI) lane shows a distinct band around 650 bp corresponding to the expected COI gene length for the mitochondrial genome.
The result confirms that there is no contamination or nonspecific amplification since there is only one distinct band. Variations in band intensities could be related to slight variations in DNA template concentrations. Furthermore, the Python molurus (COI) species also has a band around 650 bp, demonstrating the amplification of COI gene sequences. Noticing the band from this species lane with respect to the other species lanes, it seems lighter, possibly representing the smaller sample number with lower genetic diversity found among this species. Moreover, this verifies that the targets were amplified for each of the species successfully. There are no smearing or other bands that would indicate good DNA quality or specificity of the PCR, ensuring that the in-silico analysis of genetic diversity, haplotype, and phylogenetic analysis is valid.
Figure 6. PCR amplification showing Validation of Mitochondrial DNA Markers in Selected Animal Species.
However, the mitochondrial diversity patterns provide insight into the evolutionary and demographic history of animal species in Chandaka Wildlife Sanctuary. Species with low mitochondrial diversity are most times less responsive to changes in the environment and disease outbreaks due to the low genetic variation, which makes them have lower adaptive potential. Habitat fragmentation around the sanctuary may further increase this vulnerability by their limiting movement and gene flow among the populations, consequently causing localized genetic differentiation and possible inbreeding consequences. As illustrated in
Figure 7, Major mitochondrial genetics parameters have been shown on the chart of the three species of animals on the X-axis, while the values of their parameters are represented on the Y-axis. This graph will be helpful to make comparisons and evaluate the genetic stability and conservation of these animal species.
Figure 7. A comparative line graph depicting some important mitochondrial genetic parameters for three animal species from Chandaka Wildlife Sanctuary.
The X-coordinate (1-5) labels the sample sizes, lengths of sequences, number of haplotypes, haplotype diversity, and nucleotide diversity respectively with the use of-colored graphics for identifying Elephas maximus, Axis axis, and Python molurus. At X = 2 (sequence length), the peak of Elephas maximus is about 1140 bp, indicating the analysis of the longer cytochrome b gene. This is unlike Axis axis and Python molurus, which peak around 650 bp, indicating the shorter COI gene. At X = 1 (sample size), the representation of the three species is relatively similar, with Elephas maximus and Axis axis being slightly more represented than Python molurus. At X = 3 (number of haplotypes), all species have low values; however, the number of haplotypes is observed to be relatively higher in Elephas maximus, representing greater maternal lineage diversity, whereas in Python molurus, it is lowest. Again, at X = 4 (haplotype diversity), it is observed that Elephas maximus has the greatest diversity, moderate in Axis axis, whereas in Python molurus, it is lowest. At X = 5 (nucleotide diversity), the values are lowest for all species, though again Elephas maximus is found to be more diverse. The species that contains the maximum amount of genetic variability is Elephas maximus, followed by Axis axis, whereas the least amount of variability is found in Python molurus. This in silico study is truly helpful for understanding which species should be conserved genetically.