Molecular systematics of the bark lice infraorder Caeciliusetae (Insecta: Psocodea)

The phylogenetic relationships of bark lice and parasitic lice (Insecta: Psocodea) have been studied in a number of recent molecular phylogenetic analyses based on DNA sequences. Many of these studies have focused on the position of parasitic lice within the free‐living bark lice. However, fewer such studies have examined the relationships among major groups of free‐living bark lice and their implications for classification. In this study we focus on the infraorder Caeciliusetae, a large group of bark lice (?1000 species) within the suborder Psocomorpha. Using sequences of two mitochondrial and two nuclear genes, we estimated the phylogeny for relationships among the five recognized families within the infraorder Caeciliusetae. Based on the results, the sister‐group relationship and respective monophyly of Stenopsocidae and Dasydemellidae is strongly supported. Monophyly of the larger families Amphipsocidae and Caeciliusidae was not supported, although the causes of this were the placement of two distinct subfamilies (Paracaeciliinae and Calocaeciliinae). The monophyly of Asiopsocidae could not be tested because it was sampled only by one species. Based on these results and consideration of morphological characters, we propose a new classification for Caeciliusetae, recognizing six families: Amphipsocidae, Stenopsocidae, Dasydemellidae, Asiopsocidae, Paracaeciliidae and Caeciliusidae. We expect that this new classification will stabilize the higher‐level taxonomy of this group and help to identify groups in need of further work among these insects.     

Multiple origins of parasitism in lice

A major fraction of the diversity of insects is parasitic, as herbivores, parasitoids or vertebrate ectopara sites. Understanding this diversity requires information on the origin of parasitism in various insect groups. Parasitic lice (Phthiraptera) are the only major group of insects in which all members are permanent parasites of birds or mammals. Lice are classified into a single order but are thought to be closely related to, or derived from, book lice and bark lice (Psocoptera). Here, we use sequences of the nuclear 18S rDNA gene to investigate the relationships among Phthiraptera and Psocoptera and to identify the origins of parasitism in this group (termed Psocodea). Maximum-likelihood (ML), Bayesian ML and parsimony analyses of these data indicate that lice are embedded within the psocopteran infraorder Nanopsocetae, making the order Psocoptera paraphyletic (i.e. does not contain all descendants of a single common ancestor). Furthermore, one family of Psocoptera, Liposcelididae, is identified as the sister taxon to the louse suborder Amblycera, making parasitic lice (Phthiraptera) a polyphyletic order (i.e. descended from two separate ancestors). We infer from these results that parasitism of vertebrates arose twice independently within Psocodea, once in the common ancestor of Amblycera and once in the common ancestor of all other parasitic lice.     

Dramatically elevated rate of mitochondrial substitution in lice (Insecta: Phthiraptera)

Few estimates of relative substitution rates, and the underlying mutation rates, exist between mitochondrial and nuclear genes in insects. Previous estimates for insects indicate a 2-9 times faster substitution rate in mitochondrial genes relative to nuclear genes. Here we use novel methods for estimating relative rates of substitution, which incorporate multiple substitutions, and apply these methods to a group of insects (lice, Order: Phthiraptera). First, we use a modification of copath analysis (branch length regression) to construct independent comparisons of rates, consisting of each branch in a phylogenetic tree. The branch length comparisons use maximum likelihood models to correct for multiple substitution. In addition, we estimate codon-specific rates under maximum likelihood for the different genes and compare these values. Estimates of the relative synonymous substitution rates between a mitochondrial (COI) and nuclear (EF-1alpha) gene in lice indicate a relative rate of several 100 to 1. This rapid relative mitochondrial rate (>100 times) is at least an order of magnitude faster than previous estimates for any group of organisms. Comparisons using the same methods for another group of insects (aphids) reveals that this extreme relative rate estimate is not simply attributable to the methods we used, because estimates from aphids are substantially lower. Taxon sampling affects the relative rate estimate, with comparisons involving more closely related taxa resulting in a higher estimate. Relative rate estimates also increase with model complexity, indicating that methods accounting for more multiple substitution estimate higher relative rates.     

Flight loss linked to faster molecular evolution in insects

The loss of flight ability has occurred thousands of times independently during insect evolution. Flight loss may be linked to higher molecular evolutionary rates because of reductions in effective population sizes (N_e) and relaxed selective constraints. Reduced dispersal ability increases population subdivision, may decrease geographical range size and increases (sub)population extinction risk, thus leading to an expected reduction in N_e. Additionally, flight loss in birds has been linked to higher molecular rates of energy-related genes, probably owing to relaxed selective constraints on energy metabolism. We tested for an association between insect flight loss and molecular rates through comparative analysis in 49 phylogenetically independent transitions spanning multiple taxa, including moths, flies, beetles, mayflies, stick insects, stoneflies, scorpionflies and caddisflies, using available nuclear and mitochondrial protein-coding DNA sequences. We estimated the rate of molecular evolution of flightless (FL) and related flight-capable lineages by ratios of non-synonymous-to-synonymous substitutions (dN/dS) and overall substitution rates (OSRs). Across multiple instances of flight loss, we show a significant pattern of higher dN/dS ratios and OSRs in FL lineages in mitochondrial but not nuclear genes. These patterns may be explained by relaxed selective constraints in FL ectotherms relating to energy metabolism, possibly in combination with reduced N_e.     

Systematic identification and characterization of stress-inducible heat shock proteins (HSPs) in the salmon louse (<Emphasis Type="Italic">Lepeophtheirus salmonis</Emphasis>)

Salmon lice (Lepeophtheirus salmonis) are parasitic copepods, living mainly on Atlantic salmon and leading to large economical losses in aquaculture every year. Due to the emergence of resistances to several drugs, alternative treatments are developed, including treatment with hydrogen peroxide, freshwater or thermal treatment. The present study gives a first overview of the thermotolerance and stress response of salmon lice. Sea lice nauplii acclimated to 10 °C can survive heat shocks up to 30 °C and are capable of hardening by a sublethal heat shock. We searched in the genome for heat shock protein (HSP) encoding genes and tested their inducibility after heat shock, changes in salinity and treatment with hydrogen peroxide, employing microfluidic qPCRs. We assessed 38 candidate genes, belonging to the small HSP, HSP40, HSP70 and HSP90 families. Nine of these genes showed strong induction after a non-lethal heat shock. In contrast, only three and two of these genes were induced after changes in salinity and incubation in hydrogen peroxide, respectively. This work provides the basis for further work on the stress response on the economically important parasite L. salmonis.     

Gene trees,species trees and Earth history combine to shed light on the evolution of migration in a model avian system

The evolution of migration in birds has fascinated biologists for centuries. In this study, we performed phylogenetic-based analyses of Catharus thrushes, a model genus in the study of avian migration, and their close relatives. For these analyses, we used both mitochondrial and nuclear genes, and the resulting phylogenies were used to trace migratory traits and biogeographic patterns. Our results provide the first robust assessment of relationships within Catharus and relatives and indicate that both mitochondrial and autosomal genes contribute to overall support of the phylogeny. Measures of phylogenetic informativeness indicated that mitochondrial genes provided more signal within Catharus than did nuclear genes, whereas nuclear loci provided more signal for relationships between Catharus and close relatives than did mitochondrial genes. Insertion and deletion events also contributed important support across the phylogeny. Across all taxa included in the study, and for Catharus, possession of long-distance migration is reconstructed as the ancestral condition, and a North American (north of Mexico) ancestral area is inferred. Within Catharus, sedentary behaviour evolved after the first speciation event in the genus and is geographically and temporally correlated with Central American distributions and the final closure of the Central American Seaway. Migratory behaviour subsequently evolved twice in Catharus and is geographically and temporally correlated with a recolonization of North America in the late Pleistocene. By temporally linking speciation events with changes in migratory condition and events in Earth history, we are able to show support for several competing hypotheses relating to the geographic origin of migration.     

Conservative and compensatory evolution in oxidative phosphorylation complexes of angiosperms with highly divergent rates of mitochondrial genome evolution

Interactions between nuclear and mitochondrial gene products are critical for eukaryotic cell function. Nuclear genes encoding mitochondrial‐targeted proteins (N‐mt genes) experience elevated rates of evolution, which has often been interpreted as evidence of nuclear compensation in response to elevated mitochondrial mutation rates. However, N‐mt genes may be under relaxed functional constraints, which could also explain observed increases in their evolutionary rate. To disentangle these hypotheses, we examined patterns of sequence and structural evolution in nuclear‐ and mitochondrial‐encoded oxidative phosphorylation proteins from species in the angiosperm genus Silene with vastly different mitochondrial mutation rates. We found correlated increases in N‐mt gene evolution in species with fast‐evolving mitochondrial DNA. Structural modeling revealed an overrepresentation of N‐mt substitutions at positions that directly contact mutated residues in mitochondrial‐encoded proteins, despite overall patterns of conservative structural evolution. These findings support the hypothesis that selection for compensatory changes in response to mitochondrial mutations contributes to the elevated rate of evolution in N‐mt genes. We discuss these results in light of theories implicating mitochondrial mutation rates and mitonuclear coevolution as drivers of speciation and suggest comparative and experimental approaches that could take advantage of heterogeneity in rates of mtDNA evolution across eukaryotes to evaluate such theories.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻(Oryza rufipogon Griff.)线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻(O. sativa L.)的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示:普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数(Nc)分布于45.32～61.00之间,其密码子偏性较弱; Nc值仅与GC_3呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数(CBI)和最优密码子使用频率(Fop)的相关性均达到显著水平;而GC_3和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度(RSCU)、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻(O. sativa L. subsp. japonica Kato)和籼稻(O. sativa L. subsp.indica Kato)内具有同质性。     

Morphology of male genitalia in lice and their relatives and phylogenetic implications

Abstract.  Lice (Insecta: Phthiraptera) have long been considered to compose a monophyletic group of insects on the basis of external morphological characteristics. However, a recent phylogenetic analysis of 18S rDNA sequences suggested that 'Phthiraptera' have arisen twice within the order Psocoptera (booklice and barklice). The external features of lice are highly specialized to a parasitic lifestyle, and convergence may be frequent for such characters. To provide a further test between traditional and recent molecular-based phylogenetic hypotheses, a phylogenetic analysis of lice and relatives based on morphological characters that are independent from the selective pressures of a parasitic lifestyle is needed. Here, we examined the morphology of the male phallic organ in lice and relatives ('Psocoptera': suborders Troctomorpha and Psocomorpha) and detected some novel modifications that were stable within each group and useful for higher level phylogenetic reconstruction. Phylogenetic analysis based on these characters provided a concordant result with the 18S-based phylogeny. In particular, the apomorphic presence of articulations between the basal plate, mesomere and ventral plate (= sclerite on the permanently everted endophallus) is observed consistently throughout the psocid families Pachytroctidae and Liposcelididae and the louse suborder Amblycera, providing support for a clade composed of these three groups, although possible homoplasy was detected in some Ischnocera. This is the first study to provide morphological support for the polyphyly of lice.     

Factors affecting mito-nuclear codon usage interactions in the OXPHOS system of Drosophila melanogaster

Codon usage bias varies considerably among genomes and even within the genes of the same genome.In eukaryotic organisms,energy production in the form of oxidative phosphorylation(OXPHOS)is the only process under control of both nuclear and mitochondrial genomes.Although factors affecting codon usage in a single genome have been studied,this has not occurred when both interactional genomes are involved.Consequently, we investigated whether or not other factors influence codon usage of coevolved genes.We used Drosophila melanogaster as a model organism.Our χ2 test on the number of codons of nuclear and mitochondrial genes involved in the OXPHOS system was significantly different (χ2=7945.16,P<0.01).A plot of effective number of codons against GC3s content of nuclear genes showed that few genes lie on the expected curve,indicating that codon usage was random.Correspondence analysis indicated a significant correlation between axis 1 and codon adaptation index(R=0.947,P<0.01)in every nuclear gene sequence.Thus,codon usage bias of nuclear genes appeared to be affected by translational selection.Correlation between axis 1 coordinates and GC content(R=0.814.P<0.01)indicated that the codon usage of nuclear genes was also affected by GC composition.Analysis of mitochondrial genes did not reveal a significant correlation between axis 1 and any parameter.Statistical analyses indicated that codon usages of both nDNA and mtDNA were subjected to context-dependent mutations.     

The geographical segregation of human lice preceded that of Pediculus humanus capitis and Pediculus humanus humanus

In order to investigate human-louse phylogeny, we partially sequenced two nuclear (18S rRNA and EF-1 alpha) and one mitochondrial (COI) genes from 155 Pediculus from different geographical origins. The phylogenetic analysis of 18S rRNA and EF-1 alpha sequences showed that human lice were classified into lice from Sub-Saharan Africa and lice from other areas. In both clusters, head and body lice were clearly grouped into two separate clusters. Our results indicate that the earliest divergence within human pediculidae occurred between African lice and other lice, and the divergence between head and body lice was not the result from a single event.     

Evolutionary dynamics of wheat mitochondrial gene structure with special remarks on the origin and effects of RNA editing in cereals

We investigated the evolutionary dynamics of wheat mitochondrial genes with respect to their structural differentiation during organellar evolution, and to mutations that occurred during cereal evolution. First, we compared the nucleotide sequences of three wheat mitochondrial genes to those of wheat chloroplast, alpha-proteobacterium and cyanobacterium orthologs. As a result, we were able to (1) differentiate the conserved and variable segments of the orthologs, (2) reveal the functional importance of the conserved segments, and (3) provide a corroborative support for the alpha-proteobacterial and cyanobacterial origins of those mitochondrial and chloroplast genes, respectively. Second, we compared the nucleotide sequences of wheat mitochondrial genes to those of rice and maize to determine the types and frequencies of base changes and indels occurred in cereal evolution. Our analyses showed that both the evolutionary speed, in terms of number of base substitutions per site, and the transition/transversion ratio of the cereal mitochondrial genes were less than two-fifths of those of the chloroplast genes. Eight mitochondrial gene groups differed in their evolutionary variability, RNA and Complex I (nad) genes being most stable whereas Complex V (atp) and ribosomal protein genes most variable. C-to-T transition was the most frequent type of base change; C-to-G and G-to-C transversions occurred at lower rates than all other changes. The excess of C-to-T transitions was attributed to C-to-U RNA editing that developed in early stage of vascular plant evolution. On the contrary, the editing of C residues at cereal T-to-C transition sites developed mostly during cereal divergence. Most indels were associated with short direct repeats, suggesting intra- and intermolecular recombination as an important mechanism for their origin. Most of the repeats associated with indels were di- or trinucleotides, although no preference was noticed for their sequences. The maize mt genome was characterized by a high incidence of indels, comparing to the wheat and rice mt genomes.     

The human genome structure and organization

Genetic information of human is encoded in two genomes: nuclear and mitochondrial. Both of them reflect molecular evolution of human starting from the beginning of life (about 4.5 billion years ago) until the origin of Homo sapiens species about 100,000 years ago. From this reason human genome contains some features that are common for different groups of organisms and some features that are unique for Homo sapiens. 3.2 x 10(9) base pairs of human nuclear genome are packed into 23 chromosomes of different size. The smallest chromosome - 21st contains 5 x 10(7) base pairs while the biggest one -1st contains 2.63 x 10(8) base pairs. Despite the fact that the nucleotide sequence of all chromosomes is established, the organisation of nuclear genome put still questions: for example: the exact number of genes encoded by the human genome is still unknown giving estimations from 30 to 150 thousand genes. Coding sequences represent a few percent of human nuclear genome. The majority of the genome is represented by repetitiVe sequences (about 50%) and noncoding unique sequences. This part of the genome is frequently wrongly called "junk DNA". The distribution of genes on chromosomes is irregular, DNA fragments containing low percentage of GC pairs code lower number of genes than the fragments of high percentage of GC pairs.     

Mode and Tempo of Molecular Evolution in the Nematode Caenorhabditis: Cytochrome Oxidase II and Calmodulin Sequences

Through direct sequencing methods, the mitochondrial gene for cytochrome oxidase subunit two (CO II) and the single-copy nuclear gene for calmodulin were compared among strains of Caenorhabidits elegans and two other Caenorhabditis species (C. remanei and C. briggsae). In addition the CO II sequence was determined from a distantly related nematode, Steinernema intermedii. Among the 11 strains of C. elegans tested, there are four types of CO II gene, arising from two major lineages. Levels of intraspecific difference in the CO II gene are low (less than 2.0%) compared to the extraordinary divergence between congeneric species, which is about 50% when corrected for multiple hits. Concordant with the increase in divergence between taxa is a change in the pattern of substitution from a strong transition bias (24 transitions compared to two transversions) within species to a substitution pattern that appears to reflect the base composition of the mitochondrial genome when more divergent nematodes are compared. The base composition of the Caenorhabditis CO II gene is strongly biased toward A + T at all three positions of codons and appears to constrain the amino acid composition of the protein. Both the CO II and calmodulin genes show extreme conservation of amino acid sequences. When the accumulation of changes at silent sites in the two genes is compared among strains, it becomes evident that the mitochondrial gene is changing faster than the nuclear gene.     

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Silent nucleotide substitutions and G+C content of some mitochondrial and bacterial genes

Summary The G+C content of DNA varies widely in different organisms, especially microorganisms. This variation is accompanied by changes in the nucleotide composition of silent positions in codons. (Silent positions are defined and explained in the text.) These changes are mostly neutral or near neutral, and appear to result from mutation pressure in the direction of increasing either A+T (AT pressure) or G+C(GC pressure) content. Variations in G+C content are also accompanied by substitutions at replacement positions in codons. These substituions produce changes in the amino acid content of homologous proteins. The examples studied were genes for 13 mitochondrial proteins in five species, and A and B genes for bacterial tryptophan synthase in four species.In microorganisms, varying AT and GC mutational pressures, presumably resulting from shifts in the DNA polymerase system, exert strong effects on molecular evolution by changing the G+C content of DNA. These effects may be greater than those of random drift. The effects of GC pressure on silent substitutions in the systems examined are several times as great as the effects on replacement substitutions.GC pressure is exerted on noncoding as well as coding regions in mitochondrial DNA. This is shown by the close correlation (correlation coefficient, 0.99) of the G+C content of the noncoding D loop of mitochondria with the G+C content of silent positions in the corresponding mitochondrial genes.     

A comprehensive analysis of mammalian mitochondrial genome base composition and improved phylogenetic methods

Phylogenetic analysis of mammalian species using mitochondrial protein genes has proved to be problematic in many previous studies. The high mutation rate of mitochondrial DNA and unusual base composition of several species has prompted us to conduct a detailed study of the composition of 69 mammalian mitochondrial genomes. Most major changes in base composition between lineages can be attributed to shifts between the proportions of C and T on the L-strand. These changes are significant at all codon positions and are shown to affect amino acid composition. Correlated changes in the base composition of the RNA loops and stems are also observed. Following up from previous studies, we investigate changes in the base composition of all 12 H-strand proteins and find that variability in proportions of C and T is correlated with location on the genome. Variation in base composition across genes and species is known to adversely affect the performance of phylogenetic inference methods. We have, therefore, developed a customized three-state general time-reversible DNA substitution model, implemented in the PHASE phylogenetic inference package, which lumps C and T into a composite pyrimidine state. We compare the phylogenetic tree obtained using the new three-state model with that obtained using a standard four-state model. Results using the three-state model are more congruent with recent studies using large sets of nuclear genes and help resolve some of the apparent conflicts between studies using nuclear and mitochondrial proteins.     

Cophylogenetic relationships between penguins and their chewing lice

It is generally thought that the evolution of obligate parasites should be linked intimately to the evolution of their hosts and that speciation by the hosts should cause speciation of their parasites. The penguins and their chewing lice present a rare opportunity to examine codivergence between a complete host order and its parasitic lice. We estimated a phylogeny for all 15 species of lice parasitising all 17 species of penguins from the third domain of the mitochondrial 12S ribosomal rRNA gene, a portion of the mitochondrial cytochrome oxidase subunit 1 gene and 55 morphological characters. We found no evidence of extensive cospeciation between penguins and their chewing lice using TreeMap 2.02beta. Despite the paucity of cospeciation, there is support for significant congruence between the louse and penguin phylogenies due to possible failure to speciate events (parasites not speciating in response to their hosts speciating).     

Characterisation of the complete mitochondrial genome of Helice wuana (Grapsoidea: Varunidae) and comparison with other Brachyuran crabs

The mitochondrial genome (mitogenome) provides important information for phylogenetic analysis and understanding evolutionary origins. Herein, we sequenced, annotated, and characterised the mitogenome of the crab Helice wuana to better understand its molecular evolution and phylogeny. The 16,359 bp mitogenome includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region. The genome composition is highly A + T biased 68.42%, and exhibits a negative AT–skew (? 0.036) and GC–skew (? 0.269) among Brachyura crabs. Gene rearrangements were detected, as was tandem duplication followed by random loss, which explains the translocation of mitochondrial genes. Phylogenetic analysis showed that H. wuana and H. tientsinensis clustered on one branch with high nodal support values. These results confirm that the placement of H. wuana within the Varunidae family of Thoracotrematan crabs. This study will provided a better understanding for gene rearrangements and crab evolution in the further.