Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta : thysanoptera): Convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes

To help understand the mechanisms of gene rearrangement in the mitochondrial (mt) genomes of hemipteroid insects, we sequenced the mt genome of the plague thrips, Thrips imaginis (Thysanoptera). This genome is circular, 15,407 by long, and has many unusual features, including (1) rRNA genes inverted and distant from one another, (2) an extra gene for tRNA-Ser, (3) a tRNA-Val lacking a D-arm, (4) two pseudo-tRNA genes, (5) duplicate control regions, and (6) translocations and/or inversions of 24 of the 37 genes. The mechanism of rRNA gene transcription in T. imaginis may be different from that of other arthropods since the two rRNA genes have inverted and are distant from one another. Further, the rRNA genes are not adjacent or even close to either of the two control regions. Tandem duplication and deletion is a plausible model for the evolution of duplicate control regions and for the gene translocations, but intramitochondrial recombination may account for the gene inversions in T. imaginis. All the 18 genes between control regions #1 and #2 have translocated and/or inverted, whereas only six of the 20 genes outside this region have translocated and/or inverted. Moreover, the extra tRNA gene and the two pseudo-tRNA genes are either in this region or immediately adjacent to one of the control regions. These observations suggest that tandem duplication and deletion may be facilitated by the duplicate control regions and may have occurred a number of times in the lineage leading to T. imaginis. T. imaginis shares two novel gene boundaries with a lepidopsocid species from another order of hemipteroid insects, the Psocoptera. The evidence available suggests that these shared gene boundaries evolved by convergence and thus are not informative for the interordinal phylogeny of hemipteroid insects. We discuss the potential of hemipteroid insects as a model system for studies of the evolution of animal rut genomes and outline some fundamental questions that may be addressed with this system.

Novel mitochondrial gene content and gene arrangement indicate illegitimate inter-mtDNA recombination in the chigger mite, Leptotrombidium pallidum

To better understand the evolution of mitochondrial (mt) genomes in the Acari (mites and ticks), we sequenced the mt genome of the chigger mite, Leptotrombidium pallidum (Arthropoda: Acari: Acariformes). This genome is highly rearranged relative to that of the hypothetical ancestor of the arthropods and the other species of Acari studied. The mt genome of L. pallidum has two genes for large subunit rRNA, a pseudogene for small subunit rRNA, and four nearly identical large noncoding regions. Nineteen of the 22 tRNAs encoded by this genome apparently lack either a T-arm or a D-arm. Further, the mt genome of L. pallidum has two distantly separated sections with identical sequences but opposite orientations of transcription. This arrangement cannot be accounted for by homologous recombination or by previously known mechanisms of mt gene rearrangement. The most plausible explanation for the origin of this arrangement is illegitimate inter-mtDNA recombination, which has not been reported previously in animals. In light of the evidence from previous experiments on recombination in nuclear and mt genomes of animals, we propose a model of illegitimate inter-mtDNA recombination to account for the novel gene content and gene arrangement in the mt genome of L. pallidum.

Molecular mechanisms for the variation of mitochondrial gene content and gene arrangement among chigger mites of the genus Leptotrombidium (Acari : Acariformes)

The gene content of a mitochondrial (mt) genome, i.e., 37 genes and a large noncoding region (LNR), is usually conserved in Metazoa. The arrangement of these genes and the LNR is generally conserved at low taxonomic levels but varies substantially at high levels. We report here a variation in mt gene content and gene arrangement among chigger mites of the genus Leptotrombidium. We found previously that the mt genome of Leptotrombidium pallidum has an extra gene for large-subunit rRNA (rrnL), a pseudo-gene for small-subunit rRNA (PrrnS), and three extra LNRs, additional to the 37 genes and an LNR typical of Metazoa. Further, the arrangement of mt genes of L. pallidum differs drastically from that of the hypothetical ancestor of the arthropods. To find to what extent the novel gene content and gene arrangement occurred in Leptotrombidium, we sequenced the entire or partial mt genomes of three other species, L. akamushi, L. deliense, and L. fletcheri. These three species share the arrangement of all genes with L. pallidum, except trnQ (for tRNA-glutamine). Unlike L. pallidum, however, these three species do not have extra rrnL or PrrnS and have only one extra LNR. By comparison between Leptotrombidium species and the ancestor of the arthropods, we propose that (1) the type of mt genome present in L. pallidum evolved from the type present in the other three Leptotrombidium species, and (2) three molecular mechanisms were involved in the evolution of mt gene content and gene arrangement in Leptotrombidium species.

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: A case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another I I species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.

Muscle-specific regulation of tropomyosin gene expression and myofibrillogenesis differs among muscle systems examined at metamorphosis of the gastropod Haliotis rufescens

The spatial and temporal association of muscle-specific tropomyosin gene expression, and myofibril assembly and degradation during metamorphosis is analyzed in the gastropod mollusc. Haliotis rufescens. Metamorphosis of tile planktonic larva to the benthic juvenile includes rearrangement and atrophy of specific larval muscles, and biogenesis of the new juvenile muscle system. The major muscle of the larva - the larval retractor muscle - reorganizes at metamorphosis, with two suites of cells having different fates. The ventral cells degenerate, while the dorsal cells become part of the developing juvenile mantle musculature. Prior to these changes in myofibrillar structure, tropomyosin mRNA prevalence declines until undetectable in the ventral cells, while increasing markedly in the dorsal cells. In the foot muscle and right shell muscle, tropomyosin mRNA levels remain relatively stable, even trough myofibril content increases. In a population of median mesoderm cells destined to form de novo the major muscle of the juvenile and adult (the columellar muscle), tropomyosin expression is initiated at 45 h after induction of metamorphosis. Myofibrillar filamentous actin is not detected in these cells until about 7 days later. Given that patterns of tropomyosin mRNA accumulation in relation to myofibril assembly and disassembly differ significantly among the four major muscle systems examined, we suggest that different regulatory mechanisms, probably operating at both transcriptional and post-transcriptional levels, control the biogenesis and atrophy of different larval and postlarval muscles at metamorphosis.

Molecular phylogeny of Wolbachia endosymbionts in southeast asian mosquitoes (Diptera: Culicidae) based on wsp gene sequences

Wolbachia are maternally inherited intracellular bacteria that infect a wide range of arthropods and nematodes and are associated with various reproductive abnormalities in their hosts. Insect-associated Wolbachia form a monophyletic clade in the α-Proteobacteria and recently have been separated into two supergroups (A and B) and 19 groups. Our recent polymerase chain reaction (PCR) survey using wsp specific primers indicated that various strains of Wolbachia were present in mosquitoes collected from Southeast Asia. Here, we report the phylogenetic relationship of the Wolbachia strains found in these mosquitoes using wsp gene sequences. Our phylogenetic analysis revealed eight new Wolbachia strains, five in the A supergroup and three in the B supergroup. Most of the Wolbachia strains present in Southeast Asian mosquitoes belong to the established Mors, Con, and Pip groups.

wsp Gene sequences from the Wolbachia of filarial nematodes

Wolbachia endosymbiotic bacteria are widespread in arthropods and are also present in filarial nematodes. Almost all filarial species so far examined have been found to harbor these endosymbionts. The sequences of only three genes have been published for nematode Wolbachia (i.e., the genes coding for the proteins FtsZ and catalase and for 16S rRNA). Here we present the sequences of the genes coding for the Wolbachia surface protein (WSP) from the endosymbionts of eight species of filaria. Complete gene sequences were obtained from the endosymbionts of two different species, Dirofilaria immitis and Brugia malayi. These sequences allowed us to design general primers for amplification of the wsp gene from the Wolbachia of all filarial species examined. For these species, partial WSP sequences (about 600 base pairs) were obtained with these primers. Phylogenetic analysis groups these nematode wsp sequences into a coherent cluster. Within the nematode cluster, wsp-based Wolbachia phylogeny matches a previous phylogeny obtained with ftsZ gene sequences, with a good consistency of the phylogeny of hosts (nematodes) and symbionts (Wolbachia). In addition, different individuals of the same host species (Dirofilaria immitis and Wuchereria bancrofti) show identical wsp gene sequences.

Gene organization of the dnaA region of Wolbachia

The dnaA region of Wolbachia, an intracellular bacterial parasite of insects, is unique. A glnA cognate was found upstream of the dnaA gene, while neither of the two open reading frames detected downstream of dnaA has any homologue in the database. This unusual gene arrangement may reflect requirements associated with the unique ecological niche this agent occupies.

Cloning and characterization of a gene encoding the major surface protein of the bacterial endosymbiont Wolbachia pipientis

The maternally inherited intracellular symbiont Wolbachia pipientis is well known for inducing a variety of reproductive abnormalities in the diverse arthropod hosts it infects. It has been implicated in causing cytoplasmic incompatibility, parthenogenesis, and the feminization of genetic males in different hosts. The molecular mechanisms by which this fastidious intracellular bacterium causes these reproductive and developmental abnormalities have not yet been determined. In this paper, we report on (i) the purification of one of the most abundantly expressed Wolbachia proteins from infected Drosophila eggs and (ii) the subsequent cloning and characterization of the gene (wsp) that encodes it. The functionality of the wsp promoter region was also successfully tested in Escherichia coli. Comparison of sequences of this gene from different strains of Wolbachia revealed a high level of variability. This sequence variation correlated with the ability of certain Wolbachia strains to induce or rescue the cytoplasmic incompatibility phenotype in infected insects. As such, this gene will be a very useful tool for Wolbachia strain typing and phylogenetic analysis, as well as understanding the molecular basis of the interaction of Wolbachia with its host.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Evidence for persistent dysfunction of wild-type aldosterone synthase gene in glucocorticoid-treated familial hyperaldosteronism type I

Background In familial hyperaldosteronism type I (FH-I), glucocorticoid treatment suppresses adrenocorticotrophic hormone-regulated hybrid gene expression and corrects hyperaldosteronism. Objective To determine whether the wild-type aldosterone synthase genes, thereby released from chronic suppression, are capable of functioning normally. Methods We compared mid-morning levels of plasma potassium, plasma aldosterone, plasma renin activity (PRA) and aldosterone : PRA ratios, measured with patients in an upright position, and responsiveness of aldosterone levels to infusion of angiotensin II (AII), for 11 patients with FH-I before and during long-term (0.8-14.3 years) treatment with 0.25-0.75 mg/day dexamethasone or 2.5-10 mg/day prednisolone. Results During glucocorticoid treatment, hypertension was corrected in all. Potassium levels, which had been low (< 3.5 mmol/l) in two patients before treatment, were normal in all during treatment (mean 4.0 +/- 0.1 mmol/l, range 3.5-4.6). Aldosterone levels during treatment [13.2 +/- 2.1 ng/100 ml (mean +/- SEM)] were lower than those before treatment (20.1 +/- 2.5 ng/100 ml, P < 0.05). PRA levels, which had been suppressed before treatment (0.5 +/- 0.2 ng/ml per h), were unsuppressed during treatment (5.1 +/- 1.5 ng/ml per h, P < 0.01) and elevated (> 4 ng/ml per h) in six patients. Aldosterone : PRA ratios, which had been elevated (> 30) before treatment (101.1 +/- 25.9), were much lower during treatment (4.1 +/- 1.0, P < 0.005) and below normal (< 5) in eight patients. Surprisingly, aldosterone level, which had not been responsive (< 50% rise) to infusion of AII for all 11 patients before treatment, remained unresponsive for 10 during treatment. Conclusions Apparently regardless of duration of glucocorticoid treatment in FH-I, aldosterone level remains poorly responsive to AII, with a higher than normal PRA and a low aldosterone : PRA ratio. This is consistent with there being a persistent defect in functioning of wild-type aldosterone synthase gene. (C) Rapid Science Publishers ISSN 0263-6352.