Mitochondrial DNA mutations in patients with postlingual, nonsyndromic hearing impairment

Mitochondrial mutations have previously been reported anecdotally in families with maternally inherited, nonsyndromic hearing impairment. To ascertain the contribution of mitochondrial mutations to postlingual but early-onset, nonsyndromic hearing impairment, we screened patients collected from within two different populations (southern Italy and UK) for previously reported mtDNA mutations associated with hearing disorders. Primer extension (SNP analysis) was used to screen for specific mutations, revealing cases of heteroplasmy and its extent. The most frequently implicated tRNA genes, Leu(UUR) and Ser(UCN), were also sequenced in all Italian patients. All tRNA genes were sequenced in those UK patients showing the clearest likelihood of maternal inheritance. Causative mtDNA mutations were found in approximately 5% of patients in both populations, representing almost 10% of cases that were clearly familial. Age of onset, where known, was generally before adulthood, and hearing loss was typically progressive. Haplogroup analysis revealed a possible excess of haplogroup cluster HV in the patients, compared with population controls, but of borderline statistical significance. In contrast, we did not find any of the previously reported mtDNA mutations, nor a significant deviation from haplogroup cluster frequencies typical of the control population, in patients with late adult-onset hearing loss (age-related hearing impairment) from the UK or Finland.     

Clustering of tRNA genes in Paracentrotus lividus mitochondrial DNA

Summary We have determined the base sequence of the restriction fragment Bam1-2 (3,593) of Paracentrotus lividus (sea urchin) mtDNA. This fragment contains, in addition to genes previously identified (part of the 12S rRNA, ND1 and part of the ND2 mRNA), a cluster of 15 tRNA genes located between the 12S and ND1 genes. Also to be found in the tRNA gene cluster, between the tRNA^Thr and tRNA^Pro genes, is a sequence of 134 bp which constitutes the only non-coding region of this DNA so far identified. The distinctive organization of the tRNA genes and the extreme size reduction of the non-coding region suggest the existence of unique mechanisms for the regulation of gene expression in this organism.     

Phenylalanine and tyrosine transfer RNAs encoded by Tetrahymena pyriformis mitochondrial DNA: primary sequence,post-transcriptional modifications,and gene localization

We have isolated Phe and Tyr tRNAs from Tetrahymena pyriformis mitochondria and have determined that these are "native" species, encoded by the mtDNA. A single gene for the tRNA(Phe) has been positioned 12-14 kbp from the left end of the linear Tetrahymena mtDNA, while duplicate tRNA(Tyr) genes have been localized within the inverted terminal repeats of this genome. Primary sequence analysis demonstrates that the tRNA(Tyr) has all of the characteristic primary and secondary structural features of a normal tRNA; however, the tRNA(Phe) displays several atypical features, including (i) replacement of the usual T psi sequence by UC, (ii) a U.U pair in the T psi C stem, and (iii) an extra 5'-nucleotide (U).     

Enhanced ROS production and antioxidant defenses in cybrids harbouring mutations in mtDNA

It has been suggested that mutations in mitochondrial DNA (mtDNA) can produce an increase in reactive oxygen species (ROS) and that this can play a major role in the pathogenic mechanisms of mitochondrial encephalomyopathies. Many studies exist using electron transport chain (ETC) inhibitors, however there are only a few studies that examine ROS production associated with mutations in the mtDNA. To investigate this issue, we have studied ROS production, antioxidant defences and oxidative damage to lipids and proteins in transmitochondrial cybrids carrying different mtDNA mutations. Here, we report that two different mutant cell lines carrying mutations in their mitochondrial tRNA genes (A3243G in tRNA LeuUUR and A8344G in tRNA Lys) showed an increased ROS production with a parallel increase in the antioxidant enzyme activities, which may protect cells from oxidative damage in our experimental conditions (no overt oxidative damage to lipids and proteins has been observed). In contrast, cytochrome c oxidase (COX) mutant cybrids (carrying the stop-codon mutation G6930A in the COXI gene) showed neither an increase in ROS production nor elevation of antioxidant enzyme activities or oxidative damage. These results suggest that the specific location of mutations in mtDNA has a strong influence on the phenotype of the antioxidant response. Therefore, this issue should be carefully considered when antioxidant therapies are investigated in patients with mitochondrial disorders.     

Somatic mutations of the mitochondrial genome in human breast cancers

Somatic mutations in mitochondrial DNA (mtDNA) have been identified in various tumors, including breast cancer. However, their clinicopathological impact on breast cancer still remains unclear. In this study, we re-sequenced the entire mtDNA from breast cancer samples together with paired non-tumorous breast tissues from 58 Taiwanese patients. We identified 19 somatic mutations in the mtDNA coding region of 16 breast cancers. Out of these mutations, 12 of the 19 mutations (63%) are missense or frame-shift mutations that have the potential to cause mitochondrial dysfunction. In combination with our previously study on the D-loop region of mtDNA, we found that 47% (27/58) of the breast cancers harbored somatic mtDNA mutations. Among a total of 40 somatic mutations, 53% (21/40) were located in the D-loop region of the mtDNA, 5% (2/40) were in the ribosomal RNA genes, 5% (2/40) were in the tRNA genes, and 38% (15/40) occurred in mRNA genes. The occurrence of these somatic mtDNA mutations is associated with an older onset age (≥ 50-year old, P = 0.039), a higher TNM stage (P = 0.027), and a higher histological grade (P = 0.012). Multiple logistic regression analysis revealed that an older onset age (P = 0.029) and a higher histological grade (P = 0.006) are significantly correlated with patients having somatic mutations in the mtDNA in their breast cancer sample. In conclusion, our results suggest that somatic mtDNA mutations may play a critical role in the progression of breast cancer.     

A functionally dominant mitochondrial DNA mutation

Mutations in mitochondrial DNA (mtDNA) tRNA genes can be considered functionally recessive because they result in a clinical or biochemical phenotype only when the percentage of mutant molecules exceeds a critical threshold value, in the range of 70-90%. We report a novel mtDNA mutation that contradicts this rule, since it caused a severe multisystem disorder and respiratory chain (RC) deficiency even at low levels of heteroplasmy. We studied a 13-year-old boy with clinical, radiological and biochemical evidence of a mitochondrial disorder. We detected a novel heteroplasmic C>T mutation at nucleotide 5545 of mtDNA, which was present at unusually low levels (<25%) in affected tissues. The pathogenic threshold for the mutation in cybrids was between 4 and 8%, implying a dominant mechanism of action. The mutation affects the central base of the anticodon triplet of tRNA(Trp) and it may alter the codon specificity of the affected tRNA. These findings introduce the concept of dominance in mitochondrial genetics and pose new diagnostic challenges, because such mutations may easily escape detection. Moreover, similar mutations arising stochastically and accumulating in a minority of mtDNA molecules during the aging process may severely impair RC function in cells.     

Assembly of the mitochondrial membrane system. Organization of yeast mitochondrial DNA in the Oli1 region

Summary The mitochondrial DNA (mtDNA) of a cytoplasmic petite mutant (DS401) of Saccharomyces cerevisiae genetically marked for the ATPase proteolipid, serine tRNA and varl genes has been characterized by restriction endonuclease analysis and DNA sequencing. The DS401 mtDNA segment is 5.3 kb long spanning the region between 79.1 and 86.8 units of the wild type genome. Most of the DS401 mtDNA consists of A+T rich sequences. In addition, however, there are ten short sequences with a high content of G+C and two sequences that have been identified as the ATPase proteolipid and the serine tRNA genes. The two genes map at 81 and 83 units and are transcribed from the same DNA strand. Even though there are other possible coding sequences in the DNA segment, none are sufficiently long to code for a gene product of the size of the varl protein. Based on the relative organization of the G+C rich clusters and genes, a model has been proposed for the processing of mitochondria) RNA. This model postulates the existence of mitochondrial double strand specific RNases that cleave the RNA at the G+C clusters.     

Respiratory chain complex I deficiency caused by mitochondrial DNA mutations

Defects of the mitochondrial respiratory chain are associated with a diverse spectrum of clinical phenotypes, and may be caused by mutations in either the nuclear or the mitochondrial genome (mitochondrial DNA (mtDNA)). Isolated complex I deficiency is the most common enzyme defect in mitochondrial disorders, particularly in children in whom family history is often consistent with sporadic or autosomal recessive inheritance, implicating a nuclear genetic cause. In contrast, although a number of recurrent, pathogenic mtDNA mutations have been described, historically, these have been perceived as rare causes of paediatric complex I deficiency. We reviewed the clinical and genetic findings in a large cohort of 109 paediatric patients with isolated complex I deficiency from 101 families. Pathogenic mtDNA mutations were found in 29 of 101 probands (29%), 21 in MTND subunit genes and 8 in mtDNA tRNA genes. Nuclear gene defects were inferred in 38 of 101 (38%) probands based on cell hybrid studies, mtDNA sequencing or mutation analysis (nuclear gene mutations were identified in 22 probands). Leigh or Leigh-like disease was the most common clinical presentation in both mtDNA and nuclear genetic defects. The median age at onset was higher in mtDNA patients (12 months) than in patients with a nuclear gene defect (3 months). However, considerable overlap existed, with onset varying from 0 to >60 months in both groups. Our findings confirm that pathogenic mtDNA mutations are a significant cause of complex I deficiency in children. In the absence of parental consanguinity, we recommend whole mitochondrial genome sequencing as a key approach to elucidate the underlying molecular genetic abnormality.     

Mitochondrial DNA variations in patients with Type 2 (non‐insulin dependent) diabetes mellitus and a Welsh control population

Type 2 (non‐insulin dependent) diabetes mellitus may be inherited along the maternal line and a variety of mitochondrial DNA (mtDNA) variants have been implicated in the pathogenesis. We have previously reported mutations in five regions of the mitochondrial genome which encompass 11 of the 22 tRNA genes. Now we employ the technique of single stranded conformational polymorphism (SSCP) analysis to investigate a further 6 regions of the mitochondrial genome, covering the remaining 11 tRNA genes in 40 patients with Type 2 diabetes and 30 racially‐matched normal controls. A variety of homoplasmic mutations were detected in patients with diabetes and these will be of value in further population association studies. Hum Mutat 13:412–413, 1999. © 1999 Wiley‐Liss, Inc.     

Mitochondrial DNA polymorphisms in bipolar disorder

BACKGROUND: Previous studies suggested mitochondrial abnormality in bipolar disorder: (1) possible contribution of parent-of-origin effect in transmission of bipolar disorder; (2) abnormal brain phosphorus metabolism detected by phosphorus-31 magnetic resonance spectroscopy; (3) comorbidity of affective disorders in patients with mitochondrial encephalopathy; (4) increased levels of the 4977bp deletion of mitochondrial DNA (mtDNA) in the postmortem brains. We investigated mtDNA polymorphisms in association with bipolar disorder. METHODS: Twelve PCR fragments including all tRNA genes were examined by the single-strand conformation polymorphism method in 43 bipolar patients. All observed polymorphisms were sequenced. Association of these polymorphisms with bipolar disorder was examined by restriction fragment length polymorphism method in 135 bipolar patients and 187 controls. RESULTS: In total, we found 28 polymorphisms including 14 polymorphisms that have not been reported previously. The A10398G polymorphism was significantly associated with bipolar disorder (10398A genotype: 33.1% in bipolar, 22.2% in the control, P<0.05). Although this difference was not significant after Bonferroni correction, the CA haplotype of the 5178 and 10398 polymorphisms was still significantly associated with bipolar disorder (CA haplotype: 33.6% in bipolar, 16.8% in control, P<0.001). Three rare mutations substituting evolutionary conserved bases; A5539G in tRNA(Trp) gene, A5747G in the origin of L-strand replication, and A8537G in ATPase subunit-6 and -8 genes, were found in patients with family history in which maternal transmission was suspected. DISCUSSION: The 5178C/10398A haplotype in mtDNA may be a risk factor of bipolar disorder (odds ratio, 2.4). Pathophysiological significance of rare mtDNA mutations needs to be verified in the future. This finding may imply the pathophysiological significance of mtDNA in bipolar disorder.     

Variants in mitochondrial tRNAGlu, tRNAArg, and tRNAThr may influence the phenotypic manifestation of deafness-associated 12S rRNA A1555G mutation in three Han Chinese families with hearing loss

We report here on the clinical, genetic, and molecular characterization of three Han Chinese pedigrees with aminoglycoside-induced and nonsyndromic hearing loss. Clinical evaluation revealed the variable phenotype of hearing loss including severity, age-at-onset, audiometric configuration in these subjects. Penetrances of hearing loss in BJ107, BJ108, and BJ109 pedigrees are 35%, 63%, and 67%, respectively. Mutational analysis of the complete mitochondrial genomes in these pedigrees showed the identical homoplasmic A1555G mutation and distinct sets of mitochondrial DNA (mtDNA) variants belonging to haplogroups N, F, and M, respectively. Of these variants, the A14693G mutation in the tRNA(Glu), the T15908C mutation in the tRNA(Thr), and the T10454C mutation in the tRNA(Arg) are of special interest as these mutations occur at positions which are highly evolutionarily conserved nucleotides of corresponding tRNAs. These homoplasmic mtDNA mutations were absent among 156 unrelated Chinese controls. The A14693G and T10454C mutations occur at the highly conserved bases of the TpsiC-loop of tRNA(Glu) and tRNA(Arg), respectively. Furthermore, the T15908C mutation in the tRNA(Thr) disrupts a highly conserved A-U base-pairing at the D-stem of this tRNA. The alteration of structure of these tRNAs by these mtDNA mutations may lead to a failure in tRNA metabolism, thereby causing impairment of mitochondrial translation. Thus, mitochondrial dysfunctions, caused by the A1555G mutation, would be worsened by these mtDNA mutations. Therefore, these mtDNA mutations may have a potential modifier role in increasing the penetrance and expressivity of the deafness-associated 12S rRNA A1555G mutation in those Chinese pedigrees.     

Gene organization of the mitochondrial DNA of yeasts: Kluyveromyces lactis and Saccharomycopsis lipolytica

Summary Restriction fragment maps have been constructed for the mitochondrial DNA from two petitenegative yeasts, Kluyveromyces lactis and Saccharomycopsis lipolytica (Candida lipolytica). On these circular genomes, we localized the sequences homologous to the S. cerevisiae mtDNA fragments carrying known genes. The arrangement of genes for ATPase subunit proteins, ribosomal RNA and 4S RNA shows a common feature with respect to S. cerevisiae mitochondrial genome.Abbreviations bp base pairs - mtDNA mitochondrial DNA - tRNA transfer RNA - rRNA ribosomal RNA     

Novel mitochondrial DNA mutations responsible for maternally inherited nonsyndromic hearing loss

Some cases of maternally inherited isolated deafness are caused by mtDNA mutations, frequently following an exposure to aminoglycosides. Two mitochondrial genes have been clearly described as being affected by mutations responsible for this pathology: the ribosomal RNA 12S gene and the transfer RNA serine (UCN) gene. A previous study identified several candidate novel mtDNA mutations, localized in a variety of mitochondrial genes, found in patients with no previous treatment with aminoglycosides. Five of these candidate mutations are characterized in the present study. These mutations are localized in subunit ND1 of complex I of the respiratory chain (m.3388C>A [p.MT-ND1:Leu28Met]), the tRNA for Isoleucine (m.4295A>G), subunit COII of complex IV (m.8078G>A [p.MT-CO2:Val165Ile]), the tRNA of Serine 2 (AGU/C) (m.12236G>A), and Cytochrome B, subunit of complex III (m.15077G>A [p.MT-CYB:Glu111Lys]). Cybrid cell lines have been constructed for each of the studied mtDNA mutations and functional studies have been performed to assess the possible consequences of these mutations on mitochondrial bioenergetics. This study shows that a variety of mitochondrial genes, including protein-coding genes, can be responsible for nonsyndromic deafness, and that exposure to aminoglycosides is not required to develop the disease, giving new insights on the molecular bases of this pathology.     

Contrasting phenotypes in three patients with novel mutations in mitochondrial tRNA genes

We studied three patients, each harboring a novel mutation at a highly conserved position in a different mitochondrial tRNA gene. The mutation in patient 1 (T5543C) was associated with isolated mitochondrial myopathy, and occurred in the anticodon loop of tRNA(Trp). In patient 2, with mitochondrial myopathy and marked retinopathy, the mutation (G14710A) resulted in an anticodon swap (Glu to Lys) in tRNA(Glu). Patient 3, who manifested mitochondrial encephalomyopathy and moderate retinal dysfunction, harbored a mutation (C3287A) in the TpsiC loop of tRNA(Leu(UUR)). The mutations were heteroplasmic in muscle in all cases, and sporadic in two cases. PCR-RFLP analysis in all patients showed much higher amounts of mutated mtDNA in affected tissue (muscle) than unaffected tissue (blood), and significantly higher levels of mutated mtDNA in cytochrome c oxidase (COX)-negative muscle fibers than in COX-positive fibers, confirming the pathogenicity of these mutations. The mutation was also detected in single hair roots from all three patients, indicating that each mutation must have arisen early in embryonic development or in maternal germ cells. This suggests that individual hair root analyses may reflect a wider tissue distribution of mutated mtDNA than is clinically apparent, and might be useful in predicting prognosis and, perhaps, the risk of transmitting the mutation to offspring. Our data suggest a correlation between clinical phenotype and distribution of mutated mtDNA in muscle versus hair roots. Furthermore, the high threshold for phenotypic expression in single muscle fibers (92-96%) suggests that therapies may only need to increase the percentage of wild-type mtDNA by a small amount to be beneficial.     

Phylogenetic analysis of the mitochondrial genome indicates significant differences between patients with Alzheimer disease and controls in a French-Canadian founder population.

The activity of cytochrome oxidase (CO), the terminal enzyme of the mitochondrial electron transport chain, has been reported to be lower in the brains of Alzheimer disease (AD) patients. This suggests that a modification of mitochondrial DNA (mtDNA) may be responsible for this decrease of CO activity. Many mtDNA variants were found by different studies at a higher frequency in AD patients, suggesting that mtDNA variants could confer a genetic susceptibility to AD. In this study, we sequenced the entire mitochondrial genome region that encompasses the three CO genes and the 22 mitochondrial tRNA in 69 AD patients and 83 age-matched controls. We detected a total of 95 mtDNA variants. The allele frequencies of the majority of these variants were similar in patients and controls. However, a haplotype composed of three different modifications (positions: 5633, 7476, and 15812) was present in three of the 69 late-onset AD patients (4.3%) and also in 1 of 16 early-onset AD patients (6.2%) but not in control individuals. Given that one of these variants (15812) has already been shown to be associated with another neurodegenerative disease and that all three modifications are relatively conserved and their frequencies in the general population is only 0.1%, our data suggest that the presence of this haplotype may represent a risk factor for AD. We also found a significant association (P < 0.05) of two other variants at positions 709 (rRNA 12S) and 15928 (tRNA(Thr)). These two mtDNA variants are three times more frequent in control individuals compared with AD patients, suggesting that they may be protective against AD.     

Anatomy of amplified mitochondrial DNA in “Ragged” mutants of Aspergillus amstelodami: Excision points within protein genes and a common 215 bp segment containing a possible origin of replication

Summary The extranuclearly-inherited ragged growth phenotype (Rgd) of Aspergillus amstelodami is always accompanied by excision and head-to-tail amplification of mtDNA sequences. In one mutant strain (Rgd1) the amplified mtDNA segment (rgd1 DNA, monomeric length 0.9 kb) maps downstream of the large subunit ribosomal RNA gene (Region 1), whereas in all other strains analyzed the amplified sequences (rdg3-7DNA) are located in Region 2 between genes coding for cytochrome b and ATPase subunit 6. The various region 2 sequences differ in lengths (1.5 to 2.7 kb) but have in common a 215 bp sequence mapping between an. unidentified protein gene (corresponding to URF4 of human mtDNA) and an arginine tRNA gene. This common sequence may contain an origin of replication, because a looped-out hairpin structure similar to that of yeast and human mitochondrial origin sequences can be formed. Furthermore, Region 2 DNA suppresses replication of Region 1 DNA, indicating that the former group of molecules contains the more efficient origin. The nucleotide sequence of the rgd6 repeat unit starts and ends within protein genes of mtDNA, and no homologies were found between heads and tails or their flanking sequences.Abbreviations mtDNA DNA isolated from DNase — treated mitochondria - Rgd ragged mutant strain - rgdDNA highly-reiterated DNA sequences isolated simultaneously with the wild-type genome from mitochondria of ragged mutants - bp base pairs - kb kilobase pairs - URF unassigned reading frame     

Molecular analysis of mitochondrial gene mutations in Korean patients with nonsyndromic hearing loss

Mutations in mitochondrial DNA (mtDNA) are a major cause of hearing loss. In this study, we performed a systematic mutational screening of the 12S rRNA, tRNA Ser(UCN), tRNA Lys and tRNA Leu(UUR) genes in 227 unrelated patients with nonsyndromic hearing impairment for the first time in a Korean population. We found two individuals with an A1555G mutation, which is a frequency (0.9%) lower than that of other East Asians. Furthermore, two novel variants (C895T and 961-CC insertion) in the 12S rRNA gene were identified in the affected individuals, but were absent in 217 controls, indicating that they may play a role in the pathogenesis of hearing loss. Notably, 961delT and T1005C mutations were identified at similar frequencies in both patients and control subjects. Our data suggest that these variants seem to be polymorphisms rather than causes of disease. On the other hand, we did not find any of the known deafness-associated mutations in these tRNA genes. These data suggest that the 12S rRNA gene may be a hot spot for mitochondrial mutations causing hearing loss in the Korean population.     

tRNA genes in rat liver mitochondrial DNA

Summary Continuing the analysis of rat liver mitochondria) DNA, we have determined the sequences of two segments containing several tRNA genes. One of these segments also comprises the L-strand origin of replication. The characteristics of the 15 rat liver mitochondrial tRNA structures known at present are compared with those from other mammalian mitochondria.Abbreviations mtDNA mitochondrial DNA - U.R.F. unidentified reading frame (Anderson et al. 1981) - bp and Kbp base pairs and Kilobase pairs - Pu and Py purine and pyrimidine residues     

Somatic mitochondrial DNA mutations in Chinese patients with osteosarcoma

Somatic mutations in mitochondrial DNA (mtDNA) have been long proposed to drive the pathogenesis and progression of human malignancies. Previous investigations have revealed a high frequency of somatic mutations in the D‐loop control region of mtDNA in osteosarcoma. However, little is known with regard to whether or not somatic mutations also occur in the coding regions of mtDNA in osteosarcoma. To test this possibility, in the present study we screened somatic mutations over the full‐length mitochondrial genome of 31 osteosarcoma tumour tissue samples, and corresponding peripheral blood samples from the same cohort of patients. We detected a sum of 11 somatic mutations in the mtDNA coding regions in our series. Nine of them were missense or frameshift mutations that have the potential to hamper mitochondrial respiratory function. In combination with our earlier observations on the D‐loop fragment, 71.0% (22/31) of patients with osteosarcoma carried at least one somatic mtDNA mutation, and a total of 40 somatic mutations were identified. Amongst them, 29 (72.5%) were located in the D‐loop region, two (5%) were in the sequences of the tRNA genes, two (5%) were in the mitochondrial ATP synthase subunit 6 gene and seven (17.5%) occurred in genes encoding components of the mitochondrial respiratory complexes. In addition, somatic mtDNA mutation was not closely associated with the clinicopathological characteristics of osteosarcoma. Together, these findings suggest that somatic mutations are highly prevalent events in both coding and non‐coding regions of mtDNA in osteosarcoma. Some missense and frameshift mutations are putatively harmful to proper mitochondrial activity and might play vital roles in osteosarcoma carcinogenesis.     

Mitochondrial transfer RNA genes in a black fly, Simulium vittatum (Diptera: Simuliidae), indicate long divergence from mosquito (Diptera: Culicidae) and fruit fly (Diptera: Drosophilidae).

Sequences are given for nine complete genes and one partial mitochondrial tRNA gene of the black fly, Simulium vittatum (Zetterstedt). Sequenced tRNA genes were for alanine, arginine, asparagine, aspartic acid, glutamic acid, glycine, leucine(CUN), lysine, serine(AGN), and valine. Nucleotides were aligned with the same previously sequenced genes in Aedes albopictus Skuse and Drosophila yakuba Burla. A cluster of six tRNA genes, which differ in arrangement in Ae. albopictus and D. yakuba, was amplified by PCR and found to have the same position and orientation in S. vittatum as in D. yakuba. Overall, similarity with either D. yakuba or Ae. albopictus was 86%. Sequences that were common to the three insects suggest that black flies and mosquitoes are as divergent from each other as either is from Drosophila. Sequences for nine species of black flies were obtained for tRNA leucine(CUN) from DNA amplified with another primer set. Little variation occurred within the tRNA gene but, by including the flanking regions to provide 175 base pairs, a phylogeny of the nine species was obtained that was largely consistent with current classification.