Establishment of a mitochondrial DNA sequence database for the identification of fish species commercially available in South Africa

The limitations intrinsic to morphology-based identification systems have created an urgent need for reliable genetic methods that enable the unequivocal recognition of fish species, particularly those that are prone to overexploitation and/or market substitution. The aim of this study was to develop a comprehensive reference library of DNA sequence data to allow the explicit identification of 53 commercially available fish species in South Africa, most of which were locally caught marine species. Sequences of approximately 655 base pairs were generated for all species from the cytochrome c oxidase I (COI) gene, the region widely adopted for DNA barcoding. Specimens of the genus Thunnus were examined in further detail, employing additional mitochondrial DNA control region sequencing. Cumulative analysis of the sequences from the COI region revealed mean conspecific, congeneric and confamilial Kimura 2-parameter distances of 0.10%, 4.58% and 15.43%, respectively. The results showed that the vast majority (98%) of fish species examined could be readily differentiated by their COI barcodes, but that supplementary control region sequencing was more useful for the discrimination of three Thunnus species. Additionally, the analysis of COI data raised the prospect that Thyrsites atun (snoek) could constitute a species pair. The present study has established the necessary genetic information to permit the unambiguous identification of 53 commonly marketed fish species in South Africa, the applications of which hold a plethora of benefits relating to ecology research, fisheries management and control of commercial practices.     

Species identification of small fish in Xixuan Island coastal waters of Zhoushan using DNA barcoding

This study aimed to assess the fish species composition in Xixuan Island coastal waters and confirm the applicability of DNA barcoding in fish identification. Fish samples were caught with a gill net in December 2016 and in January, February, April, May and June 2017, and the fish species in these samples were analyzed and identified via both traditional morphological methods and DNA barcoding (mitochondrial 12S rRNA). The sample contained 1,357 individuals, the majority of which were Tridentiger barbatus (829 individuals). A total of 154 individuals were selected for DNA barcode analysis, and 147 sequences were obtained. A total of 35 species, which belonged to six orders and 16 families, were accurately identified. One species was identified to the genus level (Repomucenus), and only one was identified to the family level (Moridae). Most of the species belonged to Gobiidae (12 species). With the exception of those in Moridae, the identified species were commonly detected on the Zhejiang coast. The average Kimura two-parameter (K2P) distances within species, genera, families, and order and between orders were 0.15%, 4.44%, 14.96%, 20.37%, and 26.20%, respectively. The intraspecific K2P distance was markedly greater than the interspecific distance. The A + T content (55.41%) was higher than the G + C content (44.58%). Phylogenetic tree analysis showed that individuals belonging to the same species were clustered together and could be clearly distinguished. Mitochondrial 12S rRNA can be effectively used for juvenile and adult fish identification and has good application prospects in fish species identification. The investigation and species identification of small fish collected monthly from Xixuan Island coastal waters are conducive to understanding the distribution characteristics and species composition of fish in the Zhoushan coastal area.     

Molecular identification of Korean catfish (Siluriformes) based on two genetic markers

Eleven species of catfish (Siluriformes) found in Korean freshwater are economically important resources and include five endemic species. However, there are no studies on phylogenetic analysis of all catfish species in Korea at a molecular level. The species-level analysis of catfish species is usually carried out through morphological characters and controversial due to phenotypic variation. In this study, the partial sequences of 16S rRNA and cytochrome b mitochondrial genes were analyzed for species identification and phylogenetic relationships among 11 species of catfish from 10 different rivers in Korea. The nucleotide sequences of 16S rRNA and cytochrome b consisted of 587 and 441 nucleotide base pairs, respectively. Sequence analysis of both genes revealed that the 11 species fell into three distinct groups, which were genetically distinct from each other and exhibited identical phylogenetic resolution. Sequence divergences between congeneric species averaged 1.78% and 7.39% for 16S rRNA and cytochrome b, respectively. The phylogenetic relationships forming well-differentiated clades in the NJ, ML and BI trees were identical for both fragments. This research demonstrates that partial sequences of both the genes can efficiently identify the 11 species of catfish in Korea, indicating the usefulness of mtDNA-based approach in species identification.     

Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians

BACKGROUND: Identifying species of organisms by short sequences of DNA has been in the center of ongoing discussions under the terms DNA barcoding or DNA taxonomy. A C-terminal fragment of the mitochondrial gene for cytochrome oxidase subunit I (COI) has been proposed as universal marker for this purpose among animals. RESULTS: Herein we present experimental evidence that the mitochondrial 16S rRNA gene fulfills the requirements for a universal DNA barcoding marker in amphibians. In terms of universality of priming sites and identification of major vertebrate clades the studied 16S fragment is superior to COI. Amplification success was 100% for 16S in a subset of fresh and well-preserved samples of Madagascan frogs, while various combination of COI primers had lower success rates.COI priming sites showed high variability among amphibians both at the level of groups and closely related species, whereas 16S priming sites were highly conserved among vertebrates. Interspecific pairwise 16S divergences in a test group of Madagascan frogs were at a level suitable for assignment of larval stages to species (1-17%), with low degrees of pairwise haplotype divergence within populations (0-1%). CONCLUSION: We strongly advocate the use of 16S rRNA as standard DNA barcoding marker for vertebrates to complement COI, especially if samples a priori could belong to various phylogenetically distant taxa and false negatives would constitute a major problem.     

大黄鱼16SrRNA和18SrRNA基因的测定与序列分析

利用多对引物,扩增并测定出大黄鱼16SrRNA基因和18SrRNA基因的部分序列,其长度分别为1202bp和1275bp,16SrRNA基因序列的GC含量为46．12％,18SrRNA基因的Gc含量为53．oo％。将大黄鱼16SrRNA基因序列与GenBank中15种硬骨鱼类的同源序列结合,同时将其18SrRNA基因序列与GenBank中9种脊索动物的同源序列相结合,运用软件获得各自序列间差异百分比,转换和颠换数值等信息。基于这两种基因序列,利用NJ法和BI法,分别构建16种硬骨鱼类和10种脊索动物的分子系统树。18SrRNA构建的系统树包括三大支,一支为哺乳类、鸟类和爬行类共6个物种,一支为两栖类的1个物种,另一支为2种硬骨鱼类。16SrRNA构建的系统树显示大黄鱼所在的石首鱼科与鲈科和盖刺鱼科亲缘关系较近。此外还讨论了这两个基因的序列特征。     

Sequence arrangement of the 16S and 26S rRNA genes in the pathogenic haemoflagellate Leishmania donovani.

Kinetic and chemical analysis show that the haploid genome of Leishmania donovani has between 4.6 and 6.5 X 10(7) Kb pairs of DNA. Cot analysis shows that the genome contains 12% rapidly reassociating DNA, U3% middle repetitive DNA with an average reiteration frequency of 77 and 62% single copy DNA. Saturation hybridization experiments show that 0.82% of the nuclear DNA is occupied by rRNA coding sequences. The average repetition frequency of these sequences is determined to be 166. Sedimentation velocity studies indicate the two major rRNA species have sedimentation values of 26S and 16S, respectively. The arrangement of the rRNA genes and their spacer sequences on long strands of purified rDNA has been determined by the examination of the structure of rRNA:DNA hybrids prepared for electron microscopy by the gene 32-ethidium bromide technique. Long DNA strands are observed to contain several gene sets (16S + 26S). One repeat unit contains the following sequences in the order given: (a) A 16S gene of length 2.12 Kb, (b) An internal transcribed spacer (Spl) of length 1.23 Kb, which contains a short sequence that may code for a 5.8S rRNA, (C) 26S gene with a length of 4.31 Kb which contains an internal gap region of length 0.581 Ib, (d) An external spacer of average length 5.85 Kb.     

Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes.

PCR is routinely used in amplification and cloning of rRNA genes from environmental DNA samples for studies of microbial community structure and identification of novel organisms. There have been concerns about generation of chimeric sequences as a consequence of PCR coamplification of highly conserved genes, because such sequences may lead to reports of nonexistent organisms. To quantify the frequency of chimeric molecule formation, mixed genomic DNAs from eight actinomycete species whose 16S rRNA sequences had been determined were used for PCR coamplification of 16S rRNA genes. A large number of cloned 16S ribosomal DNAs were examined by sequence analysis, and chimeric molecules were identified by multiple-sequence alignment with reference species. Here, we report that the level of occurrence of chimeric sequences after 30 cycles of PCR amplification was 32%. We also show that PCR-induced chimeras were formed between different rRNA gene copies from the same organism. Because of the wide use of PCR for direct isolation of 16S rRNA sequences from environmental DNA to assess microbial diversity, the extent of chimeric molecule formation deserves serious attention.     

DNA Microarrays for Identifying Fishes

In many cases marine organisms and especially their diverse developmental stages are difficult to identify by morphological characters. DNA-based identification methods offer an analytically powerful addition or even an alternative. In this study, a DNA microarray has been developed to be able to investigate its potential as a tool for the identification of fish species from European seas based on mitochondrial 16S rDNA sequences. Eleven commercially important fish species were selected for a first prototype. Oligonucleotide probes were designed based on the 16S rDNA sequences obtained from 230 individuals of 27 fish species. In addition, more than 1200 sequences of 380 species served as sequence background against which the specificity of the probes was tested in silico. Single target hybridisations with Cy5-labelled, PCR-amplified 16S rDNA fragments from each of the 11 species on microarrays containing the complete set of probes confirmed their suitability. True-positive, fluorescence signals obtained were at least one order of magnitude stronger than false-positive cross-hybridisations. Single nontarget hybridisations resulted in cross-hybridisation signals at approximately 27% of the cases tested, but all of them were at least one order of magnitude lower than true-positive signals. This study demonstrates that the 16S rDNA gene is suitable for designing oligonucleotide probes, which can be used to differentiate 11 fish species. These data are a solid basis for the second step to create a "Fish Chip" for approximately 50 fish species relevant in marine environmental and fisheries research, as well as control of fisheries products.     

Yersinia spp. Identification Using Copy Diversity in the Chromosomal 16S rRNA Gene Sequence

API 20E strip test, the standard for Enterobacteriaceae identification, is not sufficient to discriminate some Yersinia species for some unstable biochemical reactions and the same biochemical profile presented in some species, e.g. Yersinia ferderiksenii and Yersinia intermedia, which need a variety of molecular biology methods as auxiliaries for identification. The 16S rRNA gene is considered a valuable tool for assigning bacterial strains to species. However, the resolution of the 16S rRNA gene may be insufficient for discrimination because of the high similarity of sequences between some species and heterogeneity within copies at the intra-genomic level. In this study, for each strain we randomly selected five 16S rRNA gene clones from 768 Yersinia strains, and collected 3,840 sequences of the 16S rRNA gene from 10 species, which were divided into 439 patterns. The similarity among the five clones of 16S rRNA gene is over 99% for most strains. Identical sequences were found in strains of different species. A phylogenetic tree was constructed using the five 16S rRNA gene sequences for each strain where the phylogenetic classifications are consistent with biochemical tests; and species that are difficult to identify by biochemical phenotype can be differentiated. Most Yersinia strains form distinct groups within each species. However Yersinia kristensenii, a heterogeneous species, clusters with some Yersinia enterocolitica and Yersinia ferderiksenii/intermedia strains, while not affecting the overall efficiency of this species classification. In conclusion, through analysis derived from integrated information from multiple 16S rRNA gene sequences, the discrimination ability of Yersinia species is improved using our method.     

Extreme intraspecific mitochondrial DNA sequence divergence in Galaxias maculatus (Osteichthys: Galaxiidae), one of the world's most widespread freshwater fish

Biogeographic controversies surrounding the widespread freshwater fish, Galaxias maculatus, were addressed with DNA sequence data. Mitochondrial cytochrome b and 16S rRNA sequences were obtained from representatives of six populations of this species. Substantial levels of cytochrome b (maximum 14.6%) and 16S rRNA sequence divergence (maximum 6.0%) were detected between western Pacific (Tasmania-New Zealand) and South American (Chile-Falkland Islands) haplotypes. A considerable level of divergence was also detected between Tasmanian and New Zealand haplotypes (maximum 5.1%) and within and among Chilean and Falkland Island G. maculatus (maximum 3. 8%). The phylogenetic structure of haplotypes conflicts with the accepted pattern of continental fragmentation. Molecular clock calibrations suggest that haplotype divergences postdate the fragmentation of Gondwana. These findings point to marine dispersal rather than ancient vicariance as an explanation for the wide distribution. The phylogenetic structure of South American haplotypes was not consistent with their geographic distribution. We consider factors such as population divergence, population size, dispersal, secondary contact, and philopatry as potential causes of the high level of mtDNA nucleotide diversity in this species.     

Routine DNA analysis based on 12S rRNA gene sequencing as a tool in the management of captive primates

Automated DNA sequencing of a fragment of the relatively slowly evolving mitochondrial 12S rRNA gene was used to distinguish primate species, and the method was compared with species determination based upon classical taxonomy. DNA from blood from 53 monkeys housed at the Stichting AAP Shelter for Exotic Animals, all Old World monkeys, was amplified by polymerase chain reaction (PCR) with a primer set spanning approximately 390 nucleotides of the mitochondrial 12S rRNA gene. The products were directly sequenced and compared with our database of primate 12S sequences. Many individuals were found to harbor a 12S sequence identical to one of the reference sequences. For others, phylogenetic methods were used for species estimation, which was especially informative in Cercopithecus species.     

Comparison between the complete mtDNA sequences of the blue and the fin whale,two species that can hybridize in nature

The sequence of the mitochondrial DNA (mtDNA) molecule of the blue whale (Balaenoptera musculus) was determined. The molecule is 16,402 by long and its organization conforms with that of other eutherian mammals. The molecule was compared with the mtDNA of the congeneric fin whale (B. physalus). It was recently documented that the two species can hybridize and that male offspring are infertile whereas female offspring may be fertile. The present comparison made it possible to determine the degree of mtDNA difference that occurs between two species that are not completely separated by hybridization incompatibility. The difference between the complete mtDNA sequences was 7.4%. Lengths of peptide coding genes were the same in both species. Except for a small portion of the control region, disruption in alignment was usually limited to insertion/deletion of a single nucleotide. Nucleotide differences between peptide coding genes ranged from 7.1 to 10.5%, and difference at the inferred amino acid level was 0.0–7.9%. In the rRNA genes the mean transition difference was 3.8%. This figure is similar in degree to the difference (3.4%) between the 12S rRNA gene of humans and the chimpanzee. The mtDNA differences between the two whale species, involving both peptide coding and rRNA genes, suggest an evolutionary separation of 5 million years. Although hybridization between more distantly related mammalian species may not be excluded, it is probable that the blue and fin whales are nearly as different in their mtDNA sequences as hybridizing mammal species may be. Correspondence to: Ú. Árnason     

基于12S,16S和28S rRNA基因序列的中国小毛瓢虫属系统发育分析

【目的】小毛瓢虫属Scymnus Kugelann昆虫主要捕食蚜虫、蚧虫等害虫,是一类经济上重要的天敌昆虫。目前针对小毛瓢虫属的系统发育研究尚属空白,亚属之间的系统演化关系尚不明确,为了建立合理的分类系统,亟需对小毛瓢虫属的亲缘关系进行研究和探讨。【方法】以华南农业大学馆藏的小毛瓢虫属5亚属共44种为研究对象,采用PCR技术对12S, 16S和28S rRNA基因的部分序列进行扩增;运用MEGA 7.0分析了小毛瓢虫属内12S, 16S和28S rRNA基因的碱基组成,基于K2P模型计算了小毛瓢虫属44种的种间遗传距离;采用最大似然法(maximum-likelihood, ML)和贝叶斯推断法(Bayesian-inference, BI)构建该属的系统发育树。【结果】扩增获得小毛瓢虫属44种的12S rRNA基因序列平均长度为356 bp, 16S rRNA基因序列平均长度为351 bp, 28S rRNA基因序列平均长度为315 bp;序列分析表明,12S rRNA基因的A, T, G和C平均含量分别为38.8%, 43.5%, 11.9%和5.8%, 16S rRNA基因的A, T, G和C平均含量分别为37.6%, 40.3%, 14.4%和7.7%, 28S rRNA基因的A, T, G和C平均含量分别为26.7%, 18.3%, 31.4%和23.5%;基于联合序列分析的种间遗传距离为0.004～0.276,平均遗传距离为0.115。系统发育分析结果表明,小毛瓢虫属为单系起源,而小毛瓢虫亚属Scymnus(Scymnus) Kugelann、毛瓢虫亚属Scymnus(Neopullus) Sasaji、小瓢虫亚属Scymnus(Pullus) Mulsant和拟小瓢虫亚属Scymnus(Parapullus) Yang均为并系起源。【结论】基于12S, 16S和28S rRNA基因序列的小毛瓢虫属系统发育分析显示传统的形态学分类体系与基于分子数据分析的结果部分不一致,这表明应该对该属内各亚属的鉴别特征进行全面检视,筛选并确立各亚属的形态指标,同时也表明该属内的亚属分类单元需重新厘定。     

Nuclease S1 mapping of 16S ribosomal RNA in ribosomes

Escherichia coli 16S rRNA and 16S-like rRNAs from other species have several universally conserved sequences which are believed to be single-stranded in ribosomes. The quantitative disposition of these sequences within ribosomes is not known. Here we describe experiments designed to explore the availability of universal 16S rRNA sequences for hybridization with DNA probes in 30S particles and 70S ribosomes. Unlike previous investigations, quantitative data on the accessibility of DNA probes to the conserved portions of 16S rRNA within ribosomes was acquired. Uniquely, the experimental design also permitted investigation of cooperative interactions involving portions of conserved 16S rRNA. The basic strategy employed ribosomes, 30S subunits, and 16S rRNAs, which were quantitatively analyzed for hybridization efficiency with synthetic DNA in combination with nuclease S1. In deproteinated E. coli 16S rRNA and 30S subunits, the regions 520-530, 1396-1404, 1493-1504, and 1533-1542 are all single-stranded and unrestricted for hybridization to short synthetic DNAs. However, the quantitative disposition of the sequences in 70S ribosomes varies with each position. In 30S subunits there appear to be no cooperative interactions between the 16S rRNA universal sequences investigated.     

Defining reference sequences for Nocardia species by similarity and clustering analyses of 16S rRNA gene sequence data

Background

The intra- and inter-species genetic diversity of bacteria and the absence of ‘reference’, or the most representative, sequences of individual species present a significant challenge for sequence-based identification. The aims of this study were to determine the utility, and compare the performance of several clustering and classification algorithms to identify the species of 364 sequences of 16S rRNA gene with a defined species in GenBank, and 110 sequences of 16S rRNA gene with no defined species, all within the genus Nocardia.

Methods

A total of 364 16S rRNA gene sequences of Nocardia species were studied. In addition, 110 16S rRNA gene sequences assigned only to the Nocardia genus level at the time of submission to GenBank were used for machine learning classification experiments. Different clustering algorithms were compared with a novel algorithm or the linear mapping (LM) of the distance matrix. Principal Components Analysis was used for the dimensionality reduction and visualization.

Results

The LM algorithm achieved the highest performance and classified the set of 364 16S rRNA sequences into 80 clusters, the majority of which (83.52%) corresponded with the original species. The most representative 16S rRNA sequences for individual Nocardia species have been identified as ‘centroids’ in respective clusters from which the distances to all other sequences were minimized; 110 16S rRNA gene sequences with identifications recorded only at the genus level were classified using machine learning methods. Simple kNN machine learning demonstrated the highest performance and classified Nocardia species sequences with an accuracy of 92.7% and a mean frequency of 0.578.

Conclusion

The identification of centroids of 16S rRNA gene sequence clusters using novel distance matrix clustering enables the identification of the most representative sequences for each individual species of Nocardia and allows the quantitation of inter- and intra-species variability.     

Evaluation of different partial 16S rRNA gene sequence regions for phylogenetic analysis of microbiomes

Operational taxonomic units (OTUs) are conventionally defined at a phylogenetic distance (0.03—species, 0.05—genus, 0.10—family) based on full-length 16S rRNA gene sequences. However, partial sequences (700 bp or shorter) have been used in most studies. This discord may affect analysis of diversity and species richness because sequence divergence is not distributed evenly along the 16S rRNA gene. In this study, we compared a set each of bacterial and archaeal 16S rRNA gene sequences of nearly full length with multiple sets of different partial 16S rRNA gene sequences derived therefrom (approximately 440-700 bp), at conventional and alternative distance levels. Our objective was to identify partial sequence region(s) and distance level(s) that allow more accurate phylogenetic analysis of partial 16S rRNA genes. Our results showed that no partial sequence region could estimate OTU richness or define OTUs as reliably as nearly full-length genes. However, the V1-V4 regions can provide more accurate estimates than others. For analysis of archaea, we recommend the V1-V3 and the V4-V7 regions and clustering of species-level OTUs at 0.03 and 0.02 distances, respectively. For analysis of bacteria, the V1-V3 and the V1-V4 regions should be targeted, with species-level OTUs being clustered at 0.04 distance in both cases.     

16S rRNA Gene Sequence Analysis of Drinking Water Using RNA and DNA Extracts as Targets for Clone Library Development

The bacterial composition of chlorinated drinking water was analyzed using 16S rRNA gene clone libraries derived from DNA extracts of 12 samples and compared to clone libraries previously generated using RNA extracts from the same samples. Phylogenetic analysis of 761 DNA-based clone sequences showed that unclassified bacteria were the most abundant group, representing nearly 62% of all DNA sequences analyzed. Other phylogenetic groups identified included Proteobacteria (20%), Actinobacteria (9%), Cyanobacteria (4%), and Bacteroidetes (2%). The composition of RNA-based libraries (1122 sequences) was similar to the DNA-based libraries with a few notable exceptions: Proteobacteria were more dominant in the RNA clone libraries (i.e., 35% RNA; 20% DNA). Differences in the Proteobacteria composition were also observed; alpha-Proteobacteria was 22 times more abundant in the RNA-based clones while beta-Proteobacteria was eight times more abundant in the DNA libraries. Nearly twice as many DNA operational taxonomic units (OTUs) than RNA OTUs were observed at distance 0.03 (101 DNA; 53 RNA). Twenty-four OTUs were shared between all RNA- and DNA-based libraries (OTU_0.03) representing only 18% of the total OTUs, but 81% (1527/1883) of all sequences. Such differences between clone libraries demonstrate the necessity of generating both RNA- and DNA-derived clone libraries to compare these two different molecular approaches for community analyses.     

Discrimination of Bacillus anthracis and closely related microorganisms by analysis of 16S and 23S rRNA with oligonucleotide microarray

Analysis of 16S rRNA sequences is a commonly used method for the identification and discrimination of microorganisms. However, the high similarity of 16S and 23S rRNA sequences of Bacillus cereus group organisms (up to 99-100%) and repeatedly failed attempts to develop molecular typing systems that would use DNA sequences to discriminate between species within this group have resulted in several suggestions to consider B. cereus and B. thuringiensis, or these two species together with B. anthracis, as one species. Recently, we divided the B. cereus group into seven subgroups, Anthracis, Cereus A and B, Thuringiensis A and B, and Mycoides A and B, based on 16S rRNA, 23S rRNA and gyrB gene sequences and identified subgroup-specific makers in each of these three genes. Here we for the first time demonstrated discrimination of these seven subgroups, including subgroup Anthracis, with a 3D gel element microarray of oligonucleotide probes targeting 16S and 23S rRNA markers. This is the first microarray enabled identification of B. anthracis and discrimination of these seven subgroups in pure cell cultures and in environmental samples using rRNA sequences. The microarray bearing perfect match/mismatch (p/mm) probe pairs was specific enough to discriminate single nucleotide polymorphisms (SNPs) and was able to identify targeted organisms in 5min. We also demonstrated the ability of the microarray to determine subgroup affiliations for B. cereus group isolates without rRNA sequencing. Correlation of these seven subgroups with groupings based on multilocus sequence typing (MLST), fluorescent amplified fragment length polymorphism analysis (AFLP) and multilocus enzyme electrophoresis (MME) analysis of a wide spectrum of different genes, and the demonstration of subgroup-specific differences in toxin profiles, psychrotolerance, and the ability to harbor some plasmids, suggest that these seven subgroups are not based solely on neutral genomic polymorphisms, but instead reflect differences in both the genotypes and phenotypes of the B. cereus group organisms.     

First data on the molecular phylogeography of scincid lizards of the genus Mabuya

A 487-bp fragment of the mitochondrial 16S rRNA gene was sequenced in 26 species of the circumtropical lizard genus Mabuya and used to analyze phylogenetic relationships within the genus. The species from Africa and Madagascar formed a monophyletic group relative to the included Asian and South American taxa. The Malagasy species included (M. elegans, M. cf. dumasi, and M. comorensis) did not appear as a monophylum. Combined and separate analysis of the 16S data and additional sequences of the mitochondrial 12S rRNA, ND4, and cytochrome b genes (a total of 2255 bp) in one Asian, two Malagasy, and two African species also did not result consistently in a monophyletic grouping of the Malagasy taxa. However, a monophylum containing African and Malagasy taxa was strongly supported by the combined analysis. These preliminary results indicate that Mabuya probably colonized Madagascar from Africa through the Mozambique Channel.