The Crocodylus Genus And The Gavialidae Family Biology Essay
Crocodile relationships are frequently discussed and debated to to where some groups need to be classified. Another subject is the preservation position of some of these cryptic animals. These are both subjects that are touched on in this survey which was done to find the relationship between the genus Crocodylus and critically endangered Gavialis gangeticus which was until late considered portion of the same household as Crocodylus. Next I looked with in the Crocodylus genus to see from where did these species originate as this will assist to conserve the single species and this will assist under base which species are likely able to crossbreed as this will hold grieve effects. The sequence alliance and information analyses was all done in MEGA 5 version 5 which use Genbank as its informations beginning. Genbank being a unfastened beginning it was difficult happening the right sequences. But our consequence that I obtained from these trial revieled that Gavialis gangeticus demands to be in its ain household and that there is a few Crocodylus species that do crossbreed particularly in imprisonment. So we need a much better attempt to conserve our zoology of crocodiles as many are on the IUCN ruddy list of threatened species.
Cardinal words: Crocodylus, G. gangeticus, Phylogenetic relationships of crocodiles, Neighbour-Joining, Minimum Evolution, Maximum Parsimony, UPGMA
Using molecular systematics have become popular to qualify species and their related households. It is besides a helpful tool in assisting preservation in many instances. By utilizing molecular systematics one can easy and with really high certainty classify persons in to the right household and even sub specie degree. So by agencies of these systems we are traveling to see the relationship with in the genus Crocodylus and the relationship between it and Gavialidae. By under standing the relationships one can make preservation work on each single cognizing its relationship to its relations which in most instances portion common home grounds. Therefore we had a expression at “ The IUCN ruddy list of threatened species ” , to see which of the Crocodylus genus are under force per unit area. Here follows the Scientific name ( Common name ) Status: ; Crocodylus acutus ( American crocodile ) Status: Vulnerable A1acA ver 2.3, Crocodylus intermedius ( Orinoco crocodile ) Status: Critically Endangered A1c, C2aA ver 2.3, Crocodylus johnsoni ( Johnstone river crocodile ) Status: Lower Risk/least concernA ver 2.3, Crocodylus mindorensis ( Philippine crocodile ) Status: Critically Endangered A1c, C2aA ver 2.3, Crocodylus moreletii ( Morelet ‘s crocodile ) Status: Lower Risk/conservation dependentA ver 2.3, Crocodylus niloticus ( Nile crocodile ) Status: Lower Risk/least concernA ver 2.3, Crocodylus novaeguineae ( New Guinea crocodile ) Status: Lower Risk/least concernA ver 2.3, Crocodylus palustris ( Marsh crocodile ) Status: Vulnerable A1a, C2aA ver 2.3, Crocodylus porosus ( Australian seawater crocodile ) Status: Lower Risk/least concernA ver 2.3, Crocodylus rhombifer ( Cuban crocodile ) Status: Critically Endangered A2cdeA ver 3.1, Crocodylus siamensis ( Siamese crocodile ) Status: Critically Endangered A1acA ver 2.3 ( IUCN 2012 ) . The out group position are besides available and is every bit follows ; Osteolaemus tetraspis ( dwarf crocodile ) Status: Vulnerable A2cdA ver 2.3, Gavialis gangeticus ( Gharial ) Status: Critically Endangered A2bc ; C1A ver 3.1A and diminishing in the natural state ( IUCN 2012 ) . I had a expression at the households and genus degree of these persons an specifically chose the two outgroups as one with in the household ( Osteolaemus tetraspis ) and on out side the household ( Gavialis gangeticus ) of the genus Crocodylus, Crocodylinae ( Janke 2008, Li et Al. 2007 )
Crocodile have a really interesting evolution as there is still a really hot argument on traveling to which households which species belong to ( Janke 2008, Li et Al. 2007 ) . I looked at the relationship with in the Crocodylus genus as a really of import preservation group as seven out of 11 species are either: Lower Risk/Conservation dependant, Vulnerable or Critically Endangered. This is besides taking in to consideration where these persons find themselves where most are fring their natural home ground at a really rapid rate ( Mazzotti et al. 2009 ) . Most of these animate beings where and are hunted or breed for their tegument as a leather merchandise ( Blake & A ; Loveridge 1975 ) . In some countries as the Everglades crocodiles are the flagship species of preservation and the success thereof is measured as how good these crocodiles to and if they are reproducing of course ( Mazzotti et al. 2009 ) . Other country where the preservation of crocodiles is of import is India, Indonesia, Cuba and Malaysia. In these countries the crocodiles are under high force per unit area as all of them are endangered ( IUCN 2012 ) . So how can this aid conservation even more? Well in offenses where illegal trading takes topographic point we are able to see if it is an endangered coinage or non that is being trade with. In the instance where the coinage is at danger or vulnerable the informations can be used as grounds to demo to which species the individuals are merchandising in. most states can so implement the right statute law and Begin with preservation programs. Some of these states include Zimbabwe where there was a immense success in re-populating the Nile crocodile and the on traveling attempts in India ( de Vos 1984, McGregor 2005 )
So I so looked ate the relationship of the Gavialis gangeticus in relation to Crocodylus as it is been speculated that they are non near relations and necessitate to be separate households ( Janke 2008, Li et Al. 2007 ) . I supported this premise made in past literature as some work has shown that this is true ( Janke 2008, Li et Al. 2007 ) . My work would assist to demo that this is true as I expect it to be the lowest and first divergency on any of the trees drawn. This would assist set up how large the impact is on the bing populations and engendering plans can be worked out to see a really healthy and stable population with a large assortment of allelomorphs.
Crocodiles show a really great diverseness and we should seek and maintain this as each coinage plays a really of import function in our natural ecosystem. From the big Saltwater crocs to the midget crocodiles and even their cousins the cayman and alligators.
Materials and Methods
Phylogenetic and molecular analyses of the development of the Crocodylus genus was done by utilizing MEGA version 5 ( Tamura et al. 2011 ) a powerful calculating plan. This plan, MEGA 5, uses Genbank as its familial database for seeking and happening appropriate familial sequences. These sequences are so added to MEGA 5 and so aligned with in the plan itself. MEGA 5 can be used to work out some statistical values and to build phyletic trees with the aligned informations with was collected from Genbank. The trees that were drawn in this plan were Neighbour-joining ( NJ ) , Minimum Evolution ( ME ) , Maximum Parsimony ( MP ) and UPGMA ( UP ) .
Ingroup and outgroup choice
I chose the following groups as my cliques: Crocodylus acutus ( American crocodile ) , Crocodylus intermedius ( Orinoco crocodile ) , Crocodylus johnsoni ( Johnstone river crocodile ) , Crocodylus mindorensis ( Philippine crocodile ) , Crocodylus moreletii ( Morelet ‘s crocodile ) , Crocodylus niloticus ( Nile crocodile ) , Crocodylus novaeguineae ( New Guinea crocodile ) , Crocodylus palustris ( Marsh crocodile ) , Crocodylus porosus ( Australian seawater crocodile ) , Crocodylus rhombifer ( Cuban crocodile ) , Crocodylus siamensis ( Siamese crocodile )
There was merely two outgroups chosen:
Osteolaemus tetraspis ( dwarf crocodile ) , Gavialis gangeticus ( Gharial )
These two outgroups were chosen as they are portion of different genus each and would give a really good presentation of the Crocodylus genus and their relation ship with in the genus itself. They were besides chosen as at the following systematic degree Crocodylinae is the sub-family degree to which O. tetraspis belong and G. gangeticus belong to the a household Gavialidae which is a household outside Crocodylinae, and the others that was tested for where all in the clique to which the genus is Crocodylus.
Sequence for alliance
The full chondriosome genome was chosen to acquire a much better image of the relation ship between the genus and household degrees as the all contain around 16100 base braces which could be compared to one another. This would give me much more variable sites to work with and more parsimony enlightening sites. These sequences were al taken from natural or confined animate beings therefore no Deoxyribonucleic acid debasement would hold taken topographic point.
Constructing an aligned information set for analyses
In MEGA 5 when piecing the informations it opens Genbank from within MEGA 5. To acquire the hunt traveling one demand to infix Crocodylus AND full mt-genome in the hunt check at the top of the screen. I so searched for a full mt-genome sequence from Crocodylus niloticus ( Nile crocodile ) and so added it to the information set in MEGA 5. I used the blast option within MEGA 5 to blare against the C. niloticus sequence. Then sequences near to but non precisely the same from different species were selected and added to the information set. The two outgroups I searched manually for by come ining their scientific names AND full mt-genome in the hunt check. I merely picked one sequence from every entry.
The sequence statistics was obtained by bring forthing variable sites ( V ) and Parsimony-informative sites ( Pi ) which so was recorded for the set of informations that we obtained. Some nucleotide composings were besides generated these were the mean values for the T, C, A, and G composing in the information set. ‘R ‘ statistic, transition/transversion ratio, could be generated utilizing the nucleotide brace frequence trial. p-distance was use as the theoretical account to bring forth and obtain the overall average distance of the nucleotides development.
The theoretical account choice could merely been done after all the appropriate statistical information was obtained some that the right theoretical account of sequence development could be generated and selected. This theoretical account would be used in the phonic analyses and was done by bring forthing the hierarchal likeliness theoretical account choice trial and choosing via the AICc ( Akaike Information Criterion, corrected ) . The most disposed theoretical account was the TN93+G+I for the information set at manus.
With all the obtained information in this information set we were able to build a Neighbour-joining ( NJ ) , Minimum Evolution ( ME ) , Maximum Parsimony ( MP ) and UPGMA ( Unweighted-Pair Group Method with Arrhythmic agencies ) ( UP ) tree. All the SBL values were obtained for the appropriate trees. The MP tree was constructed and the RI, CI and CRI values was obtained from MEGA 5 for the different homoplasy indexes. The molecular clock was tested by agencies of Tajima ‘s trial and the log likelihood mark was obtained this was done to prove the heterogeneousness.
When the information set was aligned and trimmed I found that there 16311 bases in length that was composed by 13 taxa of which 5734/16311 were variable sites ( V ) and 2973/16311 of variable sites that were parsimony enlightening ( Pi ) . In the nucleotide base composing consequence I found that the mean base braces were T=24.8 ; C=28.8 ; A=31.5 ; G=14.9 with a transition/transversion ratio ( R statistic ) of 2.99 and the overall average distance measured in ‘d ‘ statistic was 0.108 ( 10.8 % ) . In the consequences of the hierarchal likeliness theoretical account choice trial I used the AICc to take the best disposed theoretical account and this was accomplished by seeking for the theoretical account with the lowest value. The TN93+G+I Model of sequence development ( Tamura-Nei ) was the chosen theoretical account as it had the lowest value. The Invariant parametric quantity value was 0.374 and the Gamma parametric quantity value was 0.569. Neighbour-Joining ( NJ ) ( Fig. 1 ) , Minimum Evolution ( ME ) ( Fig. 2 ) , Maximum Parsimony ( MP ) ( Fig. 3 ) and a UPGMA ( UP ) ( Fig. 4 ) tree was drawn and all showed similar consequences with a little fluctuation in the bootstrap values. These fluctuations were really little as all the subdivisions of the trees with a bootstrap value less than 50 % were collapsed. In the NJ tree ( Fig. 1 ) one can clearly see that the Crocodylus genus is split into two major lineages which diverged in to seven bombers lineages. The two chosen out groups are still as expected on their ain ( Fig. 1 ) . The Sum of Branch Lengths ( SBL ) was 0.84424731 for the NJ tree ( Fig. 1 ) . The ME tree ( Fig. 2 ) was constructed on a SBL that was 0.84386914. This tree ( Fig. 2 ) had a really similar form and the lone difference was some bootstrap values. In the MP tree ( Fig. 3 ) there was a spot of uncertainness around the arrangement of C. moreletii but it still was placed in the same lineage as the NJ ( Fig. 1 ) an ME ( Fig. 2 ) trees had placed it. The MP tree concepts many trees and chooses the 1 with the least sum of alteration in it. I got two trees that were really similar but I chose the best fit tree by looking which had the higher bootstrap values. This tree ( Fig. 3 ) besides split the Crocodylus in to seven separate lineages which shows us the same as the NJ tree ( Fig. 1 ) and ME tree ( Fig. 2 ) . It did nevertheless trade the out groups around. I obtained the different homoplasy indexes for the MP tree and these were as follow: consistence index ( CI ) = 0.653600 ; keeping index ( RI ) = 0.471171 ; composite index = 0.307957 ( for all sites ) . For the UP tree I foremost ran a clock trial utilizing Tajima ‘s comparative rate trial which was non rejected at a p-value of 0.32302 with 1 grade of freedom therefore intending I could utilize the UP tree ( Fig. 4 ) . I ran a 2nd molecular clock utilizing the Likelihood trial for the rate of heterogeneousness ; this was rejected though therefore a molecular clock could non be carried out. I was non able to transcend 10 000 bootstrap reproductions, as the plan automatically changed it from 100 000 reproductions to 10A 000 when come ining 100 000.
Figure 1. Neighbour-Joining Tree. This is the optimum tree for the Neighbour-Joining tree and has a subdivision length amount of 0.33980092 ( Saitou and Nei 1987 ) . This tree is the consequence of 10 000 bootstrap reproductions and the per centum of trees where the associated taxa clustered together is shown by the values next to the subdivisions ( Felsenstein 1985 ) . Bootstrap values less than 50 % were collapsed. Tamura-Nei method was used to calculate the evolutionary distances ( Tamura 1992 ) . A gamma distribution ( gamma parametric quantity value = 0.57 ) was used to pattern the rate of fluctuation amongst sites. Thirteen nucleotide sequences were used in the analysis of the Neighbour-Joining tree ( Tamura et al. 2011 ) . The concluding dataset had a sum of 16311 places ( Tamura et al. 2011 ) .
Figure 2. Minimal Evolution tree. The Neighbour-Joining tree was used to bring forth the initial tree. This tree is the consequence of 10 000 bootstraps. Branchs with less than 50 % bootstrap values are collapsed. The per centum associated with the taxa are clustered together are shown following to the subdivisions. Tamura-Nei method was used to calculate the evolutionary distances ( Tamura 1992 ) . The rate fluctuation among sites was modelled with a gamma distribution ( form parametric quantity = 0.57 ) . Thirteen nucleotide sequences were used in the analysis. There were a sum of 16311 places in the concluding dataset.
Figure 3. Maximal Parsimony Tree. There were two most penurious trees, of which the above tree is tree figure one with the best bootstrap cut off values. This tree is the consequence of 10 00 bootstrap reproductions and the per centum of trees where the associated taxa clustered together is shown by the values next to the subdivisions ( Felsenstein 1985 ) . Bootstrap values less than 50 % were collapsed. Thirteen nucleotide sequences were used in the analysis of the Maximum Parsimony tree ( Tamura et al. 2011 ) . The concluding dataset had a sum of 16 311 places ( Tamura et al. 2011 ) .
Figure 4. UPGMA ( Unweighted-Pair Group Method with Arrhythmic agencies ) Tree. The Neighbour-Joining method was indirectly used to build the evolutionary history of the UPGMA tree ( Saitou and Nei 1987 ) . This tree is the consequence of 10 000 bootstrap reproductions and the per centum of trees where the associated taxa clustered together is shown by the values next to the subdivisions ( Felsenstein 1985 ) . Bootstrap values less than 50 % were all collapsed. The tree is drawn to scale measured in MYA ( Millions of Years Ago ) with intervals of 0.05 MYA. Tamura-Nei method was used to calculate the evolutionary distances ( Tamura 1992 ) . A gamma distribution ( gamma parametric quantity value = 0.57 ) was used to pattern the rate of fluctuation amongst sites. Thirteen nucleotide sequences were used in the analysis of the UPGMA tree ( Tamura et al. 2011 ) . The concluding dataset had a sum of 16 311 places ( Tamura et al. 2011 ) .
Changes occur in the evolutionary procedure and that is derived by variable sites and the addition in these sites helps the evolutionary procedure even more ( Grievink et al. 2010 ) . The more variable sites in a proportionate mode the less accurate the building of the tree become ( Grievink et al. 2010 ) . So it will profit a coinage to increase the familial diverseness within because this in bend leads to better endurance due to a better fittingness ( Groom et al. 2006 ) . The per centum variable sites were moderately high at 35.2 % therefore demoing a good familial diverseness threw the alterations in the evolutionary procedure of Crocodylus. With a really high count of parsimoniousness enlightening sites of 2973 would demo how much alteration have occurred with in the genus and the outgroup ( Stearns and Hoekstra 2005 ) . Within the transition/transversion ratio the high figure indicated that passage was happening more frequent the transversion with in this information set. This high value shows a prejudice toward passage ( Tamura et al. 2007 ) . With the overall average distance at merely 10.8 % the persons are really near in the genus degree ( Kholodova 2009 ) . In the NJ tree ( Fig. 1 ) I could see that the length of the subdivisions in the concluding tree was plenty on the topology of the tree that was drawn ( Saitu and Nei 1987 ) . Crocodylus acutus, Crocodylus intermedius, Crocodylus johnsoni, Crocodylus mindorensis, Crocodylus moreletii, Crocodylus niloticus ( , Crocodylus novaeguineae, Crocodylus palustris, Crocodylus porosus, Crocodylus rhombifer and Crocodylus siamensis where all really near in relation to each other harmonizing to the HJ tree ( Saitu and Nei 1987 ) . The bootstrap value degree was set at 50 % so to give a position to within certain relationships with in Crocodylus ( Saitu and Nei 1987 ) . There is a spot of homoplasy though as the CI value is 0.653600 which is less than one ( Lipscomb 1998 ) . This phenomenon is a consequence of the similarity of certain independent events ( Stearns and Hoekstra 2005 ) . We need to cognize that with the addition in homoplasy the opportunities are at that place that we can free familial fluctuation and therefore loose speciation finally ( Groom et al. 2006, Stearns and Hoekstra 2005 ) .
` Therefore by looking at the informations signifier all the trees ( Fig. 1-4 ) we can clearly see the relationship to each of the Crocodylus genus persons. They are besides really near linages as to how they are geographically distributed which 1 could anticipate as the evolutionary line are close to each other. The statement of some earlier literature thought that crocodiles came from Africa seem to be incorrect. As two other independent surveies besides revealed the relation ship between the African, New World and Indopacific species ( Man et al. 2011, Meredith et Al. 2011, Oaks 2011 ) . I support these writers in stating that the crocodile coinage was ab initio from Australia therefore Indopacific ( Man et al. 2011, Meredith et Al. 2011, Oaks 2011 ) . The best decision is that C. niloticus crossed to the New World and this is why it is closer related to the New World species C. moreletii, C. intermedius, C. acutus and C. rhombifer ( Meredith et al. 2011 ) . So for preservation it is of import to observe that hybridisation can happen with in some of these species and that we need to take attention other wise we might free one or even both crossbreeding species and organize a new coinage at the disbursal of two species ( Meganathan et al. 2010 ) . So the right preservation schemes are important for the endurance of some of our endangered crocodiles. Thus it is my belief that by cognizing their familial relationship we can forestall this with proper genteelness plans.
In this survey we could besides touch on another large speaking point in the Crocodilian group as to where does G. gangeticus belong? Well I would back up most recent surveies which say they are a sister group of the genus Crocodylus but even more so they therefore need to be in their ain household ( Janke 2008, Li et Al. 2007 ) . In all the trees ( Fig. 1-4 ) we could see that, merely in the MP tree ( Fig. 3 ) it was different, that the G. gangeticus was under the O. tetraspis in all the trees which suggest that it is even further out of the household as O. tetraspis has already been classified as its ain household ( Janke 2008, Li et Al. 2007 ) . So we can now state with easiness that G. gangeticus is in its ain household with Tomistoma ( Janke 2008 ) .
BLAKE, D. & A ; LOVERIDGE, J. , 1975. The function of commercial crocodile farming in crocodile preservation. Biological preservation, ( September 1974 ) , pp.261-272.
DE VOS, A. , 1984. Crocodile preservation in India. Biological preservation, 29 ( May 1982 ) , pp.183-189.
FELSENSTEIN, J. 1985. Assurance bounds on evolutions: An attack utilizing the bootstrap. Evolution 39:783-791.
GRIEVINK, L.S. , PENNY, D. , HENDY, M.D. & A ; HOLLAND, B.R. 2010. Phylogenetic tree Reconstruction truth and theoretical account tantrum when proportions of variable sites change across the tree. Society of Systematic Biologists 59 ( 3 ) : 288-297.
GROOM, M.J. , MEFFE, G.K. & A ; CARROLL, C.R. 2006. Principles of Conservation Biology. 3rd erectile dysfunction. Sunderland, MA: Sinauer Associates.
IUCN 2012. IUCN Red List of Threatened Species. Version 2012.1. & lt ; hypertext transfer protocol: //www.iucnredlist.org & gt ; . Downloaded on 30 October 2012
JANKE, A. , 2008. The importance of informations trying and analysis in deciding crocodylian relationships. Molecular phylogenetics and development, 49, pp.407-8.
KHOLODOVA, M. V. 2009. Comparative phylogeography: molecular methods, ecological reading. Molecular Biology 43 ( 5 ) : 847-854.
LI, Y. , WU, X. , JI, X. , YAN, P. , AMATO, G. , 2007. The complete mitochondrial genome of salt-water crocodile ( Crocodylus porosus ) and evolution of crocodilians. Journal of genetic sciences and genomics = Yi chuan xue bao, 34, pp.119-28.
LIPSCOMB, D. 1998. Basicss of Cladistic Analysis. George Washington University, Washington D.C.
MAN, Z. , YISHU, W. , PENG, Y. , XIAOBING, W. , 2011. Crocodilian evolution inferred from 12 mitochondrial protein-coding cistrons, with new complete mitochondrial genomic sequences for Crocodylus acutus and Crocodylus novaeguineae. Molecular phylogenetics and development, 60, pp.62-7.
MAZZOTTI, F.J. , BEST, G.R. , BRANDT, L. , CHERKISS, M.S. , JEFFERY, B.M. , RICE, K.G. , 2009. Alligators and crocodiles as indexs for Restoration of Everglades ecosystems. Ecological Indexs, 9 ( 6 ) , pp.S137-S149.
MCGREGOR, J. , 2005. Crocodile offenses: people versus wildlife and the political relations of postcolonial preservation on Lake Kariba, Zimbabwe. Geoforum, 36 ( 3 ) , pp.353-369.
MEGANATHAN, P.R. , DUBEY, B. , BATZER, M. , RAY, D. , HAQUE, I. , 2010. Molecular phyletic analyses of genus Crocodylus ( Eusuchia, Crocodylia, Crocodylidae ) and the systematic place of Crocodylus porosus. Molecular phylogenetics and development, 57, pp.393-402.
MEREDITH, R.W. , HEKKALA, E.R. , AMATO, G. , GATESY, J. , 2011. A phyletic hypothesis for Crocodylus ( Crocodylia ) based on mitochondrial Deoxyribonucleic acid: grounds for a trans-Atlantic ocean trip from Africa to the New World. Molecular phylogenetics and development, 60, pp.183-91.
MILLS, L.S. , SOULE, M.E. & A ; DOAK, D.F. 1993. The Keystone-Species Concept in Ecology and preservation. Life science 43 ( 4 ) : 219-224.
MOLLES, M.C. 2010. Ecology: Concepts and Applications. 5th Ed. McGraw Hill Higher Education, New York.
NEI, M. & A ; KUMAR, S. 2000. Molecular Development and Phylogenetics. Oxford University Press, New York.
OAKS, J. , 2011. A Time-Calibrated Species Tree Of Crocodylia Reveals A Recent Radiation Of The True Crocodiles. Evolution, pp.3285-3297.
PAINE, R.T. 1995. A Conversation on Polishing the Concept of Keystone Species.Conservation Biology 9 ( 4 ) : 962-96
SAITOU, N. & A ; NEI, M. 1987. The neighbor-joining method: A new method for retracing phyletic trees. Molecular Biology and Evolution 4:406-425.
STEARNS, S.C. , & A ; HOEKSTRA, R.F. 2005. Development: An debut. 2nd Ed. Oxford University Press, New York.
TAMURA, K. 1992. Appraisal of the figure of nucleotide permutations when there are strong transition-transversion and G + C-content prejudices. Molecular Biology and Evolution 9:678-687.
TAMURA, K. & A ; NEI, M. 1993. Appraisal of the figure of nucleotide permutations in the control part of mitochondrial DNA in worlds and Pan troglodytess. Molecular Biology and Evolution 10:512-526.
TAMURA, K. , DUDLEY, J. , NEI, M. & A ; KUMAR, S. 2007. MEGA4: Molecular Evolutionary Genetics Analysis ( MEGA ) package version 4.0. Molecular Biology and Evolution 24:1596-1599.
TAMURA, K. , PETERSON, D. , PETERSON, N. , STECHER, G. , NEI, M. & A ; KUMAR, S. 2011. MEGA5: Molecular Evolutionary Genetics Analysis utilizing Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution ( In imperativeness ) .