Morphological identification of fish taxa can sometimes prove difficult because phenotypic variation is either being affected by environmental factors, phenotypic characters are highly conserved or marker selection has been inappropriate. DNA based markers especially neutral mitochondrial DNA (mtDNA) have been used widely in recent times to provide better resolution of systematic relationships among vertebrate taxa. The Asian Arowana (Scleropages formosus) is a high value ornamental fish belonging to the family Osteoglossidae with a number of different colour variants distributed geographically across different locations around Southeast Asia. Systematic relationships among colour variants still remain unresolved. Partial sequences of the Cytochrome B (Cyt B) and DNA barcoding gene, Cytochrome C Oxidase I (COI) were used here to assess genetic relationships among colour variants and as a tool for molecular identification for differentiating among colour variants in this species. Results of the study show that in general, colour pattern shows no relationship with extent of COI or Cyt B mtDNA differentiation and so cannot be used to identify taxa. Partial sequences of the mtDNA genes were sufficient however, to identify S. formosus from a closely related species within the order Osteoglossidae.
This paper describes the development and application of multilocus sequencing typing (MLST) and multi-virulence locus sequencing typing (MVLST) methods in determining the genetic variation and relatedness of 43 Vibrio cholerae strains of different serogroups isolated from various sources in Malaysia. The MLST assay used six housekeeping genes (dnaE, lap, recA, gyrB, cat and gmd), while the MVLST assay incorporated three virulence genes (ctxAB, tcpA and tcpI) and three virulence-associated genes (hlyA, toxR and rtxA). Our data showed that the dnaE and rtxA genes were the most conserved genes in V. cholerae O1 strains. Among the 12 studied genes, transitional substitutions that led to silent mutations were observed in all, except for gmd and hlyA, while non-synonymous substitutions occurred more frequently in virulence and virulence-associated genes. Five V. cholerae O1 strains were found to be the El Tor variant O1 strains because they harboured the classical ctxB gene. In addition, the classical ctxB gene was also observed in O139 V. cholerae. A total of 29 MLST types were observed, and this assay could differentiate V. cholerae within the non-O1/non-O139 serogroups. A total of 27 MVLST types were obtained. MVLST appeared to be more discriminatory than MLST because it could differentiate V. cholerae strains from two different outbreaks and could separate the toxigenic from the non-toxigenic subtypes. Although the O1 V. cholerae strains were closely related, the combined MLST and MVLST analyses differentiated the strains isolated from different localities. In conclusion, sequence-based analysis in this study provided a better understanding of mutation points and the type of mutations in V. cholerae. The MVLST assay is useful to characterise O1 V. cholerae strains, while combined analysis may improve the discriminatory power and is suitable for the local epidemiological study of V. cholerae.
Haruan (Channa striatus) is in great demand in the Malaysian domestic fish market. In the present study, mtDNA cyt b was used to investigate genetic variation of C. striatus among populations in Peninsular Malaysia. The overall population of C. striatus demonstrated a high level of haplotype diversity (h) and a low-to-moderate level of nucleotide diversity (π). Analysis of molecular variance (AMOVA) results showed a significantly different genetic differentiation among 6 populations (F(ST)=0.37566, P=0.01). Gene flow (Nm) was high and ranged from 0.32469 to infinity (∞). No significant relationship between genetic distance and geographic distance was detected. A UPGMA tree based on the distance matrix of net interpopulation nucleotide divergence (d(A)) and haplotype network of mtDNA cyt b revealed that C. striatus is divided into 2 major clades. The neutrality and mismatch distribution tests for all populations suggested that C. striatus in the study areas had undergone population expansion. The estimated time of population expansion in the mtDNA cyt b of C. striatus populations occurred 0.72-6.19 million years ago. Genetic diversity of mtDNA cyt b and population structure among Haruan populations in Peninsular Malaysia will be useful in fisheries management for standardization for Good Agriculture Practices (GAP) in fish-farming technology, as well as providing the basis for Good Manufacturing Practices (GMP).
Data from three published studies of genetic variation at 18 microsatellite loci in water buffalo populations in China (18 swamp type, two river type), Nepal (one wild, one domestic river, one hybrid) and south-east Asia (eight swamp, three river) were combined so as to gain a broader understanding of genetic relationships among the populations and their demographic history. Mean numbers of alleles and expected heterozygosities were significantly different among populations. Estimates of θ (a measure of population differentiation) were significant among the swamp populations for all loci and among the river populations for most loci. Differentiation among the Chinese swamp populations (which was due primarily to just one population) was much less than among the south-east Asian. The Nepal wild animals, phenotypically swamp type but genetically like river type, are significantly different from all the domestic river populations and presumably represent the ancestral Bubalus arnee (possibly with some river-type introgression). Relationships among the swamp populations (D(A) genetic distances, principal component analysis and structure analyses) show the south-east Asian populations separated into two groups by the Chinese populations. Given these relationships and the patterns of genetic variability, we postulate that the swamp buffalo was domesticated in the region of the far south of China, northern Thailand and Indochina. Following domestication, it spread south through peninsular Malaysia to Sumatra, Java and Sulawesi, and north through China, and then to Taiwan, the Philippines and Borneo.
Tomato mosaic virus (ToMV), a member of the genus Tobamovirus, infects several ornamental and horticultural crops worldwide. In this study, the nucleotide sequences of the coat protein gene of worldwide ToMV isolates were analyzed to estimate the genetic structure and diversity of this virus and the involved evolutionary forces. The phylogenetic analysis showed three clades with high bootstrap support: Clade I contained three ToMV isolates from Brazil collected from pepper, Clade II comprised one Brazilian ToMV isolate from pepper, and Clade III was composed of ToMV isolates collected from different plant hosts (pepper, tomato, eggplant, lilac, camellia, dogwood, red spruce, etc.) and water (from melting ice, lakes and streams) from different countries: USA, Brazil, Korea, Germany, Spain, Denmark (Greenland), China, Taiwan, Malaysia, Iran, and Kazakhstan. With the exception of Brazil, nucleotide diversity within and between different geographic regions was very low, although statistical analyses suggested some gene flow between most of these regions. Our analyses also suggested a strong negative selection which could have contributed to the genetic stability of ToMV.
Gestational diabetes mellitus (GDM), defined as carbohydrate intolerance diagnosed for the first time during pregnancy, affects both maternal and fetal health. Possession of a specific genetic polymorphism can be a predisposing factor for susceptibility to some diseases. The aim of this study was to investigate the association between single nucleotide polymorphisms (SNP) in the promoter gene of interleukin-10 (IL-10) as well as tumor necrosis factor-alpha (TNF alpha) with the development of GDM. Two hundred and twelve consecutive series of eligible normal pregnant women (controls) and gestational diabetes mellitus women were selected based on the study's inclusion and exclusion criteria. DNA was extracted from blood and genotyped for IL-10 at three positions and TNF alpha for gene polymorphism using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Plasma levels of IL-10 and TNF alpha at different gestational periods as well as postpartum were quantified using enzyme linked immunosorbent assay (ELISA). The results of the study showed that the difference in the frequency of SNP at position -597 in the promoter of the human IL-10 gene between the control and GDM groups was statistically significant (p < 0.05). In contrast, there was no significant difference in the frequency of SNP at the other two sites in the promoter region of the human IL-10 gene (-824 and -1082) as well as position -308 in the promoter of the human TNF-alpha (p > 0.05). In addition, there was no significant difference between the two groups in terms of plasma levels of IL-10 as well as TNF alpha in different stages of pregnancy. SNP at position -597 was significantly associated with the development of GDM and shows potential for use as a predictive marker for GDM.
One hundred isolates of Phytophthora infestans collected from 10 provinces in China between 1998 and 2004 were analyzed for mating type, metalaxyl resistance, mitochondrial DNA (mtDNA) haplotype, allozyme genotype, and restriction fragment length polymorphism (RFLP) with the RG-57 probe. In addition, herbarium samples collected in China, Russia, Australia, and other Asian countries were also typed for mtDNA haplotype. The Ia haplotype was found during the first outbreaks of the disease in China (1938 and 1940), Japan (1901, 1930, and 1931), India (1913), Peninsular Malaysia (1950), Nepal (1954), The Philippines (1910), Australia (1917), Russia (1917), and Latvia (1935). In contrast, the Ib haplotype was found after 1950 in China on both potato and tomato (1952, 1954, 1956, and 1982) and in India (1968 and 1974). Another migration of a genotype found in Siberia called SIB-1 (Glucose-6-phosphate isomerase [Gpi] 100/100, Peptidase [Pep] 100/100, IIa mtDNA haplotype) was identified using RFLP fingerprints among 72% of the isolates and was widely distributed in the north and south of China and has also been reported in Japan. A new genotype named CN-11 (Gpi 100/111, Pep 100/100, IIb mtDNA haplotype), found only in the south of China, and two additional genotypes (Gpi 100/100, Pep 100/100, Ia mtDNA haplotype) named CN-9 and CN-10 were identified. There were more diverse genotypes among isolates from Yunnan province than elsewhere. The SIB-1 (IIa) genotype is identical to those from Siberia, suggesting later migration of this genotype from either Russia or Japan into China. The widespread predominance of SIB-1 suggests that this genotype has enhanced fitness compared with other genotypes found. Movement of the pathogen into China via infected seed from several sources most likely accounts for the distribution of pathogen genotypes observed. MtDNA haplotype evidence and RFLP data suggest multiple migrations of the pathogen into China after the initial introduction of the Ia haplotype in the 1930s.
Central Asia and the Indian subcontinent represent an area considered as a source and a reservoir for human genetic diversity, with many markers taking root here, most of which are the ancestral state of eastern and western haplogroups, while others are local. Between these two regions, Terai (Nepal) is a pivotal passageway allowing, in different times, multiple population interactions, although because of its highly malarial environment, it was scarcely inhabited until a few decades ago, when malaria was eradicated. One of the oldest and the largest indigenous people of Terai is represented by the malaria resistant Tharus, whose gene pool could still retain traces of ancient complex interactions. Until now, however, investigations on their genetic structure have been scarce mainly identifying East Asian signatures.
Genetic variation due to heavy metal contamination has always been an interesting topic of study. Because of the numerous contaminants being found in coastal and intertidal waters, there is always much discussion and argument as to which contaminant(s) caused the variations in the genetic structures of biomonitors. This study used a Single Primer Amplification Reaction (SPAR) technique namely Random Amplified Polymorphic DNA (RAPD) to determine the genetic diversity of the populations of the green-lipped mussel Perna viridis collected from a metal-contaminated site at Kg. Pasir Puteh and those from four relatively' uncontaminated sites (reference sites). Heavy metal levels (Cd, Cu, Pb and Zn) were also measured in the soft tissues and byssus of the mussels from all the sites. Cluster analyses employing UPGMA done based on the RAPD makers grouped the populations into two major clusters; the Bagan Tiang, Pantai Lido, Pontian and Kg. Pasir Puteh populations were in one cluster, while the Sg. Belungkor population clustered by itself. This indicated that the genetic diversity based on bands resulting from the use of all four RAPD primers on P. viridis did not indicate its potential use as a biomarker of heavy metal pollution in coastal waters. However, based on a correlation analysis between a particular metal and a band resulting from a specific RAPD primer revealed some significant (P < 0.01) correlations between the primers and the heavy metal concentrations in the byssus and soft tissues. Thus, the correlation between a particular metal and the bands resulting from the use of a specific RAPD primer on P. viridis could be used as biomonitoring tool of heavy metal pollution.
We examined allozyme variation in two camaenid tree snails, Amphidromus atricallosus and A. inversus, across two principal regions of Thailand and from Singapore, plus for A. inversus, one site in peninsular Malaysia. Using horizontal starch gel electrophoresis, 13 allozyme loci (11 polymorphic) were screened for A. atricallosus and 18 (5 polymorphic) for A. inversus. Heterozygosity was higher in A. atricallosus (Hexp=0.018-0.201, mean=0.085) than in A. inversus (Hexp=0-0.023, mean= 0.002). Genetic heterogeneity among samples was higher in A. inversus (Fst=0.965) than in A. atricallosus (Fst=0.781). Within A. atricallosus, populations were more differentiated in southern Thailand (Fst=0.551) than in eastern Thailand (Fst=0.144). The high Fst and low Hexp in populations of A. inversus suggest that this species is likely to have experienced a series of strong bottlenecks, perhaps occurring chiefly on offshore continental-shelf islands. The low Fst values of A. atricallosus in eastern Thailand suggest frequent gene flows among populations in this region. The southern and eastern samples of A. atricallosus exhibited fixed allele differences at four loci and great genetic distance (Nei's D=0.485-0.946), suggesting that these two samples may actually represent, or else be evolving into, separate species.
An 835 base pair (bp) fragment of mitochondrial DNA (mtDNA) was sequenced to characterize genetic variation within and among 1,053 samples comprising five regional populations each of longtail macaques (Macaca fascicularis) and rhesus macaques (Macaca mulatta), and one sample each of Japanese (M. fuscata) and Taiwanese (M. cyclopis) macaques. The mtDNA haplotypes of longtail macaques clustered in two large highly structured clades (Fas1 and Fas2) of a neighbor-joining tree that were reciprocally monophyletic with respect to those representing rhesus macaques, Japanese macaques, and Taiwanese macaques. Both clades exhibited haplotypes of Indonesian and Malaysian longtail macaques widely dispersed throughout them; however, longtail macaques from Indochina, Philippines, and Mauritius each clustered in a separate well-defined clade together with one or a few Malaysian and/or Indonesian longtail macaques, suggesting origins on the Sunda shelf. Longtail macaques from Malaysia and Indonesia were far more genetically diverse, and those from Mauritius were far less diverse than any other population studied. Nucleotide diversity between mtDNA sequences of longtail macaques from different geographic regions is, in some cases, greater than that between Indian and Chinese rhesus macaques. Approximately equal amounts of genetic diversity are due to differences among animals in the same regional population, different regional populations, and different species. A greater proportion of genetic variance was explained by interspecies differences when Japanese and Taiwanese macaques were regarded as regional populations of rhesus macaques than when they were treated as separate species. Rhesus macaques from China were more closely related to both Taiwanese and Japanese macaques than to their own conspecifics from India.
Invasive species are one of the main sources of the ongoing global loss of biodiversity. Invasive ants are known as particularly damaging invaders and their introductions are often accompanied by population-level behavioural and genetic changes that may contribute to their success. Anoplolepis gracilipes is an invasive ant that has just recently received increased attention due to its negative impact on native ecosystems. We examined the behaviour and population structure of A. gracilipes in Sabah, Malaysia. A total of 475 individuals from 24 colonies were genotyped with eight microsatellite markers. Intracolonial relatedness was high, ranging from 0.37 to 1 (mean +/- SD: 0.82 +/- 0.04), while intercolonial relatedness was low (0.0 +/- 0.02, range -0.5-0.76). We compared five distinct sampling regions in Sabah and Brunei. A three-level hierarchical F-analysis revealed high genetic differentiation among colonies within the same region, but low genetic differentiation within colonies or across regions. Overall levels of heterozygosity were unusually high (mean H(O) = 0.95, mean H(E) = 0.71) with two loci being entirely heterozygous, indicating an unusual reproductive system in this species. Bioassays revealed a negative correlation between relatedness and aggression, suggesting kinship as one factor facilitating supercolony formation in this species. Furthermore, we genotyped one individual per nest from Sabah (22 nests), Sarawak (one nest), Brunei (three nests) and the Philippines (two nests) using two mitochondrial DNA markers. We found six haplotypes, two of which included 82.1% of all sequences. Our study shows that the sampled area in Sabah consists of a mosaic of differently interrelated nests in different stages of colony establishment. While some of the sampled colonies may belong to large supercolonies, others are more likely to represent recently introduced or dispersed propagules that are just beginning to expand.
The species diversity and genetic structure of mosquitoes belonging to the Anopheles maculatus group in Southeast Asia were investigated using the internal transcribed spacer 2 (ITS2) of ribosomal DNA (rDNA). A molecular phylogeny indicates the presence of at least one hitherto unrecognised species. Mosquitoes of chromosomal form K from eastern Thailand have a unique ITS2 sequence that is 3.7% divergent from the next most closely related taxon (An. sawadwongporni) in the group. In the context of negligible intraspecific variation at ITS2, this suggests that chromosomal form K is most probably a distinct species. Although An. maculatus sensu stricto from northern Thailand and southern Thailand/peninsular Malaysia differ from each other in chromosomal banding pattern and vectorial capacity, no intraspecific variation was observed in the ITS2 sequences of this species over this entire geographic area despite an extensive survey. A PCR-based identification method was developed to distinguish five species of the group (An. maculatus, An. dravidicus, An. pseudowillmori, An. sawadwongporni and chromosomal form K) to assist field-based studies in northwestern Thailand. Sequences from 187 mosquitoes (mostly An. maculatus and An. sawadwongporni) revealed no intraspecific variation in specimens from Thailand, Cambodia, mainland China, Malaysia, Taiwan and Vietnam, suggesting that this identification method will be widely applicable in Southeast Asia. The lack of detectable genetic structure also suggests that populations of these species are either connected by gene flow and/or share a recent common history.
We investigated the genetic structure within and among Bornean orang-utans (Pongo pygmaeus) in forest fragments of the Lower Kinabatangan flood plain in Sabah, Malaysia. DNA was extracted from hair and faecal samples for 200 wild individuals collected during boat surveys on the Kinabatangan River. Fourteen microsatellite loci were used to characterize patterns of genetic diversity. We found that genetic diversity was high in the set of samples (mean H(E) = 0.74) and that genetic differentiation was significant between the samples (average F(ST) = 0.04, P < 0.001) with F(ST) values ranging from low (0.01) to moderately large (0.12) values. Pairwise F(ST) values were significantly higher across the Kinabatangan River than between samples from the same river side, thereby confirming the role of the river as a natural barrier to gene flow. The correlation between genetic and geographical distance was tested by means of a series of Mantel tests based on different measures of geographical distance. We used a Bayesian method to estimate immigration rates. The results indicate that migration is unlikely across the river but cannot be completely ruled out because of the limited F(ST) values. Assignment tests confirm the overall picture that gene flow is limited across the river. We found that migration between samples from the same side of the river had a high probability indicating that orang-utans used to move relatively freely between neighbouring areas. This strongly suggests that there is a need to maintain migration between isolated forest fragments. This could be done by restoring forest corridors alongside the river banks and between patches.
Molecular variations of Spiranthes sinensis Ames var. australis (R.Br.) H. Hara et Kitam. ex Kitam. in Japan were examined to evaluate the validity of the seasonally differentiated groups and a dwarf form of the species, which is endemic to Yakushima Island, Japan. Sequence differences in the plastid trnL-F locus clearly distinguished Japanese S. sinensis var. australis from S. sinensis var. sinensis collected from Ryukyu. In contrast, the trnL-F sequence of S. sinensis var. australis from Sabah, Malaysia, clearly differed from that of Japanese S. sinensis var. australis, suggesting genetic heterogeneity of Spiranthes sinensis var. australis in Asia. Moreover, a molecular analysis based on the sequences of nuclear ITS1 regions indicated that there are two major groups of S. sinensis var. australis in Japan, with a geographic distribution boundary on Kyushu Island. However, the trnL-F and ITS1 sequences did not support the genetic differentiation of the seasonally differentiated groups or the dwarf form from the other Japanese individuals. Based on these molecular data, the systematic treatment of physiological and morphological variations in the Japanese population of S. sinensis. var. australis is discussed.
The combination of two silent mutations, c.1311C>T in exon 11 and IVS11 T93C (glucose-6-phosphate dehydrogenase (G6PD) 1311T/93C), with unknown mechanism, have been reported in G6PD-deficient individuals in Asian populations including Malaysian aboriginal group, Negrito. Here, we report the screening of G6PD gene in 103 Negrito volunteers using denaturing high-performance liquid chromatography (dHPLC) and direct sequencing. A total of 48 individuals (46.6%) were G6PD deficient, 83.3% of these carried G6PD 1311T/93C with enzyme activity ranging from 1.8 to 4.8 U gHb(-1). Three novel single-nucleotide polymorphisms (SNPs), rs112950723, rs111485003 and rs1050757, were found in the G6PD 3'-untranslated region (UTR). Strong association was observed between haplotype 1311T/93C and rs1050757G, which is located inside the 35 bp AG-rich region. In silico analysis revealed that the transition of A to G at position rs1050757 makes significant changes in the G6PD mRNA secondary structure. Moreover, putative micro (mi)RNA target sites were identified in 3'-UTR of G6PD gene, two of these in the region encompassing rs1050757. It could be speculated that rs1050757 have a potential functional effect on the downregulation of mRNA and consequently G6PD deficiency either by affecting mRNA stability and translation or mirRNA regulation process. This is the first report of biochemical association of an SNP in 3'-UTR of G6PD gene and the possible role of mRNA secondary structure.
Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30× coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (<5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.
Giardia duodenalis is a protozoan parasite that can cause significant diarrhoeal diseases. Knowledge of population genetics is a prerequisite for ascertaining the invasion patterns of this parasite. In order to infer evolutionary patterns that could not be uncovered based on the morphological features, a population genetic study with the incorporation of molecular marker was carried out to access the genetic structure of G. duodenalis isolated from the Malaysian population and the global populations.
Angiostrongylus cantonensis is a zoonotic parasite that causes eosinophilic meningitis in humans. Earlier work on its mitochondrial genome was based on long polymerase chain reaction method. To date, only the mitogenome of the isolates from China has been studied. We report here the complete mitogenome of the Thailand isolate based on next generation sequencing and compare the genetic diversity with other isolates. The mitogenome of the Thailand isolate (13,519bp) is longer than those of the China isolates (13,497-13,502bp). Five protein-coding genes (atp6, cox1, cox2, cob, nad2) show variations in length among the isolates. The stop codon of the Thailand isolate differs from the China and Taiwan isolates in 4 genes (atp6, cob, nad2, nad6). Additionally, the Thailand isolate has 4 incomplete T stop codon compared to 3 in the China and Taiwan isolates. The control region is longer in the Thailand isolate (258bp) than the China (230-236bp) and Taiwan (237bp) isolates. The intergenic sequence between nad4 and cox1 genes in the Thailand isolate lacks 2bp (indels) at the 5'-end of the sequence as well as differs at 7 other sites compared to the China and Taiwan isolates. In the Thailand isolate, 18 tRNAs lack the entire TΨC-arm, compared to 17 in the China isolate and 16 in the Taiwan isolate. Phylogenetic analyses based on 36 mt-genes, 12 PCGs, 2 rRNA genes, 22 tRNA genes and control region all indicate closer genetic affinity between the China and Taiwan isolates compared to the Thailand isolate. Based on 36 mt-genes, the inter-isolate genetic distance varies from p=3.2% between China and Taiwan isolates to p=11.6% between Thailand and China isolates. The mitogenome will be useful for population, phylogenetics and phylogeography studies.
Simulium feuerborni Edwards is geographically widespread in Southeast Asia. Previous cytogenetic study in Thailand revealed that this species is a species complex composed of two cytoforms (A and B). In this study, we cytologically examined specimens obtained from the Cameron Highlands, Malaysia, and Puncak, Java, Indonesia. The results revealed two additional cytoforms (C and D) of S. feuerborni. Specimens from Malaysia represent cytoform C, differentiated from other cytoforms by a fixed chromosome inversion on the long arm of chromosome III (IIIL-5). High frequencies of the B chromosome (33-83%) were also observed in this cytoform. Specimens from Indonesia represent the cytoform D. This cytoform is differentiated from others by a fixed chromosome inversion difference on the long arm of chromosome II (IIL-4). Mitochondrial DNA sequences support genetic differentiation among cytoforms A, B, and C. The pairwise F(ST) values among these cytoforms were highly significantly consistent with the divergent lineages of the cytoforms in a median-joining haplotype network. However, a lack of the sympatric populations prevented us from testing the species status of the cytoforms.