Analysis of directional mutation pressure and nucleotide content in mitochondrial cytochrome b genes

The Mitochondrial Genome of Melipona fasciculata (Apidae, Meliponini): Genome Organization and Comparative Analyses, Phylogenetic Implications and Divergence Time Estimations

The native stingless bee Melipona fasciculata is economically and ecologically important to the Brazilian Northeast, providing a sustainable source of income to family farmers and being considered an effective pollinator in most ecosystems and crops. This study describes, for the first time, the mitogenome of the species and its phylogenetic position. The mitochondrial genome was sequenced using a MiSeq Sequencer (Illumina Inc.) and compared with other GenBank bee mitogenomes. The length of the mitochondrial DNA, excluding most of the control region, is 14,753 bp, and contains 13 protein-coding genes (PCGs), 21 transfer RNAs, 2 ribosomal RNAs (16S and 12S), and 1 AT-rich region. The GC-content of the M. fasciculata mitogenome was 13.4%. Of the 36 coding regions, 12 tRNAs and 9 PCGs were encoded in the heavy strand, and 9 tRNAs, 4 PCGs and 2 rRNAs were encoded in the light strand. The relative orientation and gene order was the same as other stingless bee mitogenomes. Phylogenetic inference produced well-resolved relationships with high statistical support for concordant branch topologies, under different optimization schemes and model parameters, within and among Melipona, Bombus, Apis, and related clades of Hymenoptera. In general, our divergence time estimates, which were based on the concatenated gene sequences (PCGs + rRNAs) from various groups, overlapped estimations captured by Bayesian analysis from different studies. The divergence time among Melipona species was estimated to occur during the Oligocene, approximately 24 Mya (95% HPD 14-36 Mya). Our results represent a valuable addition to help understanding not only the taxonomy and evolution of Brazilian stingless bee species, but also to uncover historical dispersal and isolation patterns in Meliponinae.

Distinguishing Genetic Drift from Selection in Papillomavirus Evolution

Pervasive purifying selection on non-synonymous substitutions is a hallmark of papillomavirus genome history, but the role of selection on and the drift of non-coding DNA motifs on HPV diversification is poorly understood. In this study, more than a thousand complete genomes representing Alphapapillomavirus types, lineages, and SNP variants were examined phylogenetically and interrogated for the number and position of non-coding DNA sequence motifs using Principal Components Analyses, Ancestral State Reconstructions, and Phylogenetic Independent Contrasts. For anciently diverged Alphapapillomavirus types, composition of the four nucleotides (A, C, G, T), codon usage, trimer usage, and 13 established non-coding DNA sequence motifs revealed phylogenetic clusters consistent with genetic drift. Ancestral state reconstruction and Phylogenetic Independent Contrasts revealed ancient genome alterations, particularly for the CpG and APOBEC3 motifs. Each evolutionary analytical method we performed supports the unanticipated conclusion that genetic drift and different evolutionary drivers have structured Alphapapillomavirus genomes in distinct ways during successive epochs, even extending to differences in more recently formed variant lineages.

Comparative analysis of the Liriomyza chinensis mitochondrial genome with other Agromyzids reveals conserved genome features

Liriomyza chinensis is a serious pest of onions in many countries, especially in East Asia. We sequenced the complete mitochondrial genome of this species and compared it with five other Agromyzidae species. The L. chinensis mitogenome is a double-stranded 16,175 bp circular molecule with an A + T content of 78.3%. It contains 37 genes and a control region as do the sequenced Liriomyza species. The mitogenomes of L. chinensis and other Agromyzidae species showed a clear bias in nucleotide composition with a positive AT-skew. Most PCGs used standard ATN as start codons, and TAN as termination codons. The tRNAs exhibited the typical clover-leaf structure, except for tRNASer^(AGN) and the two rRNA genes are conserved with those of other Agromyzids. The L. chinensis mitogenome control region included several conserved regions, including a poly-T, two (TA)n and one poly-A stretch, which are considered important replication and transcription. The 13 PCGs were used to study the phylogeny of L. chinensis and five related Agromyzids. Analysis by maximum likelihood, Bayesian inference and genetic distance suggest congruent phylogenetic relationships in Liriomyza spp. and provide a useful supplement to taxonomic classification by morphology.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Population Genetic Structure and Potential Incursion Pathways of the Bluetongue Virus Vector Culicoides brevitarsis (Diptera: Ceratopogonidae) in Australia

Culicoides brevitarsis is a vector of the bluetongue virus (BTV), which infects sheep and cattle. It is an invasive species in Australia with an assumed Asian/South East Asian origin. Using one mitochondrial marker (i.e., part of the cytochrome oxidase subunit I gene) and six nuclear markers, we inferred population genetic structure and possible incursion pathways for Australian C. brevitarsis. Nine mitochondrial haplotypes, with low nucleotide sequence diversity (0.0-0.7%) among these, were identified in a sample of 70 individuals from seven sites. Both sets of markers revealed a homogeneous population structure, albeit with evidence of isolation by distance and two genetically distinct clusters distributed along a north-to-south cline. No evidence of a cryptic species complex was found. The geographical distribution of the mitochondrial haplotypes is consistent with at least two incursion pathways into Australia since the arrival of suitable livestock hosts. By contrast, 15 mitochondrial haplotypes, with up to four times greater nucleotide sequence diversity (0.0-2.9%) among these, were identified in a sample of 16 individuals of the endemic C. marksi (sampled from a site in South Australia and another in New South Wales). A phylogenetic tree inferred using the mitochondrial marker revealed that the Australian and Japanese samples of C. brevitarsis are as evolutionarily different from one another as some of the other Australian species (e.g., C. marksi, C. henryi, C. pallidothorax) are. The phylogenetic tree placed four of the species endemic to Australia (C. pallidothorax, C. bundyensis, C. marksi, C. henryi) in a clade, with a fifth such species (C. bunrooensis) sharing a common ancestor with that clade and a clade comprising two Japanese species (C. verbosus, C. kibunensis).

Changes in base composition bias of nuclear and mitochondrial genes in lice (Insecta: Psocodea)

While it is well known that changes in the general processes of molecular evolution have occurred on a variety of timescales, the mechanisms underlying these changes are less well understood. Parasitic lice ("Phthiraptera") and their close relatives (infraorder Nanopsocetae of the insect order Psocodea) are a group of insects well known for their unusual features of molecular evolution. We examined changes in base composition across parasitic lice and bark lice. We identified substantial differences in percent GC content between the clade comprising parasitic lice plus closely related bark lice (=Nanopsocetae) versus all other bark lice. These changes occurred for both nuclear and mitochondrial protein coding and ribosomal RNA genes, often in the same direction. To evaluate whether correlations in base composition change also occurred within lineages, we used phylogenetically controlled comparisons, and in this case few significant correlations were identified. Examining more constrained sites (first/second codon positions and rRNA) revealed that, in comparison to the other bark lice, the GC content of parasitic lice and close relatives tended towards 50 % either up from less than 50 % GC or down from greater than 50 % GC. In contrast, less constrained sites (third codon positions) in both nuclear and mitochondrial genes showed less of a consistent change of base composition in parasitic lice and very close relatives. We conclude that relaxed selection on this group of insects is a potential explanation of the change in base composition for both mitochondrial and nuclear genes, which could lead to nucleotide frequencies closer to random expectation (i.e., 50 % GC) in the absence of any mutation bias. Evidence suggests this relaxed selection arose once in the non-parasitic common ancestor of Phthiraptera + Nanopsocetae and is not directly related to the evolution of the parasitism in lice.

The complete nucleotide sequence of the mitochondrial genome of Bactrocera minax (Diptera: Tephritidae)

The complete 16,043 bp mitochondrial genome (mitogenome) of Bactrocera minax (Diptera: Tephritidae) has been sequenced. The genome encodes 37 genes usually found in insect mitogenomes. The mitogenome information for B. minax was compared to the homologous sequences of Bactrocera oleae, Bactrocera tryoni, Bactrocera philippinensis, Bactrocera carambolae, Bactrocera papayae, Bactrocera dorsalis, Bactrocera correcta, Bactrocera cucurbitae and Ceratitis capitata. The analysis indicated the structure and organization are typical of, and similar to, the nine closely related species mentioned above, although it contains the lowest genome-wide A+T content (67.3%). Four short intergenic spacers with a high degree of conservation among the nine tephritid species mentioned above and B. minax were observed, which also have clear counterparts in the control regions (CRs). Correlation analysis among these ten tephritid species revealed close positive correlation between the A+T content of zero-fold degenerate sites (P0FD), the ratio of nucleotide substitution frequency at P0FD sites to all degenerate sites (zero-fold degenerate sites, two-fold degenerate sites and four-fold degenerate sites) and amino acid sequence distance (ASD) were found. Further, significant positive correlation was observed between the A+T content of four-fold degenerate sites (P4FD) and the ratio of nucleotide substitution frequency at P4FD sites to all degenerate sites; however, we found significant negative correlation between ASD and the A+T content of P4FD, and the ratio of nucleotide substitution frequency at P4FD sites to all degenerate sites. A higher nucleotide substitution frequency at non-synonymous sites compared to synonymous sites was observed in nad4, the first time that has been observed in an insect mitogenome. A poly(T) stretch at the 5' end of the CR followed by a [TA(A)]n-like stretch was also found. In addition, a highly conserved G+A-rich sequence block was observed in front of the poly(T) stretch among the ten tephritid species and two tandem repeats were present in the CR.

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

The mitochondrial genomes of two nemerteans, Cephalothrix sp. (Nemertea: Palaeonemertea) and Paranemertes cf. peregrina (Nemertea: Hoplonemertea)

The mitochondrial genome sequences were determined for two species of nemerteans, Cephalothrix sp. (15,800 bp sequenced, near-complete) and Paranemertes cf. peregrina (14,558 bp, complete). As seen in most metazoans, the genomes encode 13 protein, 2 ribosomal RNA and 22 transfer RNA genes. The nucleotide composition is strongly biased toward A and T, as is typical for metazoan mtDNAs. There is also a significant strand skew in the distribution of these nucleotides, with the coding strand being richer in T than A and in G than C. All genes are transcribed in the same direction except for trnP and trnT, which is consistent with that reported for Cephalothrix hongkongiensis and Lineus viridis. Gene arrangement of Cephalothrix sp. is identical to that of C. hongkongiensis, while in P. cf. peregrina it is similar to L. viridis, but differs significantly from the three Cephalothrix species in the position of four protein-coding genes and seven tRNAs. Some protein-coding genes have 3' end stem-loop structures, which may allow mRNA processing without flanking tRNAs. The major non-coding regions observed in the two genomes with potential to form stem-loop structures may be involved in the initiation of replication or transcription. The average Ka/Ks values, varying from 0.12 to 0.89, are markedly different among the 13 mitochondrial protein-coding genes, suggesting that there may exist different selective pressure among mitochondrial genes of nemerteans.

Comparative mitogenomics of Braconidae (Insecta: Hymenoptera) and the phylogenetic utility of mitochondrial genomes with special reference to Holometabolous insects

BACKGROUND

Animal mitochondrial genomes are potential models for molecular evolution and markers for phylogenetic and population studies. Previous research has shown interesting features in hymenopteran mitochondrial genomes. Here, we conducted a comparative study of mitochondrial genomes of the family Braconidae, one of the largest families of Hymenoptera, and assessed the utility of mitochondrial genomic data for phylogenetic inference at three different hierarchical levels, i.e., Braconidae, Hymenoptera, and Holometabola.

RESULTS

Seven mitochondrial genomes from seven subfamilies of Braconidae were sequenced. Three of the four sequenced A+T-rich regions are shown to be inverted. Furthermore, all species showed reversal of strand asymmetry, suggesting that inversion of the A+T-rich region might be a synapomorphy of the Braconidae. Gene rearrangement events occurred in all braconid species, but gene rearrangement rates were not taxonomically correlated. Most rearranged genes were tRNAs, except those of Cotesia vestalis, in which 13 protein-coding genes and 14 tRNA genes changed positions or/and directions through three kinds of gene rearrangement events. Remote inversion is posited to be the result of two independent recombination events. Evolutionary rates were lower in species of the cyclostome group than those of noncyclostomes. Phylogenetic analyses based on complete mitochondrial genomes and secondary structure of rrnS supported a sister-group relationship between Aphidiinae and cyclostomes. Many well accepted relationships within Hymenoptera, such as paraphyly of Symphyta and Evaniomorpha, a sister-group relationship between Orussoidea and Apocrita, and monophyly of Proctotrupomorpha, Ichneumonoidea and Aculeata were robustly confirmed. New hypotheses, such as a sister-group relationship between Evanioidea and Aculeata, were generated. Among holometabolous insects, Hymenoptera was shown to be the sister to all other orders. Mecoptera was recovered as the sister-group of Diptera. Neuropterida (Neuroptera + Megaloptera), and a sister-group relationship with (Diptera + Mecoptera) were supported across all analyses.

CONCLUSIONS

Our comparative studies indicate that mitochondrial genomes are a useful phylogenetic tool at the ordinal level within Holometabola, at the superfamily within Hymenoptera and at the subfamily level within Braconidae. Variation at all of these hierarchical levels suggests that the utility of mitochondrial genomes is likely to be a valuable tool for systematics in other groups of arthropods.

Extreme primary and secondary protein structure variability in the chimeric male-transmitted cytochrome c oxidase subunit II protein in freshwater mussels: evidence for an elevated amino acid substitution rate in the face of domain-specific purifying selection

BACKGROUND

Freshwater unionoidean bivalves, and species representing two marine bivalve orders (Mytiloida and Veneroida), exhibit a mode of mtDNA inheritance involving distinct maternal (F) and paternal (M) transmission routes concomitant with highly divergent gender-associated mtDNA genomes. Additionally, male unionoidean bivalves have a approximately 550 bp 3' coding extension to the cox2 gene (Mcox2e), that is apparently absent from all other metazoan taxa.

RESULTS

Our molecular sequence analyses of MCOX2e indicate that both the primary and secondary structures of the MCOX2e region are evolving much faster than other regions of the F and M COX2-COX1 gene junction. The near N-terminus approximately 2/3 of the MCOX2e region contains an interspecifically variable number of predicted transmembrane helices (TMH) and interhelical loops (IHL) whereas the C-terminus approximately 1/3 is relatively conserved and hydrophilic while containing conserved functional motifs. MCOX2e displays an overall pattern of purifying selection that leads to the preservation of TMH/IHL and C-terminus tail sub-regions. However, 14 amino acid positions in the MCOX2e TMH/IHL sub-region might be targeted by diversifying selection, each representing a site where there exists interspecific variation for the constituent amino acids residing in a TMH or IHL.

CONCLUSION

Our results indicate that Mcox2e is unique to unionoidean bivalves, likely the result of a single insertion event that took place over 65 MYA and that MCOX2e is functional. The predicted TMH number, length and position variability likely stems from substitution-based processes rather than the typically implicated insertion/deletion events. MCOX2e has relatively high rates of primary and secondary structure evolution, with some amino acid residues potentially subjected to site-specific positive selection, yet an overall pattern of purifying selection leading to the preservation of the TMH/IHL and hydrophilic C-terminus tail subregions. The more conserved C-terminus tail (relative to the TMH/IHL sub-region of MCOX2e) is likely biologically active because it contains functional motifs. The rapid evolution of primary and secondary structure in MCOX2e, combined with the action of both positive and purifying selection, provide supporting evidence for the hypothesis that MCOX2e has a novel reproductive function within unionoidean bivalves. All tolled, our data indicate that unionoidean bivalve MCOX2 is the first reported chimeric animal mtDNA-encoded protein.

Genetic variability in the three genomes of Nasonia: nuclear, mitochondrial and Wolbachia

Nasonia consists of three closely related species of parasitoid wasps that are all infected with the endosymbiotic bacteria Wolbachia, a reproductive parasite common in arthropods. This situation presents the opportunity to compare patterns of variation in three associated genomes, Wolbachia and the nuclear and mitochondrial genomes of its host. Furthermore, although Nasonia wasps are emerging as a model for evolutionary and genetic studies, little is known about their genetic variability. Using amplified fragment length polymorphisms (AFLPs), all three species present a relatively high level of nuclear polymorphism and have different patterns of variation, with one of the species, Nasonia giraulti, being divided into two divergent subgroups. In each species, the mitochondrial pattern of variation is different from the nuclear pattern, possibly due to genetic hitchhiking of the mitochondria during (cytoplasmically inherited) Wolbachia sweeps. Mitochondria in Nasonia show a synonymous substitution rate approximately 10-15-fold higher than nuclear genes, probably reflecting an elevated mitochondrial mutation rate that is among the highest found in insects. Finally, all three species are doubly infected with their own strains of Wolbachia, one each from the two major supergroups (A and B). Sequence analysis reveals that each of the three Nasonia species acquired their A and B bacteria independently by horizontal transfer events from other insects with the exception of B type Wolbachia in N. longicornis and N. giraulti, which were acquired prior to speciation and then codiverged with the host. This represents one of the few clear-cut examples of codivergence of Wolbachia during host speciation.

GC-biased segregation of noncoding polymorphisms in Drosophila

The study of base composition evolution in Drosophila has been achieved mostly through the analysis of coding sequences. Third codon position GC content, however, is influenced by both neutral forces (e.g., mutation bias) and natural selection for codon usage optimization. In this article, large data sets of noncoding DNA sequence polymorphism in D. melanogaster and D. simulans were gathered from public databases to try to disentangle these two factors-noncoding sequences are not affected by selection for codon usage. Allele frequency analyses revealed an asymmetric pattern of AT vs. GC noncoding polymorphisms: AT --> GC mutations are less numerous, and tend to segregate at a higher frequency, than GC --> AT ones, especially at GC-rich loci. This is indicative of nonstationary evolution of base composition and/or of GC-biased allele transmission. Fitting population genetics models to the allele frequency spectra confirmed this result and favored the hypothesis of a biased transmission. These results, together with previous reports, suggest that GC-biased gene conversion has influenced base composition evolution in Drosophila and explain the correlation between intron and exon GC content.

Evolution of base-substitution gradients in primate mitochondrial genomes

Inferences of phylogenies and dates of divergence rely on accurate modeling of evolutionary processes; they may be confounded by variation in substitution rates among sites and changes in evolutionary processes over time. In vertebrate mitochondrial genomes, substitution rates are affected by a gradient along the genome of the time spent being single-stranded during replication, and different types of substitutions respond differently to this gradient. The gradient is controlled by biological factors including the rate of replication and functionality of repair mechanisms; little is known, however, about the consistency of the gradient over evolutionary time, or about how evolution of this gradient might affect phylogenetic analysis. Here, we evaluate the evolution of response to this gradient in complete primate mitochondrial genomes, focusing particularly on A-->G substitutions, which increase linearly with the gradient. We developed a methodology to evaluate the posterior probability densities of the response parameter space, and used likelihood ratio tests and mixture models with different numbers of classes to determine whether groups of genomes have evolved in a similar fashion. Substitution gradients usually evolve slowly in primates, but there have been at least two large evolutionary jumps: on the lineage leading to the great apes, and a convergent change on the lineage leading to baboons (Papio). There have also been possible convergences at deeper taxonomic levels, and different types of substitutions appear to evolve independently. The placements of the tarsier and the tree shrew within and in relation to primates may be incorrect because of convergence in these factors.

Phylogenetic analysis of the mitochondrial cytochrome oxidase subunit I gene of five species of the Culicoides imicola species complex

The phylogenetic status of members of the Culicoides imicola Kieffer (Diptera: Ceratopogonidae) species complex of haematophagous midges is unknown, and simple means to identify the members using all life stages are unavailable. In this study, the status of three confirmed (C. imicola s.s., C. bolitinos Meiswinkel and C. loxodontis Meiswinkel) and two provisional (C. tuttifrutti Meiswinkel and C. kwagga Meiswinkel) members of the complex from South Africa was assessed using phylogenetic analysis of partial DNA and amino acid sequences of the mitochondrial cytochrome oxidase subunit I (COI) gene. The four or five individuals of each species analysed contained one or two haplotypes each. Interspecific divergence was significant and characterized by strong A <--> T transversion bias. Phylogenetic trees constructed using neighbour-joining, parsimony and maximum likelihood showed each species to be distinct. Combinations of sites for two restriction enzymes in the COI sequences were species-specific and could form the basis of a diagnostic PCR assay.

The mitochondrial genome of the primary screwworm fly Cochliomyia hominivorax (Diptera: Calliphoridae)

The complete sequence of the mitochondrial genome of the screwworm Cochliomyia hominivorax was determined. This genome is 16,022 bp in size and corresponds to a typical Brachycera mtDNA. A Serine start codon for COI and incomplete termination codons for COII, NADH 5 and NADH 4 genes were described. The nucleotide composition of C. hominivorax mtDNA is 77% AT-rich, reflected in the predominance of AT-rich codons in protein-coding genes. Non-optimal codon usage was commonly observed in C. hominivorax mitochondrial genes. Phylogenetic analysis distributed the Acalypterate species as a monophyletic group and assembled the C. hominivorax (Calyptratae) and the Acalyptratae in a typical Brachycera cluster. The identification of diagnostic restriction sites on the sequenced mitochondrial genome and the correlation with previous RFLP analysis are discussed.

Evolution and structural organisation of mitochondrial DNA control region of myiasis-causing flies

This study reports the molecular characterization of the mtDNA control region (called the A+T-rich region in insects) of five dipteran species which cause myiasis: Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, Chrysomya megacephala Fabricius, Lucilia eximia Wiedemann (Diptera: Calliphoridae) and Dermatobia hominis Linnaeus Jr (Diptera: Oestridae). The control region in these species varies in length from 1000 to 1600 bp. Two structural domains with specific evolutionary patterns were identified. These were (1) conserved sequence blocks containing primary sequence motifs, including dinucleotide pyrimidine-purine series and long T-stretches, located at the 5' end adjacent to the tRNA(Ile) gene and (2) a hypervariable domain at the 3' end characterized by increased nucleotide divergence and size variation. A high frequency of A<-->T transversions at nucleotide substitution level indicated directional mutation pressure. The phylogenetic usefulness of the insect control region is discussed.

Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15-25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

The sequence and structure of the meadow grasshopper (Chorthippus parallelus) mitochondrial srRNA, ND2, COI, COII ATPase8 and 9 tRNA genes

The nucleotide sequence of the mitochondrial ND2, COI, COII, ATPase8, srRNA and nine tRNA genes have been sequenced from two individual of the meadow grasshopper Chorthippus parallelus. Comparisons are made to other insects for which the same regions are completely sequenced. Percentage A + T is found to be relatively low in C. parallelus though consistent with that of the other Orthopteran, Locusta migratoria. The relative number of substitutions observed in the different protein-coding genes was analysed between pairs of insect species sharing different levels of relatedness. A clear change in this rate was observed between the within-genus and between-genera comparisons. This change is interpreted in terms of the functional constraints acting on these four different genes. The patterns seem to result from an early saturation of COI and COII genes with synonymous substitutions, and a tolerance of ND2 and ATPase8 function to high levels of amino acid replacements. This analysis highlights a need for further sequence studies and comparisons between taxa of different levels of divergence in order to understand the patterns of mtDNA evolution on which many evolutionary investigations are based.

Accelerated evolution and Muller's rachet in endosymbiotic bacteria

Many bacteria live only within animal cells and infect hosts through cytoplasmic inheritance. These endosymbiotic lineages show distinctive population structure, with small population size and effectively no recombination. As a result, endosymbionts are expected to accumulate mildly deleterious mutations. If these constitute a substantial proportion of new mutations, endosymbionts will show (i) faster sequence evolution and (ii) a possible shift in base composition reflecting mutational bias. Analyses of 16S rDNA of five independently derived endosymbiont clades show, in every case, faster evolution in endosymbionts than in free-living relatives. For aphid endosymbionts (genus Buchnera), coding genes exhibit accelerated evolution and unusually low ratios of synonymous to nonsynonymous substitutions compared to ratios for the same genes for enterics. This concentration of the rate increase in nonsynonymous substitutions is expected under the hypothesis of increased fixation of deleterious mutations. Polypeptides for all Buchnera genes analyzed have accumulated amino acids with codon families rich in A+T, supporting the hypothesis that substitutions are deleterious in terms of polypeptide function. These observations are best explained as the result of Muller's ratchet within small asexual populations, combined with mutational bias. In light of this explanation, two observations reported earlier for Buchnera, the apparent loss of a repair gene and the overproduction of a chaperonin, may reflect compensatory evolution. An alternative hypothesis, involving selection on genomic base composition, is contradicted by the observation that the speedup is concentrated at nonsynonymous sites.

Unbiased estimation of symmetrical directional mutation pressure from protein-coding DNA

The most generally applicable procedure for obtaining estimates of the symmetrical, or strandnonspecific, directional mutation pressure (microD) on protein-coding DNA sequences is to determine the G+C content at synonymous codon sites (Psyn), and to divide Psyn by twice the arithmetic mean of the G+C content at synonymous codon sites of a large number of randomly generated, synonymously coding DNA sequences (Psyn). Unfortunately, the original procedure yields biased estimates of Psyn and microD and is computationally expensive. We here present a fast procedure for estimating unbiased microD values. The procedure employs direct calculation of Psyn (approximately Psyn) and two normalization procedures, one for Psyn < or = Psyn and another for Psyn > or = Psyn. The normalization removes a bias sometimes caused by codons specifying arginine, asparagine, isoleucine, and leucine. Consequently, comparison of protein-coding genes that are translated using different genetic codes is facilitated.